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HAUPTAUFSÄTZE 


Strömung in gekrümmten Rohren’. 
Von M. Adler in Wembley (England). 


(Aus dem Kaiser Wilhelm-Institut für Strömungsforschung, Göttingen.) 


D: Hauptmerkmale der Strömung durch gekrümmte Rohre sind die Veränderung der 
Strömungsform und die Vergrößerung des Widerstandes gegenüber der geraden Rohr- 
strömung. Hinsichtlich der Strömungsform ist für den Krümmer das Auftreten einer 
Sekundärströmung charakteristisch, die sich 
der Hauptströmung überlagert. Ihre Ausbildung 
wird dadurch erklärt, daß sich die in der 
Mitte des Rohres mit größerer Geschwindigkeit 
fließenden Flüssigkeitsteile durch ihre größere 
Zentrifugalkraft an die Außenwand drängen 
und dabei die dortigen langsameren Partien ın 
zwei Teilströmen über oben und über unten 
gegen die Innenwand schieben. Auf diese 
Weise erhält die Sekundärströmung die in Abb. 1 
angedeutete Form eines Doppelwirbels. Durch 
Überlagerung der Sekundärströmung mit der 
Hauptströmung entsteht eine schraubenähnliche 
Bewegung der einzelnen Flüssigkeitsteile. Die 
Sekundärströmung hat zur Folge, daß die Stelle 
größter Geschwindigkeit von der Mitte nach 
außen verschoben wird und ist die Ursache der 
Widerstandsvermehrung in gekrümmten Rohren. 

Beide Erscheinungen, Veränderung der 
Strömungsform und Widerstandserhöhung, waren 
bisher Gegenstand vieler experimenteller Untersuchungen. Das Auftreten der Sekundär- 
strömung haben Eustice') und Taylor?) durch Einführung von Farbfäden in die Strömung 
sichtbar gemacht. Versuche über die Geschwindigkeitsverteilung sind von Saph und 

















*) Dissertation Technische Hochschule München. 1. Berichterstatter: Prof. Dr.-Ing. D. Thoma, 2. Berichterstatter: 
Prof. Dr.-Ing. W. Kaufmann. Die Anregung zu dieser Arbeit erhielt ieh von Herrn Prof. Dr. L. Prandt]|, dem ich 
für die stets wohlwollende Unterstützung und die wertvollen Ratschläge bei der Durehführung meinen ergebensten 
Dank ausspreche. 

1) J. Eustice: Flow of Fluids in Curved Passages. Engineering Bd. 120 (1925), S. 604. 

?)G. J. Taylor: The Criterion for Turbulence in Curved Passages. Proe. Roy. Soe., Lond. (A) Bd. 124 
(1929), 8.243. 
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Schoder°’) ausgeführt worden. Wegen der technischen Wichtigkeit von Rohrkrümmungen 
ist es verständlich, daß vor allem der Einfluß der Sekundärströmung auf den Widerstand 
und die technisch besonders stark interessierenden Widerstandszahlen, hauptsächlich in 
turbulenten Bereich, eingehend experimentell untersucht worden sind. An neueren Arbeiten 
seien erwähnt die Forschungsarbeit von Nippert*), in der die Verlustbildung in Krümmern 
mit reehteekigem und kreisförmigem Quersehnitt in Abhängigkeit von den sie beeinflussenden 
Faktoren untersucht worden ist, sowie die Arbeit von Richter?°). der auch die Laminar- 
strömung behandelt, und die im Hydraulischen Institut der Technischen Hochschule München 
ausgeführten Versuche®). 

Demgegenüber ist die Zahl der theoretischen Untersuchungen über die Strömung in 
eekrümmten Rohren noch sehr gering. Von Bedeutung ist in dieser Hinsicht die Arbeit von 


Dean’) der unter Annahme kleinen Krümmungsverhältnisses (a = Rohrradius, AR — mitt- 


> 
Li 


lerer Krümmungsradius) in sukzessiver Näherung, die bis zum 4. Schritt getrieben wird, aus 

den Navier-Stokesschen Bewegungsgleichungen eine Widerstandsformel errechnet, die 
(7 a‘ , ' ; a a 

von dem Parameter 2 2 7 (® Druckgefälle, «= Zähigkeit, » = kinematische Zähigkeit) ab- 
ev A ; 

hängie ist. Die errechnete Funktion ist aber nur im laminaren Bereich bei sehr kleinen 

Revnoldsschen Zahlen brauchbar. Jedoch ermöglicht es die Kenntnis des Parameters, auch 

außerhalb des kleinen Bereiches experimentelle Ergebnisse auf die gefundene Abhängigkeit 

nachzuprüfen. Dieser Aufgabe unterzieht sich White°), der den Deanschen Parameter auf 

p® a\s ed z are 

die Form Re| „) Re Revnoldssche Zahl, e= mittlere Geschwindigkeit, d = Rohr- 

v/ r P 


durchmesser) umreehnet. 


Aufgabe dieser Arbeit soll es sein, auf Grund der Prandtlschen Grenzschichtbetrachtung 
eine Theorie zu entwickeln, die für laminare Strömung, aber für größere Reynoldssche 
Zahlen bis zum kritischen Gebiet brauchbar ist. Außerdem sollen Widerstand und Ge- 
schwindigkeitsverteilung, im laminaren und im turbulenten Bereich, experimentell unter- 
sucht werden. 

I. Theoretischer Teil. 


1. Theorie der laminaren Strömung für kleine Reynoldssche Zahlen und kleine Krümmungen. 
/u dem von Dean erzielten Ergebnis kann man auch auf einem anderen Weg gelangen, der 
hier kurz angedeutet sein soll. Man geht dabei aus von den auf Zylinderkoordinaten um- 
verechneten Navier-Stokesschen Bewegungsgleichungen. « und ® seien die Gescehwindig- 
keiten in Richtung der reehtwinkligen Koordinaten = und y eines Rohrquerschnittes. der mit 
dem Nullquerschnitt den Winkel g um die Achse durch den Krümmungsmittelpunkt bildet. 
ı sei die Geschwindigkeit in Strömungsrichtung. Der Krümmungsradius AR des Rohres sei 


i . Bi ” . 4a } u 
oeroß gegen den Radius a des Rohrquerschnittes. Das Krümmungsverhältnis R wird infolge- 
[) 
dessen klein und sei von der Größenordnung e. 
Für die Geschwindigkeiten #, », w schreiben wir jetzt die Reihen: 
" “„.+ Wu,+... v=v, tv,+... v—w, tw, W,-+... 
, . R () 
j E&” € y | € € 
® ® b) .p ® m> - ‘ .p A 
an und vergleichen die Größenordnung ihrer Glieder mit der Größenordnung & von PB: Das 
1 


Glied »*, der Hauptgeschwindigkeit mw ist gleich der Geschwindigkeit, die vorhanden wäre, 
wenn das Rohr keine Krümmung hätte und möge von der Größenordnung 1 sein. Die 
entsprechenden Glieder von « und ® sind Null, da im geraden Rohr keine Sekundär- 
strömung auftritt. Der Einfluß der Krümmung macht sich erst von w, ab geltend. Dieses 


Glied ist von der Größenordnung der Krümmung Ri ebenso wie die Glieder «, und ®, 
i 


der Sekundärströmung. Die folgenden Glieder werden immer kleiner, von der Größenordnung 


} 


e, € usw. Die Reihenansätze (1) führen wir in die Bewegungsgleichungen ein. Berück- 


3), A. V.Saph und E. W. Scehoder: Diseussion on Flow of Water in Pipes. Trans. Amer. Soc. Civ. Eng. 
Bd. 47 (1902), S. 299. 

3, 11. Nippert: Über den Strömungswiderstand in gekrümmten Kanälen. Forseh.-Arb. Ing.-Wes. Nr. 320 (1929). 

5), H. Riehter: Der Druckabfall in gekrümmten glatten Rohrleitungen. Forsch.-Arb. Ing.-Wes. Nr. 338 (1930). 

DPD. Thoma: Der hydraulische Verlust in Formstüecken. World Power Conference, Seet. Meeting, Nr. 26, 
Tokyo 1929, 

', W.R. Dean: The Streamline Motion of Fluid in a Curved Pipe. Phil. Mag. (7). Bd. 4 (1927). S. 208 und 
Bd. > (1928), S. 673. 

‘)C.M. White: Steamline Flow trough Curved Pipes. Proe. Roy. Sce., Lond. (A) Bd. 123 (1929), 8. 645. 
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siehtigen wir dann zunächst nur die Glieder von der Größenordnung 1, so erhalten wir aus 
der dritten Bewegungsgleichung eine Differentialgleichung mit der einzigen Unbekannten ı,. 
Die Lösung dieser Differentialgleichung ist die Poiseuillesche Gleichung für die gerade 
Rohrströmung. Dann gehen wir einen Schritt weiter. Zur Berechnung der ersten Näherung 
der Sekundärströmung berücksichtigen wir nur die Glieder von der Größenordnung &e und 
erhalten aus den beiden ersten Bewegungsgleiehungen nach Elimination des Druckgliedes eine 
Differentialgleichung mit den Unbekannten «, und », und der jetzt schon bekannten 
Funktion ,. Eine der beiden Unbekannten läßt sich aus dieser Differentialgleichung be- 
stimmen, die andere ergibt sich aus der Kontinuitätsgleichung, für die natürlich auch nur 
die «Glieder anzuschreiben sind. Mit den nun bekannten Werten »*,, «, und v, kann man 
jetzt w, aus der dritten Bewegungsgleichung unter Vernachlässigung aller von der Größen- 
ordnung e verschiedenen Glieder berechnen. Durch Fortsetzung dieses Verfahrens kann man 
immer bessere Annäherungen erzielen. Diese Aufgabe scheitert jedoch sehr bald an der 
immer schwieriger werdenden Lösung der Differentialgleichungen. 

Aus der Geschwindigkeitsverteilung in Strömungsriehtung kann durch Integration über 
den Rohrquersehnitt die Durchflußmenge @ berechnet werden. PBezieht man sie auf die 
Durehflußmenge 9, des geraden Rohres, die bei dem gleichen Druckgefälle zustande kommt, 


h () . 2 ar 
so kann man den (Juotienten V als Maß für den Verlust ansehen, der durch die Krümmung 


0 


verursacht wird. Bei dem erreichten Stand der Näherung ist mit 


| Q P\ a‘ 


K 5 (>) 
24 ROy/Sw#rR 
Op a un - h ci 
(23 Druckgefälle in Strömungsrichtung, bezogen auf die Längeneinheit, «= Zähigkeit, 
04 
» - kinematische Zähigkeit): 

() en u : 
1008ER H+0HMBER! 9) 2 22. 


Über den Gültiekeitsbereich dieser Formel läßt sich folgendes aussagen: Der Verlust ist um 
so größer, je stärker bei gleichbleibender Durchflußgeschwindigkeit das Rohr gekrümmt ist, 


i I u, 
je größer also das Krümmungsverhältnis IE und damit auch X, ist. Mit wachsendem A 
D Ä Ö N 


() ’ u () a z 
muß also ) kleiner werden. Das Minimum von ER errechnet aus Gl. (3). liegt angenähert 


“0 -ov 
bei X 1,12: von da ab beginnt sein Wert mit wachsendem A wieder zu steigen, d. h., beı 
K—1,12 ist die Lösung bei dem erreichten Stand der Näherung sicherlich nicht mehr ver- 
wendbar. Für Wasser von 15°C, das durch r 
ein schwach gekrümmtes Rohr mit lem 8 und u) 
dem Krümmungsverhältnis - oo R=S m) 
fließt, entspricht diesem K-Wert ein Druck- 
vefälle von 3 - 10° * ke/em? oder 3 mm W.-S. für 
Im Rohrlänge und eine Reynoldssche Zahl 
von Re=700. Für eine stärkere Krümmung 
Ad | i , Br r m | 
RT (R = 2m) dürfte das Wasser an der 
Grenze des Gültigkeitsbereiches der Näherung 
nur mit einem Druckgefälle von 1,5: 10° * kg/em’ 
und einer Reynoldsschen Zahl von 350 
durch das Rohr fließen. 

Beim Geschwindigkeitsprofil ist inner- 
halb dieses Bereiches keine wesentliche Ab- 








weichung von dem Profil des geraden Rohres BEER ERERE VEONER WERNE: — ul 11 
vorhanden. In Abb. 2 ist ein Geschwindig- am Yun a 5 
keitsprofil für einen speziellen Fall Abb. 2. 


Op a | 
5 .10=#% kolem? 
we kgiem‘, 5 gn 


in seinen einzelnen Näherungen aufgezeichnet. 


R 4m r=0, 01145 em?/s 


9, Als Zahlenfaktor von AK? erhält Dean 0,03058. Der Untersehied erklärt sieh dureh einige Vernachlässienngen in 
der Deanschen Reehnung. 


14° 
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Man sieht also, daß vorliegende Theorie nur für ein begrenztes Gebiet kleiner Reynolds- 
scher Zahlen und kleiner Krümmungen brauchbar ist. Die Abhängigkeit vom Parameter K 
ist natürlich auch im weiteren Bereich der laminaren Strömung vorhanden, nur läßt sich 
keine Aussage darüber machen, durch welche Funktion er dargestellt ist. Wohl aber ist die 
Kenntnis des Parameters dazu geeignet, experimentelle Ergebnisse nachzuprüfen, indem man 
” über A aufträgt. Die Versuchspunkte, die verschiedenen Versuchsbedingungen, vor allem 


( 


auch verschiedenen Krümmungen entsprechen, müssen sich auf einer einzigen Kurve anordnen, 
falls die Gültigkeit des Parameters zutrifft. White!") hat eine solche Prüfung vorgenommen 
und dabei eine gute Bestätigung gefunden. Er gestaltet lediglich in formaler Weise die Be- 


() * * Y . * 

zichung ,, ft) oder D f,(K) um. Für die linke Seite kann man schreiben: 
7 c, a” O p | | ed OÖ p d / 
() C Su Roy e 64 v Ro * 4, 
i O0 

„3 
ua Op dd ai u u z 
Dabei ist 4 AR -„ definitionsgemäß die Widerstandszahl des gekrümmten Rohres (Druck- 
ı0Q C E 
Sar57 


u ö 64 .. 
oefälle bezogen auf den Staudruck der mittleren Durchtlußgeschwindigkeit) und 4, = = ] die 
u Cl 
Y 
Widerstandszahl des geraden Rohres bei der gleichen mittleren Durchflußgeschwindigkeit. 








Il op\ «a I Il ed/a\s]? 
K= Ru läßt sich umformen zu K = — |- (5) ‚so daß man schließlich 
Er Ro q | Pe) u” v" R 2 112 R 
| BE. 2. 2 B 
mit Re = die Beziehung erhält: 
5 
/ 1] 


| a 
[Re ;) 
I | fi 
Auf diese Beziehung wird im Zusammenhang mit den experimentellen Ergebnissen noch 
einzugehen sein. 


2. Theorie der laminaren Strömung für große Reynoldssche Zahlen. Im folgenden soll 
nun der Versuch gemacht werden, eine Darstellung zu finden, die für große Reynoldssche 
Zahlen im laminaren Bereich brauchbar ist. Den Grundgedanken für eine solche Lösung 
bietet die Prandtlsche Grenzschichttheorie''). Danach hat man nur innerhalb einer sehr 
dünnen Schicht an der Begrenzung der Flüssigkeit, der Grenzschicht, in der wegen des 
Haftens der Flüssigkeit an der Wand die Strömungsgeschwindigkeit auf Null abnimmt, in- 
[folge des großen Geschwindigkeitsgradienten wesentliche Reibungswirkung. Außerhalb dieser 
Grenzschicht, wo keine so großen Gesehwindigkeitsgradienten auftreten, ist die Wirkung der 
/ähiekeit nur von sehr geringer Bedeutung, so daß sich dort annähernd eine Strömung von 
der Art derjenigen in einer reibungslosen Flüssigkeit einstellen wird, die nur durch Druck- 
wirkungen bedingt ist. Nimmt man die Grenzschichtdicke als klein an, so kann man durch 
eine Größenordnungsbetrachtung die strengen Navier-Stokesschen Bewegungsgleichungen 
vereinfachen, indem man die Glieder geringerer Größenordnung vernachlässigt. Auf diese 
Weise gelangt man zu der Differentialgleichung der Grenzschicht. 

Gemäß dieser Grenzschichttheorie nehmen wir nun auch bei der Strömung im ge- 
krümmten Rohr eine Grenzschicht von der Dicke ö an, innerhalb der die Reibung zu berück- 
sichtigen ist. Genügend weit vom Rohreinlauf ist d unabhängig von der Rohrlänge und nur 
eine Funktion des Winkels y (Abb. 3). Außerhalb der Grenzschicht haben wir einen 
Strömungskern, in dem wir Reibungslosigkeit annehmen. 

Wir wollen jedoch nieht die Differentialgleichung der Grenzschicht als Grundlage ver- 
wenden, sondern eine Impulsbetrachtung für die Grenzschicht als Ausgangspunkt für die 
Reehnung nehmen, eine Methode, die v. Kärmän'?) angegeben hat und die für eine Reihe 
von Fällen vonK.Pohlhausen'*) unter vereinfachenden Annahmen für die Geschwindigkeits- 
verteilung in der Grenzschicht durchgerechnet ist. 

Wir führen folgende Bezeichnungen ein (Abb. 3): x und y bzw. r und y sind die recht- 
winkligen bzw. Polarkoordinaten des kreisförmigen Rohrquerschnittes mit dem Radıus a. 

10) C. M. White: a. a.0. 


ti) I. Prandtl: Verhandl. des Ill. Intern, Math. Kongr. 1904. (Heidelberg), Leipzig 1905, S. 484. 
12, Th. v. Karmän: Über laminare und turbulente Reibung, Z. ang. Math. u. Mech. Bd. 1 (1921), S. 233. 


13) K. Pohlhausen: Zur näherungsweisen Integration der Differentialgleichung der Jaminaren Grenzschicht, 


7. ang. Math. u. Mech. Bd. 1 (1921). S. 252. 
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Der Winkel 9 ist die dritte Koordinate. Sein Scheitel liegt im Krümmungsmittelpunkt O0 des 
Rohres, sein Bogen steht senkrecht zum Rohrquerschnitt, in Richtung der Rohrachse. Die 
Strömung erfolgt in Richtung des wachsenden 9. Die Geschwindigkeitskomponenten im 
Strömungskern sind U, V, W. W ist die Hauptgeschwindigkeit, U und V sind die sehr kleinen 
(vgl. S.263) Komponenten der Sekundärströmung in Richtung von «= und 9. Die Sekundär- 
strömung in der Grenzschicht hat die Komponenten « in Richtung von r und » senkrecht 
zu r, in Richtung des wachsenden y, die Hauptgeschwindigkeit in der Grenzschicht ist w. 
Für die Geschwindigkeit W des Kernes, die wir, ebenso wie U, nur als Funktion von .r, 
also in y-Richtung konstant, annehmen, setzen wir unter Anlehnung an Versuchsergebnisse 
eine lineare Verteilung, vom inneren zum äußeren Rand des Rohres ansteigend, an: 
W=A+Bre=A-+DBrcosy ee ae oe ee 
(A = Geschwindigkeit in der Rohrmitte z mittlere Geschwindigkeit; B=t; 
linearen Geschwindigkeitsprofils). 


J 


Neigung des 


Dieser Geschwindigkeit entspricht ein Druckanstieg durch Zentrifugalkraft (pro Längen- 
einheit): 


Ö ee x R (»>) 


(o = Dichte). Der Krümmungsradius ist dabei groß gegen den Querschnittsradius angenommen. 


Mit 2 Ry erhält man daraus durch Integration: p = n [was + (Üz und mit Gl. (4) 
l 


’ es; r A B? 
et Ca: RR A’rcosy-+ABr’cos’y+ , 
J > 


r” cos’y)+ C Ra (6). 


P a Sal 
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Abb. 3. Abb. 4 und 5. | 


Der Druckunterschied zwischen Innen- und Außenwand wirkt sich in der früher geschilderten 
Weise in einer Sekundärströmung aus, die in der Mitte nach außen, in der Nähe der Rohr- 
wand nach innen verläuft. Diese Sekundärströmung findet bei nicht allzu schwachen 
Krümmungen hauptsächlich in der Grenzschicht statt. Für die weitere Rechnung nehmen 
wir an, daß die Rückströmung von außen nach innen, von der Mitte des Rohres in sym- 
metrischer Weise nach beiden Seiten hin, nur in der Grenzschicht stattfindet. Von der Mitte 
der Innenwand aus, wo sich beide Teilströme wieder treffen, verläßt die Strömung die Grenz- 
schieht und wandert strahlenförmig durch den Kern wieder nach außen, so daß am Rande 
der Grenzschicht dauernd Flüssigkeit in die Grenzschicht einmündet. Wir denken uns nun 
ein Element der Grenzschicht herausgeschnitten, das ın Strömungsrichtung durch zweı Ebenen 
im Winkelabstand dy und senkrecht zur Strömungsrichtung durch zwei Ebenen im Winkel- 
abstand dr abgegrenzt ist. Abb. 4 und 5 zeigen das Element im Auf- und Grundriß. Auf 
dieses Element wenden wir den Impulssatz an. Nach diesem ist bei stationärer Strömung 








“ 
I. un T, 
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der Impulsfluß dureh die das Volumenelement abgrenzende Kontrolltläche gleich der Summe 
der auf die eingeschlossene Flüssiekeitsmasse wirkenden äußeren Kräfte. Wir wenden den 
Impulssatz zunächst auf die y-Riehtung an (Abb. 4). Als neue Koordinate werde der Ab- 
stand EZ von der Rohrwand eingeführt, so daß 


BE Er A 
d 
Der durch F, in der Zeiteinheit eintretende Impuls ist o Rdg \v® dS,der durch F, austretende 
0 
Impuls ıst 
) | 1) 


. wr 
„ui 53 
) [1 { w | a 


” » Q ” 


| ds= Z I da „2 ds-+o de | Fe d vd: 


1) 
h Ö ö .< 
Der Unterschied beträgt demnach: o Rdyq ” | [v: d&)dy. Der Überschuß der durch F, aus- 
( ı 2 
/ « 


IE 
f) 

tretenden über die bei F, eintretende Flüssigkeitsmasse ist oRdy | (fr dz)dy. Aus Kon- 

C AITAN 

0 
tinuitätseründen muß ebensoviel Flüssigkeit dureh die Grenzfläche zwischen Grenzschicht und 
Kern in die Grenzschieht gelangen. Der damit eintretende y-Impuls ist, wenn U’ die y- 
Komponente der Kernströmung ist und der Index 9 den Rand der Grenzschicht, £=- d, kenn- 

) 


N\ 
. rr ( - .. * r . gf ® . . 
zeichnet: o Rdy I | I» dz)dıy. An äußeren Kräften kommen in Betracht: Die Reibung 
C 


N ( 
u 
> Q Ü Y > 
an der Rohrwand Tuerakdgdy ul, | aRdgydwy, der Unterschied der Druck- 
0 c/—=0 

.. gr * . V1 .. ’ Y > Ö yp . r r .n 

kräfte auf die Flächen F, und F, -oRdy y di» und die y-Komponente der Zentrifugal- 
( ’ 
) 

om“ ’ ’ 

kraft Rdg || np dysiny)d@. Die negativen Vorzeichen erklären sich dadurch, daß alle 
[) 


u 
kräfte entgegengesetzt zur y-Richtung wirken. 
Der Impulssatz für die y-Riehtung lautet also nach Kürzung mit dy und Rdy und 
Division durch o: 


() ) () 


(was u-[var or —[P) _ 2 sinyfmrae 8 
ii DaGe= vl > | sin [wi Ar ; r . (9). 
I, y VE, o \Oy/)s R Fr 

IE " 0" 


Innerhalb der Grenzschieht herrseht näherungsweise der Druck der angrenzenden Kern- 
strömung. Man kann deshalb für p den Wert an der Stelle r =a —- önach Gl. (6) einsetzen. 

Nun wenden wir uns zum Impulssatz für die g-Riehtung: In Abb. 5 ist der Grundriß 
des betrachteten Elementes aufgezeichnet. Da sich »w mit g nicht ändert, liefert die durch 
die Flächen F, und F, fließende Flüssigkeitsmenge keinen Impulsbeitrag, wohl aber die 


an der Übergangestläche von Kern und Grenzschieht eintretende Flüssiekeit ım Betrag 
() 


oW, Rdy | (fr d£)dy und die durch die Flächen F, und F, hindurchtretende Menge 
” Oy\, 2 
0 

2) 


N . 
der Sekundärströmung im Betrag o Rdgq | | "wde)dy. An äußeren Kräften sind hier 
€ a‘. 


N . 

on 

anzusetzen: Die Reibung an der Rohrwand — Tinsady Rdyg= «| 3 1 adyRdy und 
0&J% 


u 


der Unterschied der Druckkräfte auf die Flächen F, und F, im Betrag 


3 


| dg9dadır. 
- 

04/8 

Setzt man wieder den Gesamtüberschuß des aus dem Element austretenden Impulses über 
den eintretenden gleich der Summe der äußeren Kräfte, so erhält man 


\y io na en 
h | ‘ ‘8 i ON ao (op 
rwdE Wy,,,fed: val-;) VE TRTINTE 
I (FL 0 r\o, N 


u 
ee 
=. 


eb: 


.r 


Br a 
en 
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me \lacht man jetzt für die Geschwindigkeitsprofile » und mw geeignete Ansätze, die man in die 
len ;l. (8) und (9) einführt, so hat man zwei Differentialgleichungen zur Verfügung, aus denen 
Ah- nan die Grenzschichtdicke $ = (y') sowie eine weitere in den Ansätzen für die Geschwindig- 
<eitsprofile frei zu lassende Unbekannte bestimmen kann. Um die Rechnung nicht allzu 
7 | schwierig zu gestalten, geben wir den Profilen die in Abb. 6 bezeichnete Form von Parabeln 
' ınd lassen zu, daß bei » der Übergang von der Grenzschicht in 
u ‚lie Kernströmung mit einem Knick stattfindet. Die Ansätze, die 
len Randbedingungen genügen müssen, daß die Geschwindig- | % 
keiten an der Rohrwand den Wert Null haben und am inneren N 
Rand der Grenzschieht in die Randwerte der Kernströmung | 
übergehen, lauten: 
ei, rer "6 | 
N) N) (10) N} 
IS- “ > [A = T 
z & - [AUS2.E L ß „ 
»= W,|1-(1- ‚| (11). 
( Abb. 6%. 
11 - i u & z 2) 
ist die Geschwindigkeit der Sekundärströmung » im Abstand —, von der Rohrwand. Wir 
ıd haben früher die Voraussetzung gemacht, daß innerhalb der Sekundärströmung die Rück- 
y- strömung vom äußeren zum inneren Rand des Rohres nur in der Grenzschicht stattfindet. 
N» l,egt man durch den Rohrquersehnitt, z. B. durch den Mittelpunkt, einen vertikalen Schnitt. 
so besagt die Kontinuitätsbedingung, daß die durch den Kern nach außen wandernde Menge 
1% venau so groß sein muß, wie die oben und unten in der Grenzschicht von außen nach innen 
ließende Menge, daß sie aber wegen des bedeutend größeren ‚zur Verfügung stehenden Quer- 
schnittes eine entsprechend kleinere Geschwindigkeit besitzt, deren Rückwirkung auf die 
< Grenzschicht wir vernachlässigen wollen. Das in Gl. (8) vorkommende Glied mit U’ fällt 
deshalb fort. 
Mit den Ansätzen (10) und (11) ersetzen wir die in (S) und (9) vorkommenden Glieder. 
Diese werden im einzelnen: 
© d d 
[ ’ S sp: Er tD [ 21 S sur: 
Nds =: | £ : w’ de =. ”s 
’ 15 RES N) h 15 0 
1 I) [5 
) d 
» ai . ») N . .) # 
/ . | 2 On 2 W, 
vwd&öE=--06DW;: [va: — OD: | 3 5) 
| 15 Ö > VE), Ö ee 
( 0 
Außerdem ist nach Gl. (6): 
/A N 
( P\ ) Zul 52 a ep j 
| — (A )sın vw W?: | R 
Or/s Rr di oyr/a 
| Durch Zusammenfassung der einzelnen Glieder erhält man unter Vernachlässigung von & 


veren a, so daß a -d zZ a, zwei simultane Differentialgleichungen für dd (y) und D== Diy). 


0 0) I 15 D ” a a. > 5 .) 
N) Fer (ö D’)+ > val 8 Od B sın „y (A + Ba cos yr)? —=( a EEE 
Ö ‘ 7 v = 12] Ca 
DEN ($ D)(A—+ bacos y)+ 3 0% D Ba sın y — 9 0° Ova(iA—+bacosy)- 0 (14). 
( y’ . N 


/ur Lösung machen wir für die Unbekannten d und D Reihenansätze. Wegen der vorkommen- 
den Kreisfunktionen erscheint es zunächst zweckmäßig, Fourierreihen zu verwenden, aber 
wegen der zu bildenden Reihenprodukte, deren Koeffizienten hierbei selbst wieder unendliche 
Reihen werden, ist es vorteilhafter, Potenzreihen in y zu wählen. Die Betrachtung der 
Symmetrieeigenschaften der Strömung ergibt, daß die Grenzschichtdicke d symmetrisch zur 
-Achse ist, daß dagegen die Grenzschichtgeschwindigkeit D nach der negativen y-Richtung 
hin ihr Vorzeichen ändert, also unsymmetrisch zur «-Achse ist. Berücksichtigt man diese 
Eigenschaften, so hat man in der Potenzreihe für $ nur geradzahlige Potenzen von y, in 
der Reihe für D nur ungeradzahlige anzusetzen: 


) do +69 yr + od, yr+ 
D=-D,y+ D, yw+D,y’+ 


(19), 


— | — 
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Die beiden Integrationskonstanten für die zwei Differentialgleichungen sind dann schon durch 
die in den Ansätzen erfüllten Bedingungen: y=0:9=d,, D=V0 erledigt. Auch die Kreis- 


funktionen ersetzen wir dureh Potenzreihen: 


‚s ‚> ‚2 ‚Aa 
sin j' u ‘OS | AU A 
11 / f > > + Be GO: y - 5% (' + RN 


Mit diesen Reihen, die wir bei den geraden bis zur sechsten, bei den ungeraden bis zur 
siebenten Potenz führen, bilden wir die in den Differentialgleiehungen vorkommenden Aus- 
drücke und erhalten, nach Potenzen von y' geordnet, aus Gl. (13): 


1» ! | 
vi20, D’+>vra), SR 0, (A+ ba) 
i 1» 1 Wi 1 
B u" fa) OÖ, , D, 6b O, Q, DD: z ) rd N, | au I IE | a ° | d A B 1 S B? a” 
fi os A 'LBatl 
0,0, A+ ba'Y\| 
| ne p' ir I | 


15 
w”1129°D,D,+609°D’+23200,0,D, D,+80,0,D°’+4062°D?+ 5 va D. 





da A? 61 ıcda 
\? mn! ;- 2 | ag 2? 2 z 2 ee“ > 2 
ISO R Ö, | > Ihsaıll + > Pa STR OÖ, (A Ba) 
va Br RE va | „| 
ie I? Rp 0,0, | 2 a zu ta AB-+ 2 B’ a?) iR 0,0,(A+B us 
(16) 
Fy"|169,°D, D,.+ 160,’ D, D, + 250,0, D, D, + 14 0,0, D,’+240,0,D,D, 
r s i 1 
H-1005,9,D’+120°D, D,+1W00,0,D’+zgraD, 
ı 8 „.i# i Di _. . 
sel t6taAB+, BR a) 
‘ Ad y vo . ) ‘ 22 „2 
I) Br 2 u. zu balB+ 2 Ba 
va N (A? DEU u  ı | HZ R 
[> p 99 | 5 +4daAB- > B? a?) P53 7208| 5 +4aAB+ > b?’ a 
a a | 
), 6,(4A-+ Ba) 0,0,(4+ Ba’)| =V 
FRAME - TZHAAATARN 
und aus Gl. (14): 
z (Aa | 
0, D,(A—+- ba) 9» VW ralti—+ ba) 
‘) 
2 “ BE. (a 2 
"1305 D,(1+-Ba)+1t0,0,D, (1+Ba— n: 0, D,bBa- W 0,0, +9raBa 
j . 1) ‘ . 
Fy1502° D,(1+Ba)-+80,0, D,(A+Ba)+66,0,D,(A+Ba) 
>> ,. u S 
-30, D, (+ ba) = 0, D,ba-+ 6 0, D, ba +3 0,09, D, Ba 
(a (a ) | 
10 VÖ A) Ö, ‚raba 
o o° 12 | (1m). 


178° D.(A+Ba)+12 0,6, D,(A+Ba)+100,6, D,(A+Ba)+8 6,0, D,(A+Ba) 


13 
58° D,(A+Ba)+ 80,9, D,(A+Ba)+-,,9,' D, Ba 
“ ö j . i. 


VE u Bi 
75 0, D, Ba 1S 0,0, D, Ba 6 0, D,Ba-+ > 0,0%, D, Ba 
5) cr Ca (a 
—0,0,D.Ba+ -—0,.D Ba 10 0,0O 10 N), OÖ, 
1 0 4 1 ( 2 1 Oo 0 6 o . 

| 





+ yaDa 0 
f 


») 





Igew, 
lech. 
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irch 
Peis- 


zur 
\us- 
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Diese beiden Gleichungen sind nur dann identisch Null, wenn jeder Klammerausdruck Null 
wird. Diese Bedingung liefert für die 8 unbekannten Koeffizienten der Potenzreihenansätze 
Ss Bestimmungsgleichungen, von denen jeweils zwei unter sich lösbar sind. Der Klammer- 
ausdruck von y" und derjenige von y', beide gleich Null gesetzt, bilden zwei Gleichungen 
für 0, und D,. Mit den daraus bekannten Lösungen bekommt man weitere zwei Gleichungen 
für d, und D, durch Nullsetzen der Klammerausdrücke von y” und y" usw. Die jedesmal 
mit den vorauszehenden Lösungen versehenen Gleichungen sind, mit Ausnahme derjenigen 
für Ö, und D,, die quadratisch sind, linear. 


Zur Vereinfachung führen wir noch folgende Beziehungen ein: 


Da (a R 
(A-+ Ba)’ = — a? De ne 5 
A+ba p r ) 0 
Hierin bedeuten « und 5 dimensionslose Zahlen. @ hat die Dimension einer Geschwindig- 
keit. Das Minuszeichen erklärt sich daraus, daß € negativ ist, wie man sich durch Gl. (6) 
leicht überzeugen kann, da der Druck mit wachsendem 9 abnımmt, und daß (A—+ Ba)’ 
eine positive Größe ist. a und 5 treten an Stelle der bisher gebrauchten Größen A und B. 
Der Zusammenhang ist gegeben durch: 
| | 6 
A A G(1—P); >—=ap : i ; : i ; ; i . ; ; i i ‚ (19). 
a 
Mit diesen neuen Beziehungen lauten die dimensionslos gemachten Lösungen für die Potenz- 
reihenkoeffizienten folgendermaßen: 


da, ’ = va 15\? 
Im Abkürzung für 7 u: n = Abkürzung für | Daß? a’ + a5] 
Od, | 10 gi 
A Ga 75) 
n n-+ 55 


ö,  0,014583 m — B (0.129040 m + 0,092803 n + 0,249225) 
re 0.355303 m + 2.643939 n + 0.552656 


Ö | | 


Ö,\? 
4 52727 m + 3.737121 n + 5,743629 
,  DAOT2TB m + 2.004545 n — 0,616736 5) IST m + 3,7312 +5,043629] 


Ö, 1% | BI | a | 
+3 [0,029167 m + P (0,454646 mn + 4,337626 n + 1,054194)] 
Oo 


— [0,000729 m — PB (0,009675 m — 0,052936 n + 0,051717) — P° (0,146202 an 


+0,255,4n + 0.275607) 


NV. | 
VAGSWI m + 15TA6T DI n —-0,537925 
Ö . En | a 
. ni |0,02916%7 m + P (0,506263 m + 2,0044 u + 1,117567)] 
0 


| 


Ö, OÖ, 


—_ —_ — r} Ö, ’ - - f -_— - — .- 
+ [0,3572579 nm + 7,1019485 n + 360144] | s |0,255253 m —+- 4,705267 n + 6,240276] 
0, 0, on 


u, Mr a i ce 
7 (23 [0,014583 m 9 (O,227677 m + 4,531649 + 5,517646)] 


\ ‘ 


), FERNEN Eon 
: [0.001458 m + 9 (0.078573 m + 0,604031 nn — 0,115709) 
0 
?(0,183455 m + 2,015040 n + 0,953760)] 
+,[0,000017 m + 5 (0,000707 mı + 0,003641 nn — 0,005924) 
P*? (0,0838145 m + 0,022087 n + 0,051905) 


P°(0,049683 m ++ 0,126055 n + 0,0700] 
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Din 15 
TA ee >) 
N, | NV, y „>)e)r)e)»)» ) ) ) > za Seen > . ")iı1*) | 
D 15 Fu am: "dr dd.) I] u 124242) r P(O,uauuud - 0.303030) 
n er 
N. | /d,\° ü IR N), 
- 15°) (1533333 n + 3,257579) (1,200000 72 —- 1,151515) 
D, | | OÖ, OÖ, 
H ).) 
: RE PFUND \ - ; rt 
y DEURISHRERETEE 2532525) p (VOII N 0,01 1364) + P° 0,129630 n + 0.050505) 
( 
D IP) a.6761000 + 3904769 + [2 2 54 
n 15 \»,) (1,676190 1 + 3.904762) +5.) B (n + 1,540404) 
n .) 
N), Q, ER N). m > - 
2,545 + 4559827), (1,142857 0 + 0,649350) 
N), N), OÖ, 
N), ” a 
y pP (VOI2IG mn + 0,131455) + p VI2IR6 Rn + V02T7SUYS) ! 
( 


mm 


h 
y PO SHE nn + 0,567 100) + P(0,00277S nn V,000902) 


ı 
0 


- 5° (0,027469 0 + 0,011694) 9° (0,003086 n + 0,001203) 


Bei der Auflösung der quadratischen Gleichung für D, ist das positive Vorzeichen der Wurzel 
.. ® Er A ... 7 ® \ 
eültig, da nur dann D, für Rd, d. h. für den Fall des geraden Rohres, wo keine Se- 
3 

kundärströmung auftritt, gleich Null wird. 

Um sich über die Konvergenz, die zur praktischen Verwertung der Rechnung bis y=n 
erwünscht ist, Aufschluß zu verschaffen, ist es zweckmäßig, den Verkauf von d und D in 
ihren jeweils um ein Glied vergrößerten Teilsummen für einen bestimmten Fall graphisch 
aufzutragen. Da a und 9 vorläufig noch unbekannt sind, setzen wir dafür die dem Zahlen- 
beispiel entsprechenden Versuchswerte ein. Wir nehmen folgenden Fall aus den später zu 
erörternden Versuchen vorweg: 


0 Ad | 

- .: 5 

a „em, 7 100° 

r Q p ‘ 4 1 y . . - > .. 

( 0,464 - 10 * ke/em’ = 23,2 mm W.-S. auf 50 em Rohrlänge, 
RoOg 

| mittlere Durchflußgeschwindigkeit ec 17. 0tem/s. vr — 01265 em’s. 


Aus dem über dem horizontalen Durchmesser 
des Austrittsquerschnittes gemessenen Ge- 
schwindigkeitsprofil (Abb. 7) entnimmt man 
die mittlere Neigung B=tg y= 76.25 Is. 






\ 


a=11,54, p=0,H48. 


Mit diesen Werten errechnet man: 








va d, 0.1029 em D, 4,6261 em/s 
y. 
7 | | | | % 0.109733 D; 0.202571 
7 | | | | | | | j 0.015680 I, 0013150 
68 -— 1. 
Q, - en D, 
SER. 1.  ; u 
| q Fr Een $, D. 
Abb. 7. Ermittlung der mittleren Neigung VOIIOZHA j V.VOOOSOO, 


des Geschwindigkeitsprofils. Ö, D, 
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n Abb. 8 und 9 sind die einzelnen Näherungskurven, dureh Division mit ö, bzw. D, dimen- 
ionslos gemacht, aufgetragen. Man sieht, daß die bisher berechneten Glieder der Potenz- 
eihen nieht ganz ausreichen, um mit ihnen ein genaues Resultat bis y = zu erzielen. 
Wegen der alternierenden Vorzeichen der Reihenglieder weiß man aber, daß die folgende 
\äherung zwischen die beiden vorausgehenden zu liegen kommt, so daß man einen der Wirk- 
ichkeit angenäherten Verlauf einzeichnen kann. Bei dem Verlauf von D muß man außerdem 
berücksichtigen, daß der analytische Verlauf der Kurve bei y == nicht, wie bei y==0, durch 
Null geht, wie man auf Grund der Überlegung, daß beim Zusammentreffen der beiden Teil- 
<tröme eine Umlenkung in horizontaler Riehtung stattfindet, annehmen könnte. Man hat sieh 
vielmehr die analytische Fortsetzung der Lösung unserer Differentialgleichungen (13) und (14) 
so zu denken, daß sich beide Grenzschiehtströme bei y = überkreuzen und mit endlicher 
Geschwindigkeit in der gleichen Richtung, in der sie ankommen, weiterfließen, wobei sie 
unter dem ansteigenden Druck eine Verzögerung erfahren, bis schließlich eine Umkehr und 
eine Ablösung von der Wand eintreten würde. 














y 
Abb. 8. Näherungen für die Grenzscehichtdicke. 
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Abb. 9. Näherungen für die tangentiale Geschwindigkeit Abb. 10, 
der Sekundärströmung in der Grenzschicht. 


Der zur Berechnung der noch unbekannten Werte « und 9 einzuschlagende Weg ist 
'olgender: Zur Erfüllung der Kontinuitätsbedingung innerhalb der Sekundärströmung muß, 
wie man aus Abb. 10 leicht ersehen kann, die Gleichung gelten: 

) 


da $)sın ' vrdE= 


d 


ww IV 


DD. 


Daraus: 
4 VD 


r ON. 
3 la o)sın y 


Ifierin entsprieht einem Wert y ein Wert für U an der Stelle «(a  )siny'. Die im 
reibungslos vorausgesetzten Kern gültige Eulersche Gleichung liefert, wenn man berück- 


MW I W 
A . r ( : R s r DL 
sıichtigt, daß nach Voraussetzung y Od und daß R als klein niedrigerer Größenordnung 
( N ” 
vernachlässigbar ist, die Beziehung: 
 W () 
„ Ä 
I 


0.r 
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Daraus: 





oW Ö C 
£> £ (A+Be) B= em 
und mit Gl. (20): 
n 368 RUE | 
"ir "SD mittel 
u r 
3 Ca 1 Vogt 536 2 
J 0 Ö, D, D, Ö, 3 
Er D, 6, u mittel 


Kine weitere Beziehung gewinnt man, wenn man für ein Flüssigkeitselement, welches durch 
die Rohrwand und zwei Rohrquerschnitte im Abstand Rdoy begrenzt ist, das Gleichgewicht 
der Kräfte anschreibt. Die Reibungskräfte an der Rohrwand und die Druckkräfte auf die 
beiden Querschnittsflächen müssen im Gleichgewicht stehen. Diese Bedingung ergibt: 





T (d 7 
. . “| Ö p 
2Rdg «| 7, dı [| » Ip ı r de) rdrdy. 
0 "0 
Om W Op 
Daraus wird mit 7, Reibung an der Rohrwand : m - | —2u—% und 4 dp=CRdy 
| TIOdE), Ö üR7 
(nach Gl. (6)): 
"W, A "Bla 0) Can G@’n 
) dy= h) dy nö h) ans ydy tu I{v (22). 


Ersetzt man in Gl. (21) und (22) mit den Beziehungen (18) A und B durch a und /, so stehen 
zwei Gleichungen zur Verfügung, die prinzipiell zur Ermittlung von a und 5 genügen. Prak- 
tisch treten dabei aber sehr große Schwierigkeiten auf, weil beide Unbekannte in den Koeffi- 
zienten der für 6 und D einzusetzenden Potenzreihen in sehr verwickelter Weise auftreten. 
Will man das Hauptziel vorliegender Rechnung, nämlich die Bestimmung eines Widerstands- 
gesetzes, erreichen, so ist allerdings die Ermittlung von a und 5 sehr wichtig, weil nur 
daraus der Zusammenhang zwischen der Reynoldsschen Zahl und der Rohrkrümmung einer- 
seits und dem Druekabfall andererseits gewonnen werden kann. Um den Zusammenhang 
wenigstens angenähert zu bestimmen, wollen wir einige Vereinfachungen vornehmen, die es 
ermöglichen werden, eine asymptotische Gesetzmäßigkeit zu erhalten. 
Für kleine Winkel y kann man Gl. (21) ersetzen durch 


.) 2 


+) T 


' Br 
cr U 5 EEE ee 
Dieser Wert ist natürlich nicht für das ganze Intervall von O bis = zutreffend. Aus dem 
gemessenen Geschwindigkeitsprofil (Abb. 7) erkennt man aber, daß sich die Neigung des 
Protils nieht sehr stark ändert. Der Fehler, den man begeht, wenn man den Wert für B beı 
kleinen Winkeln als Mittelwert für das ganze Intervall nimmt, wird also nicht allzu groß 
werden. Auch Gl. (22) schreiben wir für kleine y' an, oder, was das gleiche bedeutet, wir 
nehmen $, als Mittelwert für ö zwischen O und =: 
7 7 
Af Ba f Ga 26, 
dıy 4 


cos vr dıy und daraus: A 
On. Ö, / . tv 37. tv 


0 u 


(24). 


Den dabei begangenen Fehler kann man aus Abb. 8, in der der richtige Verlauf von & auf- 
getragen ist, ersehen. Nun nehmen wir noch bei den Ausdrücken für ö, und D, eine Ver- 
einfachung vor, indem wir berücksichtigen, daß a* sehr groß ist (nach den Versuchen von 
der Größenordnung 10°) und setzen angenähert: 


O, | 0a | Id“ 4,21 | 

Ad (a va ,,/15 95 (ra iS, (ra /a 
| j | ESy Tan 637 > | y | 2R" y «| R 

I 


Pa Pa «A 


|| ca 15,’ 15 | va a 
| ern ; Bo AR = 
(Gi (| te 5‘ BR" = 0564 «| 7). 










































(21). 
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Jit diesen Ausdrücken und den Beziehungen (18) und (19) wird aus Gl. (23) und (24) nach 
infachen Zwischenreehnungen: 
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Widerstandsgesetz: Damit ist man in der Lage, das Widerstandsgesetz zu berechnen. 
Nach Dimensionsbetrachtungen ist die Widerstandszahl 
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‚Druckgefälle bezogen auf den Staudruck der mittleren Durchflußgeschwindigkeit). In diese 
(leiehung hat man als mittlere Durchflußgeschwindigkeit den für A errechneten Wert aus 
(4. (25) einzusetzen und erhält damit: 
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Bezieht man die Widerstandszahl 4 auf die mit dem gleichen Druckgefälle erzielte Wider- 





64 
standszahl ,== des geraden Rohres, so erhält man: 
a 
v 
. Ad Ia]l/a 
0,1064 | Re( 7) ee En ea na Ahr 2 
vo wiz 
\lan erkennt, daß das Ergebnis vorliegender Theorie ebenfalls den bei der Deanschen 
\1/ 
a\'? r .. i s \ \ a 
kechnung vorkommenden Parameter Re (75) enthält. Wieweit die auf Grund verein- 
; l 


[achender Annahmen gewonnene Beziehung (29) mit den Experimenten im Einklang steht, 
wird im Zusammenhang mit den Messungen im nächsten Kapitel gezeigt werden. 

Es wurde versucht, Punkte der Widerstandskurve in zweiter Näherung dadurch zu be- 
rechnen, daß man mit den Gl. (21) und (22) für spezielle Zahlenbeispiele die entsprechenden 
Versuchswerte A und B durch ein graphisches Iterationsverfahren verbesserte. Das Ver- 
fahren hat sich jedoch nicht als konvergent erwiesen. | 


Il. Experimenteller Teil. 


1. Die Versuchseinrichtung. Die Anordnung der Versuchseinrichtung ist aus Abb. 11 zu 
ersehen. Das Versuehsrohr a ist mit einem Flansch an den Wasserbehälter b angeschlossen. 
/ur Vermeidung von Störungen sind die Einlaufenden der Versuchsrohre abgerundet. Von 
den beiden anderen Anschlußstellen des Behälters ist die eine (e) mit einem Überlauf e ver- 
bunden, der für einen gleichbleibenden Wasserspiegel sorgt, die andere (d) mit einer Wasser- 
standsröhre k, an der die Höhe des Wasserspiegels beobachtet werden kann. Die Zuleitung 
(des Wassers erfolgt in f von der Wasserleitung aus. 

Die Rohre, an denen die Messungen vorgenommen wurden, sind gezogene Messingrohre 
mit einem inneren Durchmesser von 1 em. Das Biegen der Rohre wurde mit einer Walze 
vorgenommen. Um dabei eine Veränderung des kreisrunden Querschnittes zu verhindern, 
wurden die Rohre vor dem Biegen mit feinem Sand gefüllt. Es wurden 3 verschiedene 

a 1 1 1 


Krümmungen untersucht, und zwar —- r — 
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Zur Messung von Druckunterschieden wurde ein gewöhnliches Flüssigkeitsmanometer und 
ein Mikromanometer der in Abb. I1 sichtbaren Bauart'’) verwendet, welches gestattet, Drucke 
zwischen 0,02 und 50 mm W.-S. zu messen. 
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\bb. fl. Versuechseinriehtung. Abb. 12. Geschwindigkeitsmeßapparäat. 


Am Austrittsende des Versuchsrohres ist ein Gewindestück befestigt, mit dem das Rohr 
an den Geschwindigkeitsmeßapparat g angeschlossen ist. Die Konstruktion dieses Apparates 
ist aus Abb. 12 näher zu ersehen. Er trägt ein Pitotrohr aus Messing mit einer Öffnung von 
0,402 mm ©, dessen Ende unmittelbar vor dem Austrittsquerschnitt des Rohres in einer 
vertikalen Ebene vorbeigeführt werden kann. Die Entfernung betrug 0,2 mm. Um die Ge- 
schwindigkeit an jedem Punkt der Ausllußmündung messen zu können, ist das Pitotrohr in 
horizontaler Riehtung auf einem gut diehtenden Schlitten durch eine Schraube zu verschieben, 
in vertikaler Riehtung durch ein Ritzel mit Zahnstange. Die Einführung des Pitotrohres dureh 
den Schlitten in das Meßgehäuse ist durch eine Stopfbüchse gegen Wasserdurchtritt abge- 
‚liehtet. Dureh ein kleines Ventil kann man die im Gehäuse befindliche Luft entweichen lassen. 

Die Durchflußmenge ist durch einen hinter dem Geschwindigkeitsmeßapparat angebrachten 
Hahn %» regulierbar und wird in einem Meßzvlinder © mıt geeichtem Wasserstandsglas gemessen. 


2. Die Messungen. a) Geschwindigkeitsverteilung. Die Geschwindigkeiten 
wurden aus Staudruekmessuneen mit Hilfe der Gleichung 


ec y2gh=44,3] h|em/s] 


ermittelt, wobei 9=9,S1 m/s die Erdbeschleunigung und » den am Manometer als Differenz 
zwischen Gesamtdruek und statischem Druck in em W.-S. abgelesenen Staudruck bedeuten. 
Der statische Druck wurde durch eine Anbohrung in der Ebene des Ausflußquerschnittes 
neben dem ausfließenden Strahl gemessen. Nach den Ermittlungen, die Nikuradse') bei 
seinen Versuchen an geraden Rohren angestellt hat, sind infolge der außerhalb des Strahles 
lieeenden Anbohrune keine Abweichungen zu bemerken. Um ein Bild der Geschwindiekeits- 
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\bb. 13. Versuchsplan. Abb. 14. Druckanbohrungen. 


14, Genaue Beschreibung in J. Nikuradse: Gesetzmäßigkeiten der turbulenten Strömung in glatten Rohren. 
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verteilung über dem Querschnitt zu erhalten, muß man an möglichst vielen Punkten die Ge: 
schwindiekeit messen. Bei gekrümmten Rohren kann man sich wegen der Symmetrie der 
Strömung ober- und unterhalb des horizontalen Durehmessers des Querschnittes auf Messungen in 
einer der beiden Hälften beschränken. Die Geschwindigkeitsmessungen wurden nach dem in 
Abb. 13 dargestellten Versuchsplan vorgenommen. Wegen des steilen Geschwindigkeits- 
anstieges an den Rohrwänden, besonders an der Außenwand, sind die Meßpunkte dort dichter 
velegt. Die der Rohrwand am nächsten liegenden Meßpunkte sind dadurch bestimmt, daß 
‚lie Pitotrohröffnung gerade noch vollkommen im Flüssigkeitsstrahl liegt. Die Entfernung des 
Meßquerschnittes vom Rohreinlauf war an allen Rohren genügend groß, so daß die Messungen 
außerhalb der Anlaufstrecke in der voll ausgebildeten Strömung vorgenommen wurden. 

b} Druckabfall. Der Druckabfall wurde durch Anbohrungen der Rohrwand gemessen. 
In einem Meßquerschnitt sind 4 Anbohrungen angebracht (Abb. 14), von denen jede mit einer 
Schlauechtülle versehen ist. Der Druckunterschied zwischen Innen- und Außenseite des Rohres 
kann gemessen werden, indem man die beiden horizontal liegenden Tüllen eines Meßquer- 
sehnittes mit einem Differentialmanometer verbindet. Bei der Messung des Druckgefälles 
längs einer Rohrstrecke wurden zum Druckausgleich die 4 Anbohrungen einer Meßstelle durch 
Schlauchstücke und Glasrohrgabeln in Verbindung gebracht, so daß im Manometer die Differenz 
der mittleren Querschnittsdrucke gemessen wurde. 

e) Menge und Temperatur. Die sekundliche Durchflußmenge wurde durch Meß- 
zylinder und Stoppuhr ermittelt. Da sich die Zähigkeit des Wassers sehr stark mit der 
Temperatur ändert, mußten genaue Temperaturmessungen vorgenommen werden. Auf dem be- 
nutzten Thermometer war die 
Temperatur in '/Jio Grad ab- “4 IN 


lesbar. Bei der Messung der 
i . r r . x 
Geschwindigkeitsverteilung ’ | N 

R ö - e 4 re © 
mußte jede mit einer be- a 
stimmten sekundlichenDurch- . Me 
lußmenge laufende Versuchs- 08: > N 
reihe bei konstanter Rey- N N 
noldsscher Zahl durchge- IN \ as 
Do. 2 an, N my ’ ua er R I } I I \ ar > N 
führt we rd n. ls war de shalb S | | 7 NE 
erforderlich, daß sich die | | N ’ 
Temperatur des Wassers | FE | N 
während dieser Zeit nicht | | N 
wesentlich änderte. Es wur- 
den höchstens Anderungen 
von 0,2°C zugelassen. 
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3. Auswertung der Ver- » 7. 
suche. a) Das Wider- i- Yr 
standsgesetz. In den a HE 1 Su 
Abh. 15. 16 und 17 igt die Abb. 15. log (110) in Abhängigkeit von log Re für R sw 
Abhängigkeit der  Wider- z 
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stab aufgetragen. Die Mes- | | Mn DS 
sungen sind bis in das tur- | IN | Don 
bulente Gebiet ausgedehnt. 
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für das gerade Rohr ıim | 
laminaren Gebiet gültige um | 
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(Gerade mit der Neigung 1:4) eingetragen. Die Widerstandskurve bei laminarer Strömung 
liegt überall höher als beim geraden Rohr, und zwar um so höher, je stärker das Rohr ge. 
krümmt ist. Mit sinkender Reynoldsscher Zahl nähert sich die Widerstandskurve derjenigen 
des geraden Rohres, was dadurch bedingt ist, daß die Sekundärbewegung bei kleinen Ge- 
schwindigekeiten weniger widerstandserhöhend wirkt, weil die Wasserteilchen einen größeren 
Winkelabstand in Riehtung der Rohrachse benötigen, um eine geschlossene Sekundärstrom- 

linie zu durchlaufen, oder 


\ . | - weil innerhalb eines ge- 
| N ; | ebenen Winkelabstandes 
| N a8 | | | | | | die Anzahl der „Schrauben- 

\ u | eänge* eeringer ist. 
| | x Na | Eine weitere Erschei- 
RN | al | VERNDEENDE GENRE RER, nung ist die Erhöhung der 
us x | | kritischen Reynoldsschen 
4 a Zahl. Der Übergang in die 
3 un Pr Sl Turbulenz findet um so später 
\ | a RR statt, je stärker die Krüm- 
N | Fe mung ist. Diese stabili- 
N I 01 | '  sierende Wirkung der Krüm- 


mung wurde auch in der 

u | | | | | | | | früher erwähnten Arbeit von 
| White festgestellt. Bei gut 

| abgerundetem Einlauf und 





3 75 Y guterberuhigungdesWassers, 
2 die dadurch gewährleistet 
» Id \ r .. 
Abb. 17. log (1004) in Abhängigkeit von log Re für 2, K war, daß das W asser übeı 
ı .«) 


ein Sieb in den Behälter 
gelangte, wurden folgende 


kritische Reynoldssche 
Zahlen beobachtet: 

A 1 | | 

R | 50 100 200 





Re, | 3620 4730 3980 
Die beim geraden Rohr be- 
sonders auffallende plötzliche 
Widerstandserhöhung beim 
Eintritt der Turbulenz wird 
mit zunehmender Krümmung 
immer mehr ausgeglichen, so 
daß allmählich ein stetiger 
Übergang in die turbulente 
Widerstandskurve stattfindet. 
Das kommt daher, daß der 
Widerstandsunterschied ge- 
genüber dem geraden Rohr im 
laminaren Bereich wesentlich 
erößer ist als im turbulenten und daß das bei zunehmender Krümmung erfolgende Anwachsen 
des Widerstandes bei laminarer Strömung rascher erfolgt als bei Turbulenz, bei der die 
Widerstandskurve nieht sehr beträchtlich von derjenigen des geraden Rohres abweicht. Diese 
Erscheinung mag gleichbedeutend sein mit einem von Millikan'") erzielten Ergebnis. Der in 
seiner Arbeit dureh näherungsweise Integration der Kärmänschen Differentialgleichung der 
Grenzschieht errechnete Reibungswiderstand für eine gekrümmte Oberfläche ist im Jaminaren 
Gebiet um 28°/o, im turbulenten Gebiet nur um 7", höher als die entsprechenden Reibungs- 
widerstände der Platte. 











Abh. IN. 


In Abb. 18 sind die Widerstandsmessungen, und zwar die Werte „ ın Abhängigkeit 
’0 
1 


a\ls. ' 
von dem Parameter Re| 7) im logarithmischen Maßstab aufgetragen. Die Meßpunkte ım 
l 


laminaren Gebiet ordnen sieh wie bei White für alle Krümmungen längs einer einzigen Kurve 


16) €‘, B. Millikan: The Boundary Layer and Skin Frietion for a Figure of Revolution. Trans. Amer. Soc. 
Mech. Eng. Applied Mech... Bd. 54 (1932), Nr. ?, S. 29. 
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Sobald die Strömung turbulent wird, weichen sie von dieser Kurve ab, und zwar um 























Ing an. 
ge. so später, je stärker die Krümmung Ist. Die Abweichungsstellen liegen jeweils bei der dureh 
en lie Krümmung bestimmten kritischen Reynoldsschen Zahl. In der Abbildung ist die von 
(Ge White experimentell ermittelte Widerstandskurve eingezeichnet. Die durch unsere Versuchs- 
ren punkte gelegte Kurve fällt mit ihr zusammen. Die Abszissenachse stellt das Verhältnis 
a ’ 1/, 
i /. : ; : Ad 2 
an i dar, das dem geraden Rohr entspricht. Bei ganz kleinen Parametern Re (7) unter- 
der ‚ i 
ge- scheidet sich die Widerstandskurve überhaupt nicht vom Widerstandsgesetz des geraden Rohres. 
/ 1 
les . a 2 er Aa | r 7 E ! 2 
Gr rst etwa be Re -) —= 12,6 (für : Re =- 250) beginnt sie davon abzuweichen. 
m. Erst etwa beı [7 R 100 | i 
Die Kurve. die das durch näherungsweise Integration der Navier-Stokesschen 
el- Bewerungsgleichungen gewonnene Widerstandsgesetz (3) 
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on | I £ " s 
‚ut 100; Re 350) überein. Von da ab weicht sie sehr rasch ab. 
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se Abb. 19 und 19a. Geschwindigkeitsverteilung und Isotachen. 2 n. Re = 1930, laminar. 
. .«) 
im 
C N i . ’ . 
h Die mit der Grenzschichtbetrachtung gewonnene Beziehung (29) 
en 
9 \1 1 
n / ' a \"/2] 2 
lie — — 0,1064 Re| ;) 
Se PR | R 
in wird in der logarithmischen Auftragung durch eine Gerade dargestellt, der sich die Wider- 
ler | a\ls , > 
en standskurve mit wachsendem Parameter Re( ,,) immer mehr anschließt. 
) 
Be Prandtl'”) hat für die Widerstandskurve die Näherungsformel 
11a] 0,36 
a\': 
0,37 [Re | 
eit FR wis 
w angegeben, die ebenfalls in Abb. 15 eingetragen ist. 


b) DieGeschwindigkeitsverteilung. In den Abb. 19 bis 23 sind Geschwindigkeits- 
ve profile, gemessen über horizontalen Linien, bei laminarer und turbulenter Strömung für 
3 Rohre verschiedener Krümmung aufgezeichnet. Die Geschwindigkeiten e sind mit der 


17), L,. Prandtl: Abriß der Strömungslehre, Braunschweig 1931, S. 121. 
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mittleren Durehflußgeschwindigkeit e, der Abstand von der Rohrwand mit dem Rohrdurch- 
messer d dimensionslos gemacht. Das Maximum der Geschwindigkeit ist überall in die 
Nähe der äußeren Rohrwand geschoben. Die turbulenten Geschwindigkeitsprofile unter- 
scheiden sich von den laminaren dadurch, daß im Kernstrom der Anstieg von innen nach außen 
viel flacher verläuft und daß der Geschwindigkeitsabfall in der Nähe der Rohrwand sehr 
stark ist. Man erkennt auch hier die ausgleichende Wirkung der Turbulenz auf die 
Geschwindigkeitsverteilung, wie es beim geraden Rohr bekannt ist. 


Kin übersichtlicheres Bild der Geschwindigkeitsverteilung über dem Rohrquerschnitt 
erhält man durch Eintragung der Linien gleicher Geschwindigkeit, die als Isotachen bezeichnet 
werden. (Abb. 19a bis 23a). Die Geschwindigkeiten sind auch hier dimensionslos gemacht. 

. . . . r 6 . » . 
Die Isotachen sind in gleichen Abstufungen 0,2 eingetragen, so daß aus ihrem Abstand 
Ä = 


die größere oder geringere Steilheit des Geschwindigkeitsanstieges zu ersehen ist. 


‚ » or o . \ 
$ y Pe. 
EINE EEE EN ae oO \ \ ? 
K- n 
" fr \ a 





. R s Aa l 2 
Abb. ?0 nnd ?da. Geschwindigekeitsverteilung und Isotachen. .Be 540, turbulent. 
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Abb. ?I und ?la. Geschwindigkeitsverteilung und Isotachen., 
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Zusammenfassung. Ein durch näherungsweise Integration der Navier-Stokesschen 
Bewegungsgleichungen für die laminare Strömung im gekrümmten Rohr gewonnenes Widerstands- 


«nr 
a\'? _ En “ 
resetz, das vom Parameter Re | [Re Reynoldssche Zahl, K rümmungsverhältnis) 


a 

R R 

abhängig ist, ist nur für sehr kleine Reynoldssche Zahlen und Krümmungen brauchbar. 
Um auch für große Reynoldssche Zahlen im unterkritischen Gebiet eine Gesetzmäßigkeit für 
den Widerstand zu erhalten, wurde eine Rechnung durchgeführt, die sich auf die Prandtlsche 
Grenzschichttheorie gründet. Die Lösungen der durch eine Impulsbetrachtung an der Grenz- 
schicht gewonnenen Differentialgleichungen ergeben für den Widerstand ein asymptotisches 
Gesetz, dem sich die experimentell ermittelten Werte in der Grenze gut anpassen. Auch bei 


" EEE a\ 
diesem Gesetz ergibt sich eine Abhängigkeit von dem Parameter Re | -,;) 
$ l 

Die Widerstandsmessungen erstrecken sich bis in das turbulente Gebiet. Außerdem 

sind die Geschwindigkeitsverteilungen über dem Rohrquersehnitt für verschiedene Krümmungs- 
verhältnisse bei laminarer und turbulenter Strömung experimentell ermittelt. 152 
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Abb. 22 und 2?a. Geschwindigkeitsverteilung und Isotachen. R 
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Abb. 23 und 23a. Geschwindigkeitsverteilung und Isotachen. = Re 220, Jaminar. 
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Plarimetrische Auswertung von Stieltjesintegralen. 
Von E. J. Nyström in Helsingfors. 
$ 1. Auswertung durch Planimetrieren unter Benutzung eines zweckmäßigen Koordinatensystems. 


Die nieht selten auftretenden Integrale der Gestalt 


b h 
‚Aaalh oder kur IFA . . > 2 e-ner 


wo fit) und A(t) irgendwie gegebene Funktionen von f, ferner a und b Konstanten sind, 
heißen Stieltjesintegrale'). Für ihre angenäherte Auswertung ist es im allgemeinen 
nieht vorteilhaft, den Übergang zu Integralen der Form 
h 
a a 


a 


dureh direkte Approximation des Produktes 
BET 8 ern ar er er 


zu vollziehen, weil in den Anwendungen die auszuwertenden Integrale meistens primär als 
Stieltjesintegrale auftreten, sondern man approximiert besser den einen Faktor der rechten 
Seite von (9) ”). . 

Bei planimetrischer Auswertung von (1) kann aber der Übergang zu (2) unmittelbar durch 
die Darstellung der Funktion fif) in einem rechtwinkligen Koordinatensystem besonderer Art 
erfolgen. 

Wir gehen von der Darstellung der Funktion A (f) durch ihre Funktionsleiter?) aus. 
Diese entsteht, wenn man von einem Punkt O einer F-Achse Strecken % (t,), h (t,), ... abträgt 
und die Endpunkte mit den Argumentwerten #,, f,, ... beziffert. In Abb. 1 ist die Funktions- 
leiter für A (f)=yt angegeben. Mit der benutzten Längeneinheit wird auf der y-Achse eine 
eleichförmige Teilung konstruiert. In dem so definierten Koordinatensystem kann irgend- 
eine gegebene Funktion f(f) graphisch dargestellt werden. 

Wir behaupten, daß die in dem genannten Koordinatensystem gezeichnete Kurve y=f(h) 
zusammen mit ihren Ordinaten in den Endpunkten a, b des Intervalls und dem dazwischen 
liegenden Stück der f-Achse eine Figur bildet, deren Flächeninhalt eben den Wert (1) hat. 
Dies zu beweisen ist sehr einfach. Wenn x die Entfernung des Punktes ? der t-Achse vom 
Ursprung 0 ist (Abb. 1), so sind die Koordinaten irgendeines Punktes der gezeichneten Kurve 


BEN. REF: - 5. © 0 wu en An te 


und die bekannte Formel für den Flächeninhalt ergibt nun direkt den Ausdruck (1). Das 
h 


Integral (1) zeigt somit, als Flächeninhalt gedeutet, gegenüber dem Integral \f(f)dt keinen 


a 
anderen Unterschied als den, daß auf der f-Achse die Funktionsleiter Ah (f) angebracht ist statt 
einer eleiehförmigen Teilung. Wir können also sagen: Um das Stieltjesintegral (1) 
planimetrischzu finden, versieht man die f-Achse mit der Funktionsleiter Äk(f) 
und verfährtdann,alsob es sich umeingewöhnliches Integral der Funktion 
fi) handelte. 

Die Funktionswerte f(#) werden als Ordinaten gewissen Punkten der t-Achse zugeordnet, 
da diese Zuordnung durch die Funktion k(f) bestimmt wird, so kann man im Vergleich mit 
der gewöhnlichen Interpretation von einer Abszissentransformation sprechen. 

Sind f(f) und h(t) tabellarisch gegeben, so erfordert die betreffende graphische Dar- 
stellung nicht wesentlich mehr Zeit, als die gewöhnliche in rechtwinkligen Parallelkoordinaten. 
Wenn aber eine oder beide Funktionen graphisch gegeben sind, entsteht die Frage, ob sich 
nieht ein Umzeichnen der Figur lohnt, da die Methode sich sonst durch große Einfachheit 
auszeichnet). In den Anwendungen dürfte insbesondere der Fall oft vorkommen, daß Ahif) 
bei einer Reihe verschiedener Funktionen f(f) stets dieselbe ist, und dann vereinfacht sich 
natürlich das Verfahren, da man dasselbe Koordinatensystem mehrfach anwenden nnd vielleicht 
im voraus herstellen kann. 

ı\ Ein bestimmtes Integral der Form (1) kann als Grenzwert einer Summe definiert werden, falls /(t) in (a,b) 


stetig und A dt) von besehränkter Schwankung ist (vgl. z. B. E. Pascal: Repertorium der höheren Analysis, 2. Aufl., 
Bd. I, 3, Leipzig und Berlin 1929, S. 1328). 

®\ Hellerieh und Sehmidt: Numerische Auswertung von Stieltjesintegralen (Journal für die reine und 
angewandte Mathematik, Bd. 164, S. 243 bis 255 (1951)). Orloff: Über die angenäherte Auswertung von Stieltjes 
integralen (Commun. Soe. Math. Kharkow et Inst. Sci. math. Ukraine, IV, S. 6, 47 bis 59 (1935). 

3) V\el. z. B.Schwerdt: Lehrbuch der Nomographie, Berlin 1924, S. 15. 
%) Betr. das praktische Umzeichnen vgl. man die Fußnote ®). 
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Als Beispiel betrachten wir das Fresnelsche Integral’) 

z t t 
. 
sın 

| sın 2 | 2 2. / 5 
S (2) d2=7> —dti= \sın tayt eur © 

I > 2 \t - 

2a) ya 


Wenn wir bis z=r7,d.h. t=2, integrieren, kommt es auf die Ermittlung des ersten 
\aximalwerts der Funktion S (z) an. 


Der Bogen 0<t=2 der Sinuslinie y= sin, f wurde vom Verfasser mit der Leiter der 


Funktion k(h)=yYt dargestellt, wodureh sich die in Abb. 1 sichtbare Kurve ergab; die Längen- 
einheit war 10 em. Bei der Planimetrierung der von der Kurve und der t-Achse begrenzten 
Fläche wurde der Wert S(z) mit einer Genauigkeit von etwa 0,3°/, erhalten. Der auf drei 
Stellen abgerundete Wert ist S(z)=0,714. Die Abweichung ist natürlich sowohl auf Un- 
senauigkeiten der Kurvendarstellung wie auf solche der Planimetrierung zurückzuführen. 
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$ 2. Ein Spezialinstrument zur Auswertung von Stieltjesintegralen. 


Da das Umzeichnen der Kurve y==f(f) mit einer gewissen Mühe verbunden ist, erscheint 
es angemessen, nach Methoden zu fragen, die das Stieltjesintegral (1) direkt durch Planı- 
metrieren zu bestimmen gestatten, falls die Kurvenbilder f und h der in (1) auftretenden 
Funktionen schon gezeichnet, bzw. photographiert vorliegen. Um ein solches Verfahren zu 
finden, liegt es nahe, die in $ 1 besprochene Abszissentransformation durch eine Bewegung 
des Zeichenblattes zu ersetzen, und eben diese Idee liegt den folgenden Betrachtungen zugrunde. 

Die erwähnte Bewegung sowie die sonstige Anordnung kann mit recht einfachen Mitteln 
realisiert werden, und zwar ist die ganze Apparatur in Abb. 2 ziemlich genau im Anschluß 
an die vom Verfasser vorgenommenen Versuche dargestellt. Vielleicht läßt sich dieselbe noch 
vereinfachen, bzw. zwecks Erlangung größerer Genauigkeit abändern, worauf wir noch zurück- 
kommen werden. 

Der Einfachheit wegen möge angenommen werden, die benutzte Längeneinheit sei für 
beide Kurven dieselbe. Später machen wir eine Bemerkung über den Fall, wo diese Be- 
dingung nicht erfüllt ist. 

Die beiden Zeichenblätter mit den Kurven f und %k seien so miteinander verbunden, bzw. 
auf einem dritten Blatt = befestigt, daß die f-Achsen parallel sind und die y-Achsen aufein- 
anderfallen, olıne daß die Kurven sich schneiden. In Abb. 2 sind die in $ 1 besprochenen 
T 
») 


abgebildet, wegen Platzersparnis sind die positiven Richtungen der y-Achsen einander ent- 
gegengesetzt gewählt. 

Das Zeichenblatt = wird zwischen ein (horizontal liegendes) Reißbrett und eine feste 
Glasscheibe y gebracht, und zwar so, daß die t-Achsen einer Kante k des Reißbretts parallel 
werden (Abb. 2). Damit das Papier 7 mit Hilfe einer Klemmvorriehtung « und einer Lauf- 


Funktionen f(f=sin , t und hf) yt auf diese Weise in zwei Koordinatensystemen 0 (t, y) 


5) Vgl. Jahnke-Emde: Funktionentafeln, Leipzig und Berlin 1909, S. 231. (Die 1935 gedruckte neue Auflage 
ist an der betreffenden Stelle weniger ausführlich. 
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schiene in der Riehtung von k verschoben werden kann, ist die Glasscheibe y von zwei randlich 
zelezgenen dieken Papierstreifen o gestützt. Auf der Glasplatte befindet sich ferner ein in 
dem Reißbrett befestigtes Zelluloidlineal /, dessen eine Kante ! um 45° gegen k und die 
Riehtunz der f-Achsen geneigt ist, sowie ein bewegliches rechtwinklig-gleichschenkliges 
Zelluloiddreieck m. Dieses kann mit seiner einen Kathete längs der Kante ! von 4 gleitend so 
bewegt werden, daß seine Hypotenuse gegen k stets rechtwinklig bleibt. 

Das nun beschriebene Gerät in Verbindung mit einem gewöhnlichen 
Planimeter gestattet die Auswertung von Stieltjesintegralen. 

Wir betrachten die geradlinige Bewegung des Zeichenblatts z, nachdem dasselbe in 
der Klemmvorriehtung u befestigt und der Parallelismus der t-Achsen und %k kontrolliert 
worden ist. Bezeichnen wir mit FH einen Schnittpunkt der Kurve k mit der Linealkante 1, 
so hat die Verschiebung des Papiers = zur Folge, daß FH die Kurve h stetig durchläuft. 
Wir denken uns nun das Dreieck ® so bewegt, daß der eine Endpunkt seiner Hypotenuse 
stets mit dem Punkt F/ zusammenfällt; die -parallel bleibende Hypotenuse bestimmt dann 
auf der Kurve f einen Punkt F, und eben diesem beweglichen Punkt soll der Fahrstift des 
Planimeters folgen. 

Um den nieht gezeichnet vorhandenen Ort g des Punktes F (in der Abb. 2 gestrichelt) 
zu untersuchen, bedienen wir uns eines ruhenden Koordinatensystems (x, 9), dessen x-Achse 
mit der zur Kurve f gehörigen t-Achse zusammenfällt und dessen Ursprung O0 als Orthogonal- 
projektion des Sehnittpunktes € der Geraden ! und der « Achse der Kurve Ah bestimmt ist. 
In diesem System ist die Abszisse von F gleich «= 0A=C0B=BH=h(f), wo t die gemein- 
sıme Abszisse der Punkte HM und F in den auf z gezeichneten beweglichen Koordinaten- 
systemen ist. Die Ordinate von F im System (“, y) ist offenbar „= AF=f(h). Der Punkt F 
beschreibt somit in bezug auf die ruhende Ebene, das Reißbrett, bzw. die Glasscheibe y die 
Kurve g mit den Gleichungen 





DRAN ET ur tete te 


die ja mit der Kurve (4) identisch ist. 

Aus dem eben Gesagten erhellt die Wirkungsweise des Instruments. Es wird in der 
Tatdieins$1betrachteteFigurplanimetriert, nur liegt die krummlinige Begrenzung 
derselben jetzt nicht gezeichnet vor, sondern die beweglichen Ebenen = und & bestimmen in 
jedem Augenblick die richtige Lage des Planimeterfahrstiftes®). 

Unser Apparat kann zwar nicht als kompliziert bezeichnet werden, wohl aber seine 
Handhabung. Es müssen in der Tat zwei Personen bei seinem Gebrauch tätig sein; die eine 
hält den Planimeterstift im Punkte F, während die andere mit der linken Hand für die 
Parallelbewegung des Papiers sorgt und mit der rechten eine Winkelspitze des Dreiecks 
die Kurve A entlang führt. Beide Führungen sind wesentlich schwieriger zu bewerkstelligen 
als ein einfaches Entlangführen des Planimeterstifts längs einer Kurve, und sie erfordern 
einige Übune. 

Nach den Versuchen des Verfassers ist bei Auswertung von Stieltjesintegralen (1) eine 
(‚senauigkeit von etwa 1°/, verhältnismäßig leicht zu erreichen. 

Eine gewisse Erleichterung in dem Gebrauch des Instruments wird erzielt, wenn man 
durch Schlitze oder andere Anordnungen verhindert, daß der Planimeterstift und die Winkel- 
spitze ]/ des Dreiecks ® sich von den ihnen vorgeschriebenen Geraden überhaupt entfernen 
können. Der naheliegende Gedanke, die Parallelbewegung des Papiers automatisch zu be- 
werkstelligen, scheint von zweifelhaftem Nutzen zu sein, wenn man beachtet, daß es ein Vorteil 
ist, die Geschwindigkeit des Papiers nach Belieben regeln zu können. Besser scheint die Idee, 
das Planimeter auf einem in der #-Richtung beweglichen Wagen aufzustellen, dieselbe kann 
aber nicht mit ganz einfachen Mitteln realisiert werden. 

Über die Führung ist des näheren folgendes zu bemerken. Die Hypotenuse des Dreiecks » 
verschiebt sich bei monoton verlaufender Funktion h stets in demselben Sinn, dagegen ändert 
sich im allgemeinen der Sinn der Verschiebung des Papiers z in jedem Punkt D, wo die 
l,inealkante ! die Kurve Ah berührt. Man ist geneigt, eine solche Stelle als eine besonders 
schwierige anzusehen, da der Punkt MH als Schnittpunkt von h und I überhaupt in einer mehr 
oder weniger weiten Umgebung von D schlecht bestimmt ist. Zu beachten ist aber, daß die 
lage D des Punktes MH einer Umkehrstelle in der Verschiebung des Zeichenblattes 7 ent- 
spricht, wo seine Geschwindigkeit gleich Null ist. Der Planimeterstift wird also längs einer 
Kurve auf nahezu ruhender Unterlage verschoben, was natürlich leicht ist: die Kurven g und h 
decken sich in der Umgebung der betreffenden Stelle. 

%, Um die Kurve g wirklich punktweise zu zeichnen, braucht man nur ein Stück durehsiehtiges Millimeterpapier 


auf die Glasscheibe y zu legen und von beliebig gewählten Punkten H der Kurve h nach Punkten F der Kurve / 
heraufzuloten. 













































ngew., 
Mech. 


dlich 
n in 

die 
liges 
dl so 


hen 


e in 
liert 
be |, 
iuft, 
1USe 
lann 
des 


ıelt) 
:hse 
nal- 
ist. 
ein- 
ten- 
t F 
die 


(6). 


ler 
Ing 
in 


ine 
ine 
die 

() 
ren 
rn 


ne 


an 
e]- 
en 
)e- 
>] 
’e, 
ın 


ie 


Band 14, Heft 5 


Oktober 1934 279 


Weinig. Beitrag zur Theorie dünner schwach gewölbter Tragflügelprofile 





Schwierigkeiten treten dagegen auf, wenn die Kurve h nahezu parallel der y-Achse ver- 
läuft, also wenn die Funktion A (ft) ein sehr heftiges Steigen bzw. Fallen aufweist. 

Wären die Kurven f und A mit verschiedener Längeneinheit der y-Achsen gezeichnet, 
so brauchte man nur die Neigung des Lineals / gegen k zu ändern. Man könnte auch durch 
mechanische Hilfsmittel auf eine eventuelle Verschiedenheit der Längeneinheit der t-Achsen 
sowie auf eine schräge Lage der Kurvenbilder Rücksicht nehmen, hierauf wollen wir aber 
nicht eingehen. 

Der Apparat kann insbesondere zur Bestimmung von statischen Momenten und Träg- 

b b 
heitsmomenten [tfil) dt bzw. | fit) dt in bezug auf eine y-parallele Gerade dienen. Hierbei 
da Ad 

p E.2 s 4 m 
ist ersichtlich A) = 5 F, bzw. h()- — F zu setzen. Auch könnte man damit Fourier- 

_ +.) 
koeffizienten bestimmen. Werden die entsprechenden Führungen mittels Schlitze bewerk- 
stelligt, so kann der Apparat von einer Person gebraucht werden. 

Schließlich ist zu bemerken, daß der Apparat nicht nur zur Auswertung bestimmter 


Integrale brauchbar ist, sondern, wie auch das einfache Planimeter, als Integrator’), 
! 


so daß die Planimeterrolle stetig den Wert der Funktion \f(f) dh (t) entweder direkt oder nach 


da 
Hinzufügen einer unmitttelbar abzulesenden Korrektion anzeigt. Diese Anwendungsmöglich- 
keit beruht darauf, daß wir über die Aufstellung des Planimeters keine Voraussetzungen nötig 
rehabt haben, außer der selbstverständlichen, daß es auf dem Reißbrett für seine Bewegung 
veenügend Platz finden muß. 148 


Beitrag zur Theorie dünner schwach gewölbter Tragflügel- 
profile. 
Von F. Weinig in Berlin-Charlottenburg. 
Aus dem Institut für technische Strömungsforschung, Technische Hochschule Berlın. 


I. Die Abbildungsfunktion. Durch die Beziehung 


wird ein Streckenprofil der z=x#+iy-Ebene in einen Kreis der {=S-+in-Ebene ab- 
gebildet. Das Streckenprofil fällt hierbei in die «-Achse; Anfangs- und Endpunkt haben die 
Abszissen +1 und 2=-—1. Der Bildkreis hat seinen Mittelpunkt in &=0, 7 =0 und 
hat den Radius r=1. Für Punkte des Bildkreises ist also Z(k) = ei«, wenn a den Winkel 
zwischen dem Fahrstrahl nach dem betreffenden Punkte und der &-Achse bedeutet. 


Zwischen den Punkten des Streckenprofils 2 .x und denen des Bildkreises (= ei« be- 
steht die Beziehung: 


e=—fea+e-ie)=cosa. 


Bildet man nun ein dünnes schwachgewölbtes Profil, dessen Endpunkte mit denen des 
oben beschriebenen Streckenprofils zusammenfallen, in gleicher Weise ab, so wird die ent- 
sprechende Kontur k’ um so weniger vom Bildkreis des Streckenprofils abweichen, je weniger 
vekrümmt und je dünner das Profil ist. Im Grenzfalle werden also diese Abbildungen zu- 
sammenfallen. 


7) Man sehe die Arbeit des Verfassers: Anwendung des Planimeters als Integrator (Soe. Seient. Fenniea, Comm. 
Phys.-Math. VII. 10), Helsingfors 1934, 25 8. 


Problem, das schon Munk (Nat. adv. comm. f. aeronauties 


*, Anmerkung. Die Arbeit betrifft das gleiche 
Der Vorteil der hier mitgeteilten Ergebnisse 


Rep. 142 und 191) und Birnbaum (ZAMM 1923, S. 190) behandelt haben. 
besteht darin, daß man bei der angegebenen Entwieklung nur die ersten Glieder (vo, "1, ?») zur Kenntnis der Haupt 
eigenschaften des Profils (Anstellriehtung für glattes Umströmen, Auftrieb, Moment) zu ermitteln braucht. Die Ent- 
wieklung selbst ist eine in Fourierreihen, wo nur die eos na Glieder stehenbleiben, so daß auch die Koeffizienten- 
Ermittlung einfacher wird als bei Birnbaum oder Munk. Als Nebenergebnis erhält man in einfachster Weise auch 
die Gesehwindigkeitsverteilung an den Profilen. Auf anderm Wege kam Glauert (Tragflügel- und Propellertheorie) zu 
ähnlichen Ergebnissen. Das hier angewandte Verfahren ist aber sehr einfach und allgemein anwendbar. Die Arbeit 
konnte deshalb als Ausgang für Doppeldecker-, Gitter- und Gleitflächen-Probleme dienen, die sich hiermit in recht ein 
facher Weise lösen ließen. Hierüber wird besonders berichtet werden. Die Anwendung einer ähnlichen Methode auf 
endlich dieke Profile ist in Vorbereitung. 
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II. Die Isoklinen und Isotachen der Strömung. 7(z)=y-+iyr stellt bekanntlich eine 
Strömung in der z-Ebene dar, wenn man als das Potential, » als die Stromfunktion auffaßt. 


dy(z) nr A ; ; : £ 2 
- ID t+id ist dann der an der y-Achse gespiegelte Wert der Geschwindig- 
dz 
-- 2 20 dyız) en = 
keit der Strömung. Damit wird In Fr In— w=Inw-+i(r vr). Hierbei ist w der Absolut- 
wert der Geschwindigkeit, » die Richtung, den diese Geschwindigkeit mit der «#-Achse 
einschließt. Die Linien In w==const sind Isotachen, die Linien »=const Isoklinen der 


Strömung (2). 
Für die foleenden Betrachtungen ist es zweekmäßieer, 


ls (z 
In — wi 


Ep »+ilnw 


u(z)= ila 


einzuführen. 

Die Strömungsrichtungen sind also hierdurch als Realteil, die Logarithmen der Ge- 
schwindigkeiten als Imaginärteil einer Funktion « der komplexen Veränderlichen 2 dargestellt. 

Für die Parallelströmung in Richtung der positiven «-Achse #= - z um ein Strecken- 
profil ist natürlich „= const—0. Für andere Strömungen ist aber u(z) von z abhängig. 
Bildet man nun die Isoklinen und Isotachen der Strömung in der z-Ebene in die [-Ebene ab, 
so ist das Bild «({) natürlich nicht mehr das Isoklinen-Isotachenfeld der Strömung z in der 
£-Ebene. 

Man kann sich aber in der &-Ebene irgendein im Äußeren des Einheitskreises reguläres 
u(£) geben und nach dem Profil suchen, das in der z-Ebene gelegen ist und a(z) dort als 
dy 
da’ 
ganz allgemein ein Verfahren entwickeln, um zu in der Z-Ebene gegebenem u(£) Profile in 
der z-Ebene zu finden. 

Wir wollen uns hier jedoch nur mit dünnen, wenig gewölbten Profilen abgeben. Solchen 
Profilen können natürlich nur ganz gewisse 1. ({) entsprechen. Zunächst wollen wir noch an- 
nehmen, daß die Profile glatt umströmt werden, d.h. daß Anfangs- und Endpunkt dieser 
Profile die Verzweigungspunkte der Strömung werden. Bei dieser Strömung ist die Riehtung 
der Geschwindigkeit auf beiden Seiten des Profils die gleiche. Je schwächer die Profile 
eekrümmt sind, desto mehr werden aber auch die Potentiale auf beiden Seiten der Profile 
einander gleich, nämlich 7°... Am Bildkreis ist p(k) = cosa. Die Richtung v ist bei ver- 
schwindender Dieke und Wölbung nur von cos «a abhängig, ist also, da cos - a=cosa, für 
positive und negative a die gleiche. 

In Abhängigkeit von a als Parameter muß sich demnach »v in der Form 


Isoklinen-Isotachenfeld besitzt. Es ıst nun allgemein te» Man kann auf diese Weise 


. | .» r > wu: 
vcf v, cosa+tr,cos2a - 1,00s9a-+. 


darstellen lassen, wobei natürlich ebenso wie vr auch »,, #,.”,... nur sehr kleine reelle Größen 
sein können. 


ill. Profile, deren Isoklinen-Isotachen in der Bildkreisebene durch Doppelquellen dargestellt 
werden. Die Darstellung von » zeigt, daß v» der Realteıl von 


! y ’ P, | 
u=(r,tilnwz)+ - +5 -5+. re] 
ist. 


. . 1’ ’ £ > . x .. 
Die einzelnen Glieder —, 3, 3: -:--- stellen jedes eine Strömung aus Doppelquellen 


erster und höherer Ordnung dar, die in ©=0 gelegen sind. Es ist deshalb interessant zu- 
nächst zu untersuchen, welehe Profile sich aus Isoklinen-Isotachenfeldern ergeben, die in der 
„-Ebene durch solche Doppelquellen dargestellt werden. 

Da die v; nur sehr kleine Werte haben sollen, so kann man 


tov Y’ 


setzen. Dann erhält man die gesuchten Profile, entsprechend 


x 
y=|tgvdce+C, z=cosa, dx sınada 
| 
«aus 
\(r ' cosa-tr,cos?2 v. cos» | \sı la+Ü 
VE ’„, — y,c0sarYy,Ccosza ‚„cosvoa-—. ...)/)SsıNaqa- . 
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ne Man kann natürlich y in seine einzelnen Anteile 
Bt. . 
12 Yı F\r,cosnasinada+(, 
0 
". zerlegen. Soll aber y,(—- 1) = y„(+1)=0 sein, so muß man auch noch v, in solehe Anteile 


zerlegen, daß dieser Bedingung genügt wird 


er v, “ "ni + V”o9 "3 + "oa ’ 
Dann kann man schreiben 
Ya Fl/(r,ntr„eosna)sinada+(y,. 
0 
Durch Integration ergibt sich 
(? 
| v, 1) v, BD) 
%. Y,- - > (1 cos” a) E= > (I x”) 
n- 
B. & v, Ban (in - Va cos(n + 1a 2 | ' 357 | 
y r =D 1... 
b, In 2 n — 1 n—+ | n? 1 
Tr 
v„[cos(in—1)a cos(n+1la 2cosa i 
Un | > In ö, 4,6. | 
S ' 2 n 1 n + 1 n? —1 
s Die Anteile von v, sind 
e „md fürns=1357... 
. Ya 
a für 1 =2 46... 
N 
- Mit Hilfe der Beziehung 
4 
m m p} IL 2 | m n) 2 Im A 
X COS MACOS A | 9 ‚a COS” a) COS "@ + )d COS” a)” cos A E. RERAR 
P ” 
ü oder 
cosa=ı cosda= br - Wr?’ +5x 
l cos2a=2r? 1 cos ba = 32 x" 8x’ +-18.r? | 
cosdsa=4ır? —3r cos 7 a=64 x” — 112 x’ +56 x? — 7 x 
cosda=8r-- sr’ +1 cos Sa = 128 a? — 256 ® + 160 x’ - 320° +1 
| läßt sich y(x) darstellen: 
” Yv 16 
R 1 2 A ) m... ) _.,4 . 
| Y, er (1 — x?) Y, ER 2-+5a 3a 5“ 
V, = | | ”. - > u A v 6 
IE » (1 e) 3 d UR .. 3 15.430. . 16 .r ) “ > . 
Yv 


NY — + > ( 1+3x0° —2x°) 


IV. Bestimmung der Koeffizienten der Doppelquellen durch harmonische Analyse. Ist ein 
dünnes, schwach gewölbtes Profil gegeben, so könnte man durch Darstellung dieses Profils 
als ganze rationale Funktion und entsprechende Zerlegung die Koeffizienten »,,r,.... finden. 


F dy er 

Einfacher ist es jedoch, y’ = j,„ 7u bestimmen und dann y' in Abhängigkeit von a = are cos 
N . » 

darzustellen. Die harmonische Analyse, die bei unserem Profil natürlich nur cos n «Glieder 

ergeben kann, liefert dann sehr schnell neben »,, »v,, v,.... auch v,. Dieser Wert ist nämlich 


der Mittelwert der y’ im Intervall =. 


In vielen Fällen wird es genügen, nur », und v, zu bestimmen. Am einfachsten be- 
kommt man bei einem Profil zweiter Ordnung », und r, aus den Anfangs- und Endneigungen. 
Ks ist 

dy, dy. 


‘) 
d x , d Av r | 3) 
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Ks ıst also 


! | 
v(+J1) +3? vr(—1) "vr, t7% 
und hieraus 
v(-1)-»r(+1) > | | 

v, 5 v. zu HrHrl). 

Diesen Werten entspricht 
d | | 
’=Vvo > ar ‚+-r(+1)). 


V. Anstellungsrichtung für glattes Umströmen. Wir haben verlangt, daß das Profil glatt 
umströmt werde. Da die Sehne des Profils in die »-Achse fällt, wird im allgemeinen die 
kichtung der Strömung im Unendliehen, die Anströmriehtung, die glattem Umströmen ent- 
spricht, nieht in diese Achsenrichtung fallen. Diese Richtung finden wir aus 


un, tiunw.)-=+5-5+ 


Für sehr große { wird der Realteil von «, welcher der Strömungsrichtung entspricht, 


yz=y 


Dies ıst also der Anstellwinkel für elattes Umströmen: 


. . "nn » 
ZI, — | (n EP 


Der Imagıinärteil von « ist für sehr große 
In » = In w, 
Die Geschwindiekeit im unendlichen ı", ist also frei wählbar. 


FL 


VI. Geschwindigkeits- und Druckverteilung bei glattem Umströmen. Aus dem Imaginärteil 
von 1 erhält man die Geschwindigkeitsverteilung »", bei glattem Umströmen auf dem Profil: 


I, ® = 
a » N ») 
In vsina- v,sin2a-tr,sında 
N 
P4 
m, 2 Be 
e’i sın ( "» si 2 1 Y3 sin DEZ 
m, 


Aus der Bernoullischen Gleichung ergibt sich die Druckverteilung 


) U ' a u 
| | | | e:bı sine — v3sin2«e+ v3 sın da Tee. 
4 m, 


Da nun voraussetzungszemäß die », nur klein sind, so erhält man entsprechend 


„ Mi Mr ge? J 
e I4 5, Pa t37°*+ e >1+9 
a: .). >) 
m, | u Mr 
lI+-r sina »,sin?da+vr,sında 

m e 

FT. 

N Ip. si „sin?2a+tr sind t 

q ö(v, sına—v,sınza v„snoa—. ee De 


VII. Auftrieb bei glattem Umströmen und auftriebfreier Anströmrichtung. Aus der Druck- 
verteilung läßt sich durch Integration der Auftrieb bestimmen bzw. deı Auftriebsbeiwert. 
Da der Unterdruck auf der einen Profilseite gleich dem Überdruck auf der anderen ist, ist 
der Druekunterschied auf beiden Seiten des Profils 


ly» 
Wi 


I(v,sina - r,sin?2a+r,sinda .+...). 
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Da die Gesamtlänge des Profils !=2 ist, findet man: 


Im 
Durch Einsetzen erhält man 
u 
Ca = 2\er,sina—v,sin@2a+r,sinda—.+...)sinada. 
u 
In diesem Integral verschwinden sämtliche aus v,,r,..... herrührenden Teile, und es ver- 


hleibt als Auftriebsbeiwert: 


Sam», T. 
Zählt man die Anstellung von der ersten Achse des Profils, der Riehtung der auftriebs- 
[reien Anströmung aus, mit Ö, so ist allgemein 


£, =8nrksind. 


Bei einem Streckenprofil ist k = 1, und bei dünnen Profilen schwacher Wölbung ist k ebenfalls 
»>1. Somit findet man 


i v 
sad=—, 
oder, da v, nur klein ist, 
v 
=. 


Sind alle v„n>1)=0, so ist das Profil ein Kreisbogen, der an seinen Endpunkten den 
Winkel v, mit der Sehne einschließt. Die erste Achse des Kreisbogenprofils halbiert also 
diesen Winkel. Dies Resultat stimmt mit den Ergebnissen exakter Berechnung überein. 


VIII. Moment bei glattem Umströmen und Bedingung für festen Druckpunkt. Aus der 
Druckverteilung kann man auch das Moment des Auftriebs berechnen. Zweckmäßig berechnet 
man zunächst das Moment bzw. dessen Beiwert, bezogen auf den Mittelpunkt (0) des Profils 


Durch Einsetzen erhält man 


JT 
Im \o,sina- »,sin2a+tr,sinda -—+...)ecosasinada 
0 


7 
ni | u oe | En 
m=5 [rer sna—r,snda+ v,sında—.+...‚sındada. 
0 


In diesem Integral verschwinden sämtliche mit anderen Koeffizienten als v, behafteten Glieder. 
Man findet somit 
2 | 


“NG N V. T. 
Bekanntlich hüllen die Angriffslinien des Auftriebs eine Parabel ein. Bezogen auf den Brenn- 
punkt dieser Auftriebsparabel ist das Moment des Auftriebs konstant. Wird die Brennweite 
dieser Parabel null, so gehen alle Angriffslinien des Auftriebs durch den Brennpunkt F dieser 
ausgearteten Parabel. Dieser Punkt wird dann zum festen Druckpunkt. 


Ein Streckenprofil hat einen solchen festen Druckpunkt. Er liegt in = 9: 5, bei der 


Strömung von links und rechts. Bei schwacher Wölbung ist natürlich der Druckpunkt im 
allgemeinen nicht mehr fest, wie man schon am Kreisbogenprofil sehen kann. Aber der 
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Bi 
> , Je dünner und schwächer ge- 


wölbt das Profil ist. Für sehr dünne. schwachzewölbte Profile erhält man also ım Grenzfalle 
einen festen Druckpunkt, wenn das Moment bezüglich dieses Punktes © verschwindet. Be- 
züglich dieses Punktes muß also Zu, verschwinden. Es ist 


Brennpunkt der Auftriebsparabel liegt um so näher an .«: 


n R E. 
“ll „mo | „(+ 
Oder mit oben errechneten Werten 
ö 7 JS T 
“Ni E r.yp | JT v, | (iv, v.) . 


Man erhält also einen festen Druckpunkt, wenn 
Y v 


2 1 


wird. Die übrigen v,(n >2) können frei gewählt werden. Sind diese v,(n>2)= 0, so ist 


| 7 | 
Y < (i — =) 3 2 [ri )=+ n ,,; r(+)=1 3 rı 


das entsprechende druckpunktfeste Profil. Seine Anstellung für glatte Umströmung ist 


r ; : . . » 
„ verschieden, so daß die Richtung der auf- 


Die Riehtune der ersten Achse ist von dieser um 


triebfreien Anströmung um 6 von der Protilsehne abweicht. 
1) 


IX. Geschwindigkeitsverteilung bei beliebiger Anströmung. Bei beliebiger Anströmung 
ändert sich die Lage des vorderen Verzweigungspunktes. Es sei d die Anstellung gegenüber 
der Anströmrichtung der glatten Umströmung. Dann wandert am Bildkreis die Lage des 
vorderen Verzweigungspunktes um 25. Lag dieser also unseren Voraussetzungen entsprechend 
in az, so liegt er nach weiterer Anstellung in a=r+26. Die Richtung der Strömung 
ist nun zwischen dem alten (E,) und dem neuen (E,) Staupunkt um 7 von der Richtung der 
glatten Umströmung verschieden. Dem ursprünglichen Isoklinen- und Isotachenbild u(£) der 
Protilströmung in der Bildkreisebene muß also noch ein weiteres u,({) überlagert werden, 
von dem zwei »-Linien, die um z verschiedene Argumente haben, mit dem Bildkreis zu- 
sammenfallen. 


Man sieht leicht ein, daß «s(£) dureh die Potential- und Stromlinien zweier Wirbel 
von der Zirkulation /'r, 2a bzw. IE; = +27 gebildet wird, die ihren Sitz in E, bzw. Ey 


haben. Hieraus läßt sich In  ° der zusätzlichen Strömung bestimmen. Nach einigen Um- 
®, 


formungen findet man hiermit die Geschwindigkeitsverteilung bei der Anstellung ö gegenüber 
der Geschwindigkeitsverteilung bei glatter Umströmung: 





| \l/g 
ms m, vos h) 1 > | te 2) 
oder mit dem für », gefundenen Werte: 
‘ | A 
n | 
2 A + r,sına—vr,sında cosöll+tg ztg 9). 
m, ’ 











Als Vorzeichen der Wurzel wählt man auf der Saugseite +, auf der Druckseite —, ent- 
sprechend der Vermehrung bzw. Verminderung der Geschwindigkeit auf den entsprechenden 
Profilseiten. Als Beiwerte für Auftrieb und Moment findet man hiermit 


Cu =2asnö+tn:v, =2nsin(ö+?v)), 
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Die zweidimensionale Wärmeströmung des Beharrungszu- 
standes im rechteckigen Querschnitt geblätterter Eisenkörper 
bei flächenhaft, unstetig oder stetig verteilten Wärmequellen‘). 


Von Herbert Buchholz in Berlin. 


ie vorliegende Arbeit behandelt einige technische Beispiele für die zweidimensionale 
Wärmeströmung in anisotropen Körpe ın mit rechteckigem Querschnitt. Über die 
in allen Querschnitten gleichartige Verteilung der Wärmequellen wird in Übereinstimmung 
mit den wichtigsten praktise hen Fällen angenommen, daß sie nur nach einer der beiden 
Hauptrichtungen variabel ist, dann aber auch unstetig, während sie nach der andern 
Richtung konstant ist. Als Grundlösung für die verschiedenen Beispiele dient die Beziehung 
für die Temperaturverteilung, die eine einzelne, zu einer der Seitenflächen parallele, lächenhafte 
(Juelle in dem rechteckigen Querschnitt des Körpers erzeugt. Diese Lösung wird mittels kom- 
ple xer Integration hergestellt und aus ihr die verschiedenen für die numerische Berechnung ge- 
eieneten Reihenentwie ee 'n abgeleitet. Zum Schluß wird noc h ein Zahlenbeispiel angegeben. 
- Die Berechnung der stationären Temperaturverteilung in dem geblätterten Eisenkörper 
elektrischer Maschinen führt nach Vornahme einiger naheliegender Vereinfachungen auf eine 
Aufzabe der zweidimensionalen Wärmeströmung in einem rechteckigen (Juersehnitt. Abge- 
sehen von der Anisotropie ist für sie vor allem charakteristisch das gleichzeitige Auftreten 
flächenhafter und stetig verteilter Wärmequellen. 


















































Um die Daten der Aufgabe möglichst deutlich m FR gu) 
vor Augen zu führen, mögen sie an Hand von Abb. 1 
kurz besprochen werden. — 

Das Bild stellt einen Meridianschnitt durch | 
ein mittleres Eisenpaket des Ständers von der — == Kühlschlitz 
Breite 2b dar. Die Richtung der y-Achse im Bilde “ = mit uy 
verläuft dabei parallel zur Maschinenachse. Die as | Akypferstab .a, PR 
‚-parallele Schraffur deutet die Blätterung des | 1% 
Eisens an. Am linken Rande des Querschnitts -@ en E 
auf der dem Luftspalt zwischen Ständer und Läufer %Y | 
zugewandten Seite ist der als nicht unterteilt EmER | Mx2 
vorausgesetzte Kupferstab einer Ständernut einge- Az) L ” 
zeichnet. Die an den vier Kanten stehenden Buch- i_ı te = yühlschlitz 
staben «u bringen die Gleichheit oder Verschiedenheit se = mit 1, 
der an den zugehörigen Flächen geltenden Wärme- | 
abgabeziffern zum Ausdruck. Die Diagramme über ur; 7 e 
und links neben jedem Querschnitt geben die an- a [A025 


genommene Verteilung der Wärmequellen in der 
x“- oder y-Richtung wieder. 

Mit diesen Daten enthält die Aufgabe verschiedene andere, einfachere Probleme als 
Spezialfälle. Unter diesen verdient besondere Erwähnung das Beispiel der Temperaturver- 
teilung im Querschnitt der Kernsäule eines Transformators. Die Vereinfachungen der obigen 
Problemstellung betreffen in diesem Falle einmal die mit konstanter Ergiebigkeit über den 
Querschnitt verteilten Wärmequellen, zum anderen die an allen Kanten praktisch gleichgroße 
Wärmeabfuhr. Die mathematische Formulierung der Aufgabe wird natürlich nur die anfangs 
besprochene, kompliziertere Aufgabe berücksichtigen. Die notwendigen Vereinfachungen 
werden erst hinterher an der zugehörigen Lösung vorgenommen werden. 

In dem allgemeinen Falle darf wegen des in der y-Richtung konstanten Maenetfeldes 
die Ergiebigkeit der Wärmequellen q als von y unabhängig angesehen werden. Im Innern 
eines jeden Querschnitts hat dann die Temperatur 9 (x, 1), die auf die überall als gleich- 
groß vorausgesetzte Temperatur des kühlenden Mittels bezogen werden möge, zunächst die 
Bedingung !) zu erfüllen: 


Abb. 1. 


.. 9° 09 Re) 
tere... nen 


*, Wissenschaftliche Mitteilung der AEG. 

1) Die Aufgabe könnte nach der gleichen Methode auch noch in dem Falle gelöst werden, daß für das Kupfer 
eine endliche, vom Eisen verschiedene Wärmeleitfähigkeit in Rechnung gestellt wird. Die Beziehungen würden dann 
allerdings sehr kompliziert werden. 

Die gleiche Aufgabe wurde unter der Voraussetzung Ay” * als eindimensionales Problem behandelt von 
0. Böhm: Vorausbereehnung der Erwärmung elektrischer Maschinen, ETZ. Bd. 43, 1922, S. 810. — E. Roth: Etude 
analytique de l’öchauffement des machines eleetriques, Bull. de la Soe. france. des Eleetriens, Ser. 4, Vol. 7, 1927, 
pP. 840... 966. 
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Hierin bedeuten 4, und 4, die beiden Wärmeleitfähigkeiten in W/em °C in der «- und y-Richtung. 
An den zur «-Achse parallelen Begrenzungslinien des Querschnitts sind die beiden Wärme.- 
übergangszahlen «, in W/em? € gleichgroß. Die Temperaturverteilung ist dann gewiß zur 
x-Achse symmetrisch, und es muß daher gelten: 


Pi 
u (0) (21<Za,y Be eV. y ©, TE KO ARITENE (la), 
£ 0 () 
Hy ) (X, x b) + Ay h y 0 (ITI<QA, „ + b) P A : : ; \ i : (1b). 


In bezug auf die y-Achse ist dagegen die Temperatur sicherlich unsymmetrisch verteilt. An 
der rechten äußeren Kante muß nämlich die Grenzbedingung (lc) erfüllt sein. 


‚009 
U Ita, y)TAxrı, ed yvzbehe). . ip. Ubi ck ng re 


An der linken Querschnittskante lautet aber die entsprechende Grenzbedingung wegen der 
Anwesenheit des selbst wärmeerzeugenden Kupferstabes ganz anders. Um weitläufige 
Rechnungen zu vermeiden, kann für den Kupferstab wegen seiner nur kleinen Abmessung 
in der «-Richtung und wegen der großen Wärmeleitfähigkeit des Kupfers eine örtlich kon- 
stante Temperatur ©, angenommen werden. Bezeichnet dann @ die Ergiebigkeit der Wärme- 
quellen im Kupfer in W/em® und d in em die Breite des Kupferstabes, so folgt aus dem 
Energieprinzip zunächst die Gl. (2a). Ist andererseits e die gegenüber d sehr geringe Dicke 
der Isolation zwischen Eisen und Kupfer und 4, deren resultierende Wärmeleitfähigkeit 
wiederum in W/em ’C, so besteht nach der Definition von 4 die Gl. (2b): 


00 

() 2 d Ux 1 J, Ihx £ hi ” (za), 
ö . 009 
I, I, — 9(-a,y)l;e=— 4; 3. (gl: b,.r We RD 4 


Kliminiert man aus beiden Formeln die vorläufig nicht interessierende Kupfertemperatur, so 
gilt für den Wärmestrom im Eisen an der dem Luftspalt zugewandten Seite als vierte 
Grenzbedingung die Gl. (ld): 


’ | Ö ie) ’ 
Mxı O(—-a,y)—4;r- "wer 0 .d (yi<b,x Be ar 
Ka Uxıl e Ux;| 
’ A | ; = . ° . . . . . . ° . ° lo r 
Q u Ir | 110) 


Die Isolation zwischen Kupfer und Eisen setzt also einmal für das Eisen die Möglich- 
keit der Wärmeabgabe an die Luft herab, andererseits verringert sie aber auch die Er- 
wärmung des Eisens durch die Wärmequellen im Kupfer, zwei Wirkungen der Isolation, die 
durchaus plausibel sind. 

Zwecks Lösung dieser Gleichungen empfiehlt sich zunächst die Zerlegung der physika- 
lisch einheitlichen Temperatur des (@uerschnitts in zwei Komponenten ©, und ©,. Die 
erste Komponente ©, soll allein von den Wärmequellen im Eisen herrühren. Für diese 
Temperatur werden demnach die Gl. (1... 1d) bestehen müssen mit der einzigen Abänderung, 
daß in (ld) 9-0 zu setzen ist. Die zweite Temperaturkomponente möge dagegen lediglich 
durch die Wärmequellen im Kupfer erzeugt werden, Für ©, gelten mithin die Gl. (l...1d) 
mit 4(#)==0. Die beiden Temperaturfelder unterscheiden sich also im wesentlichen nur da- 
durch, daß in bezug auf das Eisen die Erzeugung des Feldes ©, durch räumlich verteilte 
(Juellen, die des Feldes ©, dagegen durch eine flächenhafte Quelle verursacht wird. 

Die Ausdrücke für die beiden Temperaturfelder ©, und ©, lassen sich aus der Lösung 
einer einzigen neuen Aufgabe herleiten. Dieses Ersatzproblem betrifft die Berechnung 
der Temperaturverteilung 9 in dem gegebenen rechteckigen Querschnitt unter den Bedingungen 
(l...1d), wenn in (1) und (ld) zwar q(&)=0 und Q=0 sind, dafür aber an der beliebigen 
Stelle e=s mit —a<s<-+a eine schichtförmige, flächenhafte Quelle mit der Ergiebigkeit 
I liegt. Zu den schon bestehenden Bedingungsgleichungen (1...1d) für die Temperaturver- 
teilung, die nur entsprechend dem Übergang von © zu d und dem Verschwinden von q(&) 
und @ zu verändern sind, treten dann noch die beiden folgenden hinzu: 


VIs:s V,) a Fe (le), 


2 1.09) a rn a 7 


Q SX—s—(\ Oo ix 
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Das Svstem (1) von Bedingungsgleichungen für die Temperatur 9 (x, y) kann noch be- 
‚uemer formuliert werden durch Benutzung eines anderen Koordinatensystems. Setzt man 


nämlich 
. 4/3 a . 4/2 3 ° 2 2 2 s 
z=6&. YVAhx irre ' YENYhyliıx -» - (3b), Am If CT Fo 2° 
4/53 75 } — mn. . 
Uy=lUy'YtxlAy- - - 8), Aur=ur- VYiylde: » » (80), 


so hat das neue Gleiechungssystem für 9 die folgende, für die Lösung bequemere Gestalt: 


2, 9%Dd 


— + wi , , ie AT I en (4). 
0° 07 

A}, g . 0% ; 

Be er ee er EEE ar 0 (n=+P). (4b), 
04 
OU ö 5. { Vu j 

Uxs'0(0; ra,n)-+ Ay DE V (& ra) (de). u zı’9iG; a,n) Am'yE 0 (E a) (dd), 

s O7, 0 nn 
J)l0:0 ‚) (0,0 t-, 29) Be N: Am’ Er, inne), nr Aula (tf). 


In diesen Gleichungen stehen entsprechend (3a, b) für a,b, s die Zeichen «a, P,6. Die vor- 
enommene Koordinatentransformation hat die Bedeutung einer Streekung oder Schrumpfung 
aller Abmessungen in der x- und y-Richtung im Verhältnis } ;y| 4x oder ] /x/ A, unter 
Wahrung des gleichen Flächeninhalts für den Querschnitt. Die Ergiebigkeit der Quellen 
bleibt dann bei dem Übergang von dem einen zum anderen Koordinatensystem unverändert. 
Die Streekung dagegen verändert sowohl den Zahlenwert der Wärmeabgabeziffern als auch 
den der Wärmeleitfähigkeiten, und zwar diese letzteren absichtlich in dem Sinne, daß beide 
für die «- und y-Richtung des Querschnitts gerade gleich 4, werden. 

Die vollständige Lösung des Gleichungssystems (4) wird aus zwei Funktionen bestehen, 
von denen die eine 9 (0; £&,n) für &< und die andere 9 (o,&:;n) für o<E& gültig ist. Hat 
im Koordinatensystem der &, y die bei #=s gelegene Wärmequelle die Ergiebigkeit 4 (s)- ds, 
so ist offenbar die Lösung des Systems (4) im Koordinatensystem der &,n mit q(s)- ds 
y(o)-do:- Yix/A, zu multiplizieren. | 

Aus der Lösung von (4) kann sofort das Verteilungsgesetz für die vordem eingeführte 
Temperaturkomponente ©, gewonnen werden, denn diese sollte ja selbst von einer flächen- 
haften Wärmequelle herrühren, die nur nicht wie die Quelle 4(s) an einer beliebigen Stelle 
des Querschnitts liegt, sondern speziell an dessen linkem Rande. In diesem Falle lassen 
sich aber die drei Gl. (4d...4f) zu der einzigen neuen Grenzbedingung (4h) zusammenziehen: 


9 
a. | 5 | BR 
Uxı Vla;a—,n) Amy VYAylkx (nl<Pp,$ o).: x»... fh) 
Man braucht also nur die Lösung von (4) für den speziellen Fall = — «a noch mit 9. d zu 


vervielfachen, um gemäß (6) die gesuchte Temperaturkomponente ©, zu erhalten: 
9,(2,y=W-d-d(—a;f,n) . . . Et nd 21 


Wird die Erwärmung wie bei der Komponente 9, durch ein räumliches Quellgebiet 
(x) verursacht, so läßt sich auch die hierzugehörende Temperaturverteilung sofort aus der 
Lösung (4) gewinnen, und zwar durch eine einzige weitere Integration. Berechnet man 
nämlich ©, (£, 7) nach der Vorschrift (5) 
4 


5 t 


4 


wer 
Ar 
9. (,n)= | Pu \ (0: &,n)-qlo)-do+ \ 
/ y R 


| to; &,)-4(0)- do Erg (»), 
hy 
a = 

so erfüllt diese Funktion Jedenfalls die entsprechend umzuschreibenden Bedingungen (4a... 4d) 
und in Rücksicht auf (de, f) statt der Laplaceschen Gl. (4) die Poissonsche Gl. (5a), die 
nach Übergang zum System der (#, y) mit (1) übereinstimmt. Mithin befriedigt die so her- 
gestellte Funktion ©, (x, y) 

E 0’ I, . "Oo, 2 € 

Am’ DE? TAm' on” gl) v ER BE SE, Sue SZ EZ Ze u 


tatsächlich alle Bedingungen der Aufgabe. Den Gl. (5) und (6) zufolge sind also die ge- 
suchten beiden Temperaturen ©, und ©, bekannt, sobald die Lösung des Gleichungssystems (4) 
vorliegt. 






































er. } = 2 ’ 2 u x Zisehr. f. angew 
IS8 3uchholz, Zweidimensionale Wärmeströmung in geblätterten Eisenkörpern Math. nnd Mech. 





Die Lösung von (4) kann am einfachsten mit Hilfe der komplexen Integration?) erzielt 
werden. Als Lösungsansätze mögen die beiden folgenden Ausdrücke dienen: 





) [ id -+B | 3 ned I (£ \.(W@ 
7(0,5,0)= % Alu)-cosucsH+ (u)-sın ugs, "er Sn: ‚au en. va), 
10 u...” ' koiußptuk, Sinuf 
(2L) 
) Er; + Diu)- si > un I (£ \. (A | 
VI0,5,N) e wu)-cosuwz-r Diu) - Ssınu SS,’ zz. : pr + ou e > os. ib), 
ae 77 J s | ' Kojuptuk,  Sinup 
2 L 
f worin k,, gemäß (Sa) definiert ist. (2 .L) bedeutet den 
wr in Abb. 2 dargestellten Integrationsweg der u-Ebene, 
| der aus zwei nach beiden Seiten ins Unendliche 
[72 7, 45 
| "Ze 79E 5 Pe E . ; 9 
290 7% U FF 34° 
K,; ki | Hi . . . . . . (Se) 





auslaufenden und etwa zur x- und y-Achse symme- 
trischen Kurven besteht. Selbstverständlich müssen 
die beiden in (7) vorkommenden Integrale, wenn 
der Ansatz überhaupt einen Sinn haben soll, zu- 
nächst als konvergent vorausgesetzt werden, und 
erst später kann nachgeprüft werden, ob diese Vor- 
iz aussetzung auch tatsächlich zutrifft. Aus demselben 
| Grunde darf keine der geschweiften Klammern eine 
gerade Funktion enthalten, weil sonst das Integral 
identisch verschwände. Die bisher allein angeh- 
baren Pole der beiden Integranden liegen an den Stellen w= +i-r,, wo r, die positiven 
Wurzeln der Gl. a) sind. Uber deren Lage zu den Schnittpunken des Weges (2L) mit der 
y-Achse kann und soll jedoch vorläufig nichts ausgesagt werden. 











Abb. 2. 


ro ky hy - (9: ) Q 
Ceg T, B Tı „ or ee KR), [En . SP ES Er | | ° 


Die beiden Ansätze (7) erfüllen von den sieben Bedingungen, die das Gleichungs- 
system (4) enthält, vorerst auf jeden Fall die Gl. (4) und (4a), die Gl. (4b) aber nur dann, 
falls z. B. für (Ta) die 

l 5 i : u 
= tat) cosus ra mut) . 2. rn 
ni. 

(2L) 
Beziehung (7) bestehen sollte. Drei weitere Bedingungen können aber sehr leicht ohne 
Einschränkung dadurch befriedigt werden, daß man jeweils nach Zusammenfassung aller 
Glieder die geschweifte Klammer gleich null setzt. Damit die so entstehenden Gleichungen, 
die als Bestimmungsgleichungen für die vier unbekannten Funktionen A (u) usw. aufzufassen 
sind, möglichst kurz angeschrieben werden können, mögen schon an dieser Stelle die folgenden 
vier Abkürzungen?) eingeführt werden: 


n„(u,a)=cosua u: k’,.,-sinua. (10a), n.,(u,a)=sinua+u-k’,,-cosua . (10e), 
N".,(t,a) =cosua  u-k,,:Ssinua . (10b), n,,(u,a) =sinua+ u:k,,-cosua . (10d). 


Die Berücksichtigung der drei Bedingungen (4e...e) führt mit Benutzung dieser Abkürzungen 
zu den folgenden drei Gleichungen: 


a) a.) Ban) a. a ded, et re er ae 
Cl) m.) — Die) nm. wm. d=0 . 28 Era are a 
Alu)-cosuo+ B(u)-siınuo (C(u)cosuwo— D{u)snuo=V. . . . (le). 


*), Über die erstmalige Anwendung dieser Methode vgl. man 1. T. Bromwich: Proe. London Math. Soe., 
2 Ser., Vol. %0, 1930, p. 165 .... 177. 

3) Von den beiden Zeigern des Buchstabens » stimmt immer der erste mit dem Zeiger von kx überein, der 
zweite charakterisiert die Funktion als gerade (0) oder ungerade (1). 
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I'm auch die letzte der Grenzbedingungen zu erfassen, müßte nach (4N) und (7a, b) für jedes 
;, die Beziehung gelten: 


YEPR FE. [ u: du-Kof u 


hm 2ai JKojup+u-k, Sinup 
(2L) 


I A(u) sin uo + B(u)cos uo + C(u)sinuo - D(u)cos uo!. 


Setzt man hierin die geschweifte Klammer versuchsweise gleich const | u”, so entsteht zunächst 


Vizji,_ const f Col u du 


Am Zni J Gojuptu:.k, Sinup u 
(2L) 


Da in diesem Ausdruck der Integrand für alle ||) </ und alle Re (x) 0 absolut konver- 
viert, so darf in Abb. 2 der Weg (2L) etwa durch die beiden die “-Achse schneidenden 
Kreisbögen von sehr großem Radius / zu einem geschlossenen Kurvenzug ergänzt werden, 
ohne daß sich für R>»x der Wert des obenstehenden Integrals verändert. Wird dann ein 
für allemal die Verabredung getroffen, daß abgesehen von «0 die Pole + i-r, des obigen 
Integranden alle oberhalb des oberen und unterhalb des unteren Wegteils von (21) liegen, 
mit anderen Worten der Weg (21) die y-Achse zwischen den beiden kleinsten Polen +: r, 
schneidet, so läßt sich nach dem Residuensatz die obige Bedingung erfüllen, wenn const 

Viuldz 

Äys 

Funktionen A (uw) usw. noch die weitere, vierte Gl. (11d) besteht: 


gesetzt wird. Die Forderung (4f) ist demnach befriedigt, falls zwischen den vier 


4 - En 

Pr 

A(u)-sın eo + Blu) - cos wo + ( (u) -sın wo D(u)-cosuo n; — . % 18) 
"m 


Die Gl. (Ila...d) reichen hin, die vier unbekannten Funktionen A... D zu berechnen. 
Setzt man die errechneten Ausdrücke sofort in den Ansatz (7a, b) ein, so erhält man als 
Integraldarstellung für die Temperaturverteilung infolge der bei £=6 gelegenen 
lächenhaften Wärme«uelle die beiden Beziehungen: 


3 'n,,(u,a—o):n,,(u,a+%£) Gofan n 
"(0:5,%) ne B ä | er 7 ‚du (Z<o) (l2a), 
5 ni. u” A, (u) bojupßtuk,- Sinup 
(21) 
Io: &,.n) Fat, 3 fer I). n,,(u.a-+ 0) Go] us, In (& \ (1%bı 
Vo: £.ı = =. ze = .—— = er . du (E>o >h). 
An ni u” A, (u) bojup+uk,. cntup 
(2L) 


worin zur Abkürzung die ungerade Funktion A, (uw) gemäß (12) bestimmt ist: 


IA, en,.lu, a). nn, (u,a)+tn,.,(u,a).n,,(t,«) 


.) 
— (1 Krk. WW) -sin2ua+(k’., +k,,) uw-cos2ua . 
Die beiden Lösungsintegrale (12a, b) des Gleichungssystems (4) sind tatsächlich ent- 
sprechend der eingangs erwähnten Voraussetzung auf dem gewählten Integrationsweg absolut 
konvergent. Ebenso ist auch die Bedingung (7) erfüllt, denn die beiden Wegbestandteile von 
(21) können in diesem Fall, wie in Abb. 2 angedeutet ist, nach oben und unten durch zwei 
Kreisbögen von sehr großem Radius zu geschlossenen Wegen ergänzt werden, in deren Innern 
der zugehörige Integrand singularitätenfrei ist. Schließlich befriedigen die Lösungen (12a, b) 
auch die Differentialgleichung (4), denn durch das Einsetzen von (12a, b) in (4) entstehen 
identisch verschwindende Integranden. Die Gl. (12) stellen demnach tatsächlich die richtige 
lösung des durch das Gleichungssystem (4) definierten Temperaturfeldes dar. Nach einer 
bekannten Beweisführung ist sie dann zugleich auch die einzige Lösung. 


Aus der Integraldarstellung (12) kann sehr leicht eine für die numerische Berechnung 
geeignete Reihenentwicklung gewonnen werden. Die Möglichkeit dazu bietet der Residuen- 
satz. Seine Benutzung hätte in dem allgemeinen Falle von (12), wo die Temperaturquelle an 
einer beliebigen Stelle o liegt, in ganz derselben Weise zu erfolgen wie in den gleich zu be- 
handelnden speziellen Fällen. Aus diesem Grunde soll auf die Herleitung der Reihen- 
19 
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entwieklungen aus (12) für eine beliebig gelegene Temperaturquelle verzichtet und diese Be- 
ziehung sofort auf die Bestimmung der Temperaturkomponente 6, (x, y) angewendet 
werden. In Rücksicht auf (6) und (12b) ergibt sich zunächst wegen n,, (u,0) = uk',, 


V.d 1 [" (u,a—£) Col un 


At )E Euer * — 7 ar 
J xı 2ni J u-4,(W) boluptu:-k,- Sinup 
(?L) 


Bu 


[Ei 


Um von (13) aus zu einer geeigneten Reihenentwicklung zu gelangen, ist nur erforderlich, 
den Intererationswer in (13) auf die reelle Achse zusammenzuziehen und die Residuen für die 
dort lieeenden Pole zu bereehnen. Für den Pol «=0 ist das Residuum abgesehen von dem 
ersten Faktor vor dem Integralzeichen von (13) durch den Ausdruck (14) gegeben: 


7 Er x Ad “— Rh. vr 

0 y kx; 7 h. 2 2a + hx; T h. 2 (14), 

hr’, N A, / np r ä . ° . z . ° (1 tat). Bi ER | Ay . R ° . A R A . c (14y), 
PR » | 4x = ; . . P F ö : ° (149), h,=4y| 1, ! £ j : i . { y (A 16). 


Der Integrand von (13) besitzt aber außerdem noch eine unendliche Folge von rein reellen 
Polen an den Stellen u = + o, mit p=1,2,3...., deren Lage die Gl. (9b) festlegt. Aus der 
bekannten graphischen Darstellung beider Gleichungsseiten zur angenäherten Bestimmung 
der Wurzeln ergibt sich für die Lage der Wurzeln die Ungleichung (p — 1)-r/2<0o,<p:-a[/2 
mit der genaueren Angabe 0,>(p— D):-r/2 für py>x. Andererseits besteht 


’ 
h’x; hr: 
Ad Ad | | 
ctg 20, = 0p° 1 ea u Ga. :...:+ &) 
WM, „ 8 O5 By, ‚Pi 
a da a Ad 


für die erste Ableitung von A,(v#) nach « an einer Nullstelle die Beziehung (15): 


u Ba a A in > sin ?o, 
‚ | x x2 / ) x x9 > / - 
iu. a+ N (0,) a=| r - + 3.11 . 0, |° = (15). 
1 / ) | y .) - I v ‘) 
\a a /sınZzo, a a 0 





lın ganzen erhält man damit aus (13) für die Tremperaturverteilung 9,(&#, ), die durch eine 
am linken Rande gelegene Wärmequelle in dem rechteckigen Querschnitt verursacht wird, 
die folgende Reihendarstellung (13a): 


| E.,; h Y ü i 1 [1 h a 2 | N) 
Pe sın 0 + 9"+0,'C0S 0 | 
):.d a a ” a, a ” ’ a, 
Ir, My) ‚ 2. N {a 
x i .) h 2. hx M —— On (0) 
_ ö ) ] 
7 m ’ 
. . (la). 
6 vv b u | 
DD | :G 
h a iu 
h / h h / 
£ x M - m ir - 
bo] d_ 5.) - „om a O) 
dA Ay dl \ A Ay 





Geht in (13a) b>x oder A,>0, so verschwindet in beiden Fällen der zweite Bruch jedes 
Summengliedes. Die erste Zeile von (13a) gibt also die Temperaturverteilung an, die durch 
eine flächenhafte Wärmequelle am linken Rande des rechteckigen Querschnitts erzeugt wird, 
wenn entweder die Höhe des Rechtecks parallel zu dieser Kante unendlich groß ist oder aber 
die Wärmeleitfähigkeit in dieser Richtung verschwindet. 


Die Berechnung der Temperaturkomponente ©,(#, y), die nach den früheren Aus- 
führungen von den Wärmequellen im Eisen selbst herrühren sollte, möge unter der speziellen 
Annahme erfolgen, daß die Ergiebigkeit q(6) der Wärmequellen in (5) im ganzen Querschnitt 
konstant gleich 4 ist. Geht man mit den Ausdrücken (12) in (5) ein, so ergibt sich zunächst 
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| | u I) I n.w,a+gH,+n,.,(u,a+g-[1—n,.lu, a —&)] 


I NE 
‚Ss, 9) u?» A, (u) 


iu SE. 
(2 L) 
Co] um 


Gojup-+uk, Sinup 


Aus dem Zähler des ersten Bruches im Integranden läßt sich jedoch in Rücksicht auf (10) 
und (12) das Glied A,(w«) herausziehen. Für die Temperaturverteilung 9,(£, 7), die von einem 
gleichmäßigen Wärmequellenbelag des rechteckigen Querschnitts herrührt, gelangt man danach 
zu der folgenden Integraldarstellung: 


7 | Kugler D)+n,u,a+5) ‚\ Cojun che (17) 
. . i). 


Am Zniı J | AI,iu) |Kojup+uk, Sinup uw 
(21) 


I (u, M) 


Die Umwandlung von (17) in eine für die numerische Berechnung geeignete Reihe kann 
auf die gleiche Weise vorgenommen werden wie beim Integral (13). Auch in (17) ist der 
Nullpunkt nur ein Pol erster Ordnung. Sein Residuum ist mit Einschluß des Faktors q/A 
gleich der ersten Zeile von (17a). Über die Lage der übrigen Pole gilt das gleiche wie 
früher. Alles in allem entsteht damit für die Temperaturverteilung ©, (xy) innerhalb des 
gleichmäßig beheizten rechteckigen Querschnitts die folgende Reihendarstellung: 


u a1 Pape Ti - Fe h..\ 





are ge \? a? l a a 
I (2, y)= 4, 3 | 2» ‚ 
| Sr a °F Wo. , Aa 
4 7 
a a 
N i (lva), 
: y b RE 
4 u e & Cof| f z | - ö,| 
\’n,(0,a—&)+n,(ona+ 8) b a Ay 
v y. — 3 r —_— a ” 
, 02 * N (0,) x h 2, h Zn aa: I 
== F r Soil / 0,)+ "5, Sin| | he 5, 
a Ay a a Ay 


worin im vorliegenden Falle gemäß (10e,d) für die beiden Zählerfunktionen in der unend- 
lichen Reihe geschrieben werden kann: 


) = sin 0 i + u... B ; PR. 17a) 
n.,(i0,. a )=sın 0 )- —.0,.C0OSo | . ; A , . Ilva), 
1 (9 p a da p p a 

| ER x Pe 2 2 
n,,(0,@ +5) sın o, 1 + - H- + 5. cos0,[1+ .-) N, ER (Ep). 


Die erste Zeile von (17a) allein bedeutet wie vorher die Temperaturverteilung, die sich 
bei örtlich konstanter Ergiebigkeit der Quellen einstellen würde, falls der Körper entweder 
in der Richtung der »-Achse unendlich lang wäre oder in dieser Richtung keine Wärne- 
leitfähigkeit besäße. Die in Wirklichkeit von null verschiedene Wärmeleitfähigkeit A, ver- 
ringert natürlich überall die auftretende Temperatur. Tatsächlich ist auch in den Reihen (13a) 
und (17a) das erste und größte Reihenglied negativ, da o, kleiner ist als =/2. 

Von besonderem Interesse sind die Beziehungen für den Wärmefluß, der in den 
behandelten beiden Fällen aus dem linken Rande des Querschnitts, d.h. aus dem Sitz der 
Flächenquelle des Temperaturfeldes ©,, aus dem Eisen ausströmt. Da an dieser Stelle die 
beiden Wärmeflüsse im allgemeinen entgegengesetzte Richtung haben, ist es für manche 
Fragestellungen der Praxis von Interesse, die Richtung des resultierenden Wärmeflusses zu 
kennen. Im Falle der flächenhaften Wärmequelle ist nach (13a) der parallel der .-Achse 
aus dem Eisen austretende Wärmefluß 





Rn 
09, l h., 1 +2 \’cos20, - (hx,/a)-0,-sin2o, 
‚> * y , 5 4 u‘ .6— e 
" A) x ji da fe) nn h xı wu hx A N (0,,) 
nu a a ; 
Co] 4 ‚ % ae 5.) 
a Ay | 
2075 ER ET 2 
r Pi — Zn ai BR rn 
GCoil--1/ 7-0 op: (h„la)- Zn I 2.8 
j Aa / 2, p + P ( m| ) = a hy ID 
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und nach (17a) derselbe Fluß im Falle gleichmäßig verteilter Quellen 





B.; Mi r 
/ 0m, +10 | A dl > \yl cos 20,+(hx.|a) - 0,:sin 20, 
\ . N ° . Z - - 
dx I u B.; Rn — 3° + N (6,) 
a "T” T Dy= 
Ad Ad 
(1Sb). 
nn h > 
Cof| J wi 5, 
I A Ay 
= - (“ «) 
ei en = 
Goil ,) Ho, (hula)- Zin| | x 5, 
ay Äi, ay ih, 


Nach Gl. (2a) ist außerdem die um Q-d vermehrte algebraische Summe beider Wärmeflüsse 
gleich dem Produkt aus der Wärmeabgabeziffer «,, und der Temperatur ©, im Kupfer. Die 
Berechnung der beiden Wärmeflüsse (15a, b) liefert also zugleich die Temperatur des 
Kupferstabes. 

Die Gl. (13) und (17) nebst ihren Abwandlungen vereinfachen sich wesentlich, wenn 
auch für die zur “„-Achse senkrechten Flächen die Wärmeabgabeziffern einander gleich 
werden. In Rücksicht auf das wiehtige Beispiel des Transformators, bei dem es sich 
obendrein nur um gleichmäßig verteilte Quellen handelt, brauchen deshalb die verschiedenen 
dadurch bewirkten Änderungen nur an der Gl. (17) weiter verfolgt zu werden. Für k,, kr, 
ergibt sich sofort aus (10) und (12) 


n.,(na Dthn,(u,a+! = 2-.cosuf-(snwa—tu:k,-cosua) 


Rz 


| (v)=2-(cosuwa- uk,-sinwa)- (sinwatnu:k,-cosua). 


Für die Temperaturverteilung ım rechteekigen (Wuerschnitt eines Trans- 
[ormators besteht mithin die folgende Integraldarstellung: 


| 


(| An) \ cosu: \ Col u, du 
lcosuwa u-k.-sinua I\Kojup-+u.k, Snup u 


(I) (.r., ) I ; (19). 
eh) 


Ihr entspricht, wenn gemäß (20a) + o, die Wurzeln des ersten Nenners in (19) sind, die 
keihendarstellung') (20). 


h, n A: . . \ 
et oO, ee Öj a a 
( 
r COS = 0 | 
PPARRER ERBE AS BE 2 | in 
i2. 4 | Fk = —. . 
2 Äx | a Ad — On .- + (a/h)+ (hx| d)-0,°% COS 6, 
[7 
20). 
U 5 en 
bo] | . ne ö,) 
ba Ay 
Gof h m wur J Si Aa Be 
) | 9 0, (n.la ‚Sin| za. g 
| A ki 7 n m l A In 





Sie kann aus (19) in ganz derselben Weise gewonnen werden wie z. B. (13a) aus (13). 
Selbstverständlieh gibt die Reihe (20) z. B. auch die Temperaturverteilung im Innern einer 
unsymmetrisch aufgebauten stromdurehflossenen Spule rechteckigen Querschnitts an. 

Außer den schon früher berücksichtigten Spezialisierungen D>x und /,>0 kann im 
Anschluß an (20) noch der Fall «u, >x mit in Betracht gezogen werden. Physikalisch be- 


°%) Die gleiche Aufgabe behandelt auch Dr. E. Roth in seiner in Fußnote !) zitierten Arbeit auf S. ®1 ... m». 
Er lost sie jedoch anf einem mehr elementaren Wege mit Hilfe doppelter Fourierreihen. Die im vorstehenden 
benutzte Methode besitzt den Vorzug, daß sie von vornherein nur zu einer einfach unendliehen Reihe führt. Zwar 
ließe sieh aueh die Lösung von E. Roth dureh geschickte Zusammenfassung der Glieder der Doppelreihe auf eine 
einfach unendliche Reihe reduzieren, aber die große Schwierigkeit besteht eben darin, diese Summationsmögliechkeit 
anch wirklich zu erkennen. 

Bei der Durehreehnung des von E. Roth auf 8. 926 ff. erörterten Beispiels mittels der Gl. (90) der vorliegenden 
\rbeit genügt zwecks Ermittlung der Höchsttemperatur bereits die Berücksichtigung des ersten Reihengliedes, um zu 
dem gleichen numerischen Wert zu gelangen. 
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deutet der Grenzübergang u,>®, daß an den zur „-Achse senkrechten Flächen wegen der 
daselbst äußerst kräftigen Kühlung die Temperatur gleich der des kühlenden Mittels bleibt. 
Mit u. > geht h.>0, und zwar wegen (20a) so, daß cos o,/|h.> (sin o„[a): 0, 2n—|]) 
ı/2-(-- 1)" 1.1/a geht. Aus dieser Bedingung läßt sieh die für «.>x aus (20) hervor- 
eehende Grenzform leicht ermitteln. 


Wiehtiger ist der andere Spezialfall, bei dem etwa die rechte Vertikalkante des (Quer- 
sehnitts mehr und mehr nach außen rückt, so daß schließlich aus dem Rechteck ein einseitie 
ins Unendliche sieh erstreekender Streifen der Breite 2b wird. Um die zugehörige Temperatur- 
verteilung zu finden, muß in (19) vor dem Grenzübergang «>» zunächst £ — «a statt & gesetzt 
werden. Die linke, von der Verschiebung nicht betroffene Rechteekskante fällt dann in die 
‚Achse des neuen Koordinatensystems. Wird das Integral in (19) entsprechend der darin 
vorkommenden Differenz in zwei Teile zerlegt, so kann das allein von », abhängende Integral 
in der sehon mehrfach beschriebenen Weise durch Zusammenziehen auf die reelle Achse 
sofort ausgewertet werden. In dem noch übrigbleibenden Integral werde der obere und 
untere Teil des Integrationsweges für sich mit Z, und 7, bezeichnet. Da längs /,, bzw. I, 
stets /m(u) =O ist, so läßt sich für den ersten Bruch des zugehörigen Integranden im Falle 
+ >» entsprechend schreiben: 


| | | 
.cosu&.eruwieaz —..si E,pFu-i-a 
) COS ME e + >; sın ME e u 
| | rm 
ä .oru.i.a Fo;hr erute 1 u ik, 


Setzt man daraufhin — « statt « in dem Integral längs /,,, so entsteht aus Z, gemäß Abb. 2 
sowohl in bezug auf den Weg als auch hinsichtlich des Integranden gerade das Integral Z,.. 
Im ganzen gewinnt man demnach für die 'Temperaturverteilung in dem einseitig ins Unend- 
liche gehenden Streifen der Breite 2b, der von der örtlich konstanten Wärmequelle y beheizt 
wird, die Integraldarstellung: 


(21). 


bh? \e . .) ) 2 y Uries af 
I (ar, y) > ur Er rs BuRLT TE =] Aulsz.\. e bolun du) 


"FE b = h ib" ih, Kojuptu:k,Siınup uf 
(Lo) - 


Das in (21) auftretende Integral hat völlig die Form der bekannten Hakenintegrale. Es ließe 
sich wie in den früheren Fällen etwa dadurch auswerten, daß der Integrationsweg von links 
und rechts her an die positiv imaginäre Achse herangezogen wird und auf die so entstehende 
Schleife der Residuensatz angewendet wird. Auf die Wiedergabe der sieh dadurch ergebenden 
Reihenentwicklungen soll jedoch verzichtet werden. Auch ohne Reehnung läßt sich aber 
aus (21) erkennen, daß die ersten drei Glieder in der geschweiften Klammer der 'Temperatur- 
verteilung in dem weit entfernten Teil des Streifens angeben, da für £>x das Integral 
verschwindet. 

Mittels der Gl. (5) und (12a,b) kann selbstverständlich auch der Fall eines mit x ver- 
änderlichen 4 behandelt werden, sobald das analytische Gesetz bekannt ıst, das die Ergiebig- 
keit der Wärmequellen in ihrer Abhängigkeit von der “-Koordinate festlegt. Der Verlauf 
von q(x) dürfte sogar sprunghafte Änderungen aufweisen. 

Alle bisher aufgestellten Reihenentwicklungen sind den Gl. (13a), (17a) und (20) zufolge 
dadureh ausgezeichnet, daß ihre erste Zeile jeweils nur die .-Koordinate enthält. Dieser 
Bevorzugung der x-Riehtung entspricht physikalisch der Umstand, daß wegen 4,< 4x die 
Wärme hauptsächlich in Richtung der x-Achse abströmt. Aus der Integraldarstellung der 
angeführten Gleichungen kann aber jederzeit auch eine Reihenentwicklung hergeleitet 
werden, die in dem angegebenen Sinne umgekehrt auch die y-Riehtung bevorzugt erscheinen 
läßt. Zu diesem Zwecke wäre z. B. bei der Gl. (13) in der Weise zu verfahren, daß der 
obere und untere Wegbestandteil von (27) mit seinen beiden Enden nach oben und unten 
an die imaginäre Zahlenachse herangeklappt wird. Die Integrale über die so entstehenden 
beiden Schleifen, die man sich im Unendlichen geschlossen denken kann, lassen sich wieder- 
um durch Anwendung des Residuensatzes auswerten, indem nunmehr die Residuen der durch 
die Gl. (9a) festgelegten Pole berechnet werden. 


Die Verwendung dieser Reihenentwicklungen für die numerische Berechnung der 
Temperaturverteilung wäre aber bei dem nun einmal angenommenen Größenverhältnis 
zwischen 4, und 4, äußerst unzweekmäßig, da ihre Konvergenz weitaus schlechter ist als 
diejenige der im Text angegebenen Reihen. Falls nämlich wie im vorliegenden Falle 4. > 4, 
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ist und entweder 5b größer ist als a oder doch nur wenig kleiner, genügt bei diesen Reihen 
ın der Tat allein die Berücksichtigung des ersten, höchstens aber des zweiten Gliedes, um 
eine für praktische Bedürfnisse ausreichende Genauigkeit zu erzielen. Die Güte der Reihen- 
konvergenz ist natürlich nicht für alle Punkte des (@Querschnitts die gleiche. Am besten ist 
sie wegen 4-0 für Punkte auf der x-Achse, am wenigsten gut für die x-parallelen Ränder 
des Querscehnitts. Aber auch hier führt schon die Hinzunahme des zweiten Reihengliedes 
zu einer ausreichenden Genauigkeit. 


Für ein Zahlenbeispiel zu dem zuerst behandelten Fall (le) mögen die Zahlenan- 
gaben benutzt werden: 


2a=8cm, 2b=-6cem, ix:=0,32 Wlem’C, /„=0,01 Wlem’C, 
iin. IH mW/em?’’C, A, =8mW/cm?’C, 4,=20mW/em?’C, 
/,»=0,SmW/em’C 


, e—(1S5em. 


Nach Gl. (1) berechnet sıch daraus eine infolge der Isolation zwischen Kupfer und Eisen 


reduzierte Wärmeabgabeziffer u'., = 432 mW/em?’C. Da damit infolge der Isolation 
Mx, < 4x, geworden ist, so wird die Höchsttemperatur im 
in 7 | Kisenkörper nach dem linken Rande zu verschoben sein. 
u 128 Für die angegebenen Zahlenwerte lauten die ersten 
L--T TI | beiden Wurzeln der Gl. (9b): 
| 6 
a———— 5, = 0,2732, 0,— 1,6183. 


In Abb. 3 sind mit diesen Zahlenwerten für den Längs- 
schnitt „=0 in Abhängigkeit von x die Beträge der beiden 
eckigen Klammern in den Gl. (13a) und (17a) aufgetragen. 

72 Die Kurven I geben die Temperaturverteilung nach (13a) an, 
wenn nur eine einzelne Flächenquelle am linken Rande des 
(Juerschnitts vorhanden ist, während die Kurven II dem anderen 
Falle der Gl. (17a) entsprechen, wo die Wärmequellen gleich- 
5 45 3 mäßig über den Querschnitt verteilt sind. Die beiden ge- 
IT Tre strichelten Kurven in dem Bilde zeigen den Verlauf jeweils 
| l42 der ersten Zeile der erwähnten Klammerglieder allein. Sıe 
bilden mithin ein Maß für die im Querschnitt sich einstellende 
SE 7 nr” Temperatur, falls keine Wärme nach den Seiten in Richtung 
| der y-Achse abströmen könnte. Die voll ausgezogenen Kurven 
Abb. 3. berücksichtigen dagegen die Temperaturverminderung, die 
durch diesen seitlichen Wärmestrom bewirkt wird. Bei der 
Genauigkeit der graphischen Darstellung genügt in beiden Fällen bereits die Berücksichtigung 
allein des ersten Reihengliedes. 





| 
| 
| 
| 
| 


Nach der Kurve Il liegt tatsächlich aus dem schon eben erwähnten Grunde die Maximal- 
temperatur etwas links von der Mittellinie. Der Unterschied der beiden Randtemperaturen 
beträgt in diesem Falle etwa 4,5°/,, und die Maximaltemperatur liegt gar nur um 1,9°/o über 
der höheren der beiden Randtemperaturen. Die Kurve I dagegen ergibt zwischen der Höchst- 
temperatur am linken Rande und der niedrigsten Temperatur am rechten Rande einen Unter- 
schied von 26,4 "/, bezogen auf die Höchsttemperatur. 


Für die beiden Wärmeflüsse, die am linken Rande aus dem Eisen ausströmen, erhält 
man für das gewählte Beispiel 


im Falle der flächenhaften Wärmequelle 0,212.0-.d, 


im Falle der gleichmäßig verteilten Quellen +0,347-2a:gq. 


Die Wärme fließt also vom Eisen zum Kupferstab, sobald 2a -q > 0,61: @-d ist. Die Temperatur 
des Kupferstabes erfährt unter dieser Bedingung durch den Wärmestrom aus dem Eisen eine 
zusätzliche Steigerung. 425 
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Die maximale Feldstärke und die Kapazität 
der Drehstromkabel. 


Von J. Friedrich Böttcher in Dresden. 
A. Einleitung. 


isher liegen für die Bestimmung der maximalen Feldstärke in Mehrleiterkabeln nur wenige 
B rechnerische Unterlagen vor, die den eharakteristischen Bedingungen des Problems an- 
gepaßt sind. Die bekannte Formel für zwei parallele Zylinder bzw. für Zylinder und Ebene, 
die wohl vielfach unter Vereinfachung der Aufgabe angewendet wird, kann nur sehr grob an- 
eenäherte Werte liefern, besonders wenn sie zur Berechnung der Feldstärke zwischen den 
leitern dient. Wegen der mathematischen Schwierigkeiten, die sich einer allgemeinen Dar- 
stellung der Feldstärke als Funktion der geometrischen Abmessungen des Kabelquerscehnittes 
entgegenstellen, wurden mehrfach experimentelle Ermittlungen an elektrolytischen Modellen 
ausgeführt, u. a. von R. W. Atkinson'), der das Wesentliche seiner Meßergebnisse in einer 
empirischen Formel zusammenfaßte. Verhältnismäßig leichter läßt sich die Kapazität der 
Mehrleiterkabel angenähert berechnen. Werden die Ladungen als in den Leitermitten kon- 
zentriert betrachtet, so ergeben sich auch bei Berücksichtigung des Bleimantels einfache 
Gleichungen ?). 

Man wird also wenigstens in der Frage nach der Feldstärke die Vorteile einer exakten 
Berechnung zugeben können. Die vorliegende Arbeit behandelt nun als wichtigsten Fall das 
Drehstromkabel (Abb. 1), und zwar nach einer von E. Trefftz angegebenen Methode, die 
in etwas anderer Form auch von R. Liebold in seiner Dissertation’) verwendet wird. 
Hierbei werden wir uns auf die Höchstwerte der elektrischen Beanspruchung beschränken. 
Für die Feldstärke stellt sich sofort ein Minimalproblem ein: Die Abmessungen des Kabels 
so zu bestimmen, daß die größte im normalen Betrieb auftretende Feldstärke möglichst niedrig 
bleibt. Bei der Lösung dieser Aufgabe ergeben sich die Kapazitäten. Es wird sich hierbei 
zeigen, daß die nach dem bekannten Näherungsverfahren berechneten Teilkapazitäten erheh- 
liche Abweichungen von den exakten Werten aufweisen können. 

Voraussetzung für die folgenden Ableitungen ist, daß das Isoliermaterial des Kabels ın 
bezug auf seine Dielektrizitätskonstante durchweg homogen ist. Eine weitere notwendige 
Einschränkung ist die Idealisierung der Leiter durch glatte Kreiszylinder, wodurch der Um- 
stand vernachlässigt wird, daß die Leiter aus einzelnen Drähten zusammengesetzt sind. 


B. Hauptteil. 
I. Allgemeiner Teil. 


Bekanntlich ist die Feldstärke definiert als der in Richtung des größten Gefälles ge- 
bildete negative Differentialquotient des Potentials 9. Nun sind die beim Verseilen der Kabel 
üblichen Schlaglängen gegen die Querschnittsdimensionen so groß, daß wir praktisch mit dem 
Verschwinden der axialen Feldkomponente rechnen können. Dann werden also die Verhält- 
nisse durch Betrachtungen in einer zur Kabelachse senkrechten Ebene vollkommen festgelegt. 
Insbesondere wird 9 im Quersehnitt als diejenige ebene Potentialfunktion bestimmt, die im 
Innern des Gebietes, das vom Isolierstoff erfüllt ist, der Differentialgleichung 


a ee ee 


genügt und an den Rändern die vorgeschriebenen Werte annimmt, d.h. also, daß g am Blei- 
mantelkreis verschwinden, an den Leiterkreisen in den jeweils durch den zeitlichen Spannungs- 
verlauf des Drehstroms gegebenen Potentialwert übergehen soll. 


Die Erledigung der gestellten Aufgabe wird dadurch erheblich gefördert, daß der Ansatz 
zur Lösung des Potentialproblems bereits fertig vorliegt in der erwähnten Arbeit von R. Lie- 
bold, die sich die zeichnerische Darstellung der Potentiallinien eines Drehstromkabels mit 
gegebenen Abmessungen zur Aufgabe stellt. Da die hierin enthaltene Potentialfunktion «den 
Ausgangspunkt für die folgenden Ableitungen bildet, soll im nächsten Abschnitt kurz auf ihre 
Entwicklung eingegangen werden. Zugleich werden wir den Ansatz um die höheren Potenzen 
erweitern, um auch für größere Leiterdurchmesser eine hinreichende Genauigkeit zu erzielen. 


I) R. W. Atkinson: Proc. A. 1.E.E. Bd. 38, 1919, S. 815 (s. a. ETZ 1921, S. 108). 
2), L,. Liehtenstein: ETZ 1904, S. 106. 
3) R. Liebold: Diss. Dresden 1931. 
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I. Die Reihenentwicklung der zu berechnenden Potentialfunktion.e Zu Berinn wird die 
Symmetrie des betrachteten Kabels dazu benutzt, das Problem auf den Spezialfall zu redu: 
zieren, bei welchem 


der Leiter 1 das Potential -1\ 

alle übrigen Ränder das Potential 0/J 
besitzen; denn es ist klar, daß an irgendeiner Stelle des Feldes erstens das vom Leiter | 
herrührende Potential proportional dem Potential des Leiters 1 ist, und daß sich zweitens 
ım allgemeinen Fall das Potential aus den von den drei Leitern stammenden Anteilen additiv 
zusammensetzt. 


Kine Vereinfachung der Lösung wird durch die folgende Abbildung erzielt. Liegt der 


Nullpunkt der Ebene 2 «—+iy im Kreismittelpunkt (Abb. 1), so wird dureh die Funktion 
® R - 2 

- a ») 

\ ih R Ze Se IE SS Be 12 SE Zr Ze 2 (3) 


die z-Ebene derart auf die Ebene C=&-+in abgebildet, daß das Innere des Kabels in die 
obere, das Außere in die untere Halbebene übergeht, während der Bleimantelkreis als 
Trennungslinie beider Bereiche der &-Achse entspricht (Abb. 2). 


Y z- fbene 
?R 


ei 











7 2 
Ye 
ie 
I 4 = 
wur 6 5 
Abb. 1. Abb. 2. 


Da die gewählte Abbildungsfunktion der Gruppe der linearen Transformationen angehört, 

so werden Kreise wieder in Kreise, bzw. im Grenzfall in Geraden übergeführt. Dabei sind 
allgemein die Mittelpunktskoordinaten £&,, 7, und Radıen 4 der Ö-Ebene durch die folgenden 
Gleichungen’) mit den entsprechenden Größen .,, 9, und o der z-Ebene verknüpft: 
2: Yy KR’ R-(R’+0’— x, — Y) j 2.0: R? ’ 
TR at PTR) te ’ TR’ tm oe: 
Der Leiter 1 hat in der Ö-Ebene seinen Mittelpunkt auf der „-Achse, die Leiter 2 und 3 
liegen symmetrisch zur 7-Achse und erhalten gleiche Radien. Für die Mittelpunktskoordinaten 
werden gemäß Abh. 2 die Bezeichnungen 4, h, Ik eingeführt. 

Die Anwendung des Abbildungsverfahrens wird dadurch ermöglicht, daß das gesuchte 
Potential p auch als Funktion der neu eingeführten Variablen der Gleichung 


O7, 0°g 


— 0 
0 £ 


r 


ö Ö 7" 
eenügt. Die dureh Gl. (2) einander zugeordneten Punkte der z- und der Ö-Ebene haben gleiches 
Potential. Daher ist die Aufgabe prinzipiell als gelöst zu betrachten, wenn es gelingt, das 
Potential 9 in der Ö-Ebene darzustellen. Um der Bedingung, daß 9 auf der &Achse ver- 
schwinden soll, zu genügen, werden die Leiterkreise an der &Achse gespiegelt, wodurch die 
Kreise 4. 5, 6 neu hinzukommen, denen die Potentialwerte 1, 0, © beizulegen sind. Das 
Potential 9 ist also in der Ö-Ebene so zu bestimmen, daß es folgende Randwerte annimmt: 


auf dem Kreise I +1 | 
m { l 7. 
u 2 359.6 0 | 


ı), R. Liebold: a.a.0d. Gl. 19 und ». 
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Nach einem grundlegenden funktionentheoretischen Satz erfüllen sowohl Realteil als 
Imaginärteil einer analytischen Funktion w = «+ ir der komplexen Veränderlichen © £+ip 
die für Potentialfunktionen eharakteristische Gl. (1). Außerdem sind beide Teile miteinander 
verknüpft durch die Beziehungen: 


ou on ou or 


. i - N). 
0 07 7 08 14) 


Ks gilt jedoch auch die Umkehrung dieses Satzes: Jedes Potential 9 ist Realteil einer analy- 
tischen Funktion, deren Imaginärteil y durch die Gl. (4) bestimmt ist. 


Im vorliegenden Problem wird diese komplexe Funktion in der Form: 


D 6 u 
y Y Y & (h) 
p+iy= \ AU) .In(i—- m)+ \ \ Zn u (9) 
. SZ 7 (Ä Sn) 
int h=1 n=I 


angesetzt, wobei Z, die komplexen Koordinaten der sechs Kreismittelpunkte sind und 4% 
reelle, 7, komplexe Konstanten darstellen. Unter Berücksichtigung der beiden Symmetrie: 
bedingungen bezüglich der Koordinatenachsen: 


plE,n): o(&£,—n) 
el,n)= PI-E,n) 


vereinfacht sich der Ansatz durch Verminderung der Zahl der Konstanten. Wenn die 
Größen 7',® in Real- und Imaginärteil zerlegt werden und allgemein die zu den Kreisen | 
und 4 gehörenden Konstanten durch große Buchstaben, die übrigen dureh kleine ausgedrückt 
werden, so sind wegen der Symmetrie nur einfache, auf den Grad des Gliedes in der Potenz- 
reihenentwieklung bezügliche Indizes notwendig. Der Ansatz (5) erhält dadurch folgende 
(Gestalt: 











C—L)U—L, Ü—{L 
po+tiy: +a,-In.; E)(E HA, ln) ” 
S >2) \> “3 > “j 
— | ( >" Y l ( 1)” 
na \ (a Hab) - €,” (£ I \ (a, Tan En (£ L,y" \ (6) 
B—=} - ö N] 
hs | | ke | | 
r \ an »__fı\2n PP r g \ ib, il __?F ı2n ri n | 
u (Ä GC) (Ä GC) en (Ä Gı) (< G,) 
— nn} 


Dieser Ansatz unterscheidet sich von der Form (32) der Lieboldschen Arbeit nur darın, 
daß in der letzteren die Potenzen der Radien 4 aus den Konstanten herausgetrennt sind. Die 
Konvergenz der Entwicklung (6), die natürlich die Voraussetzung aller Reihendasstellungen 
bildet, ist durch bekannte funktionentheoretische Sätze gewährleistet. 


2. Die Entwicklung des Potentials in der Umgebung der Kreismittelpunkte. Nachdem die 
prinzipielle Lösung für das Potential vorliegt, haben wir uns mit der Ermittlung der Kon- 
stanten a,, A,,.. aus den Randbedingungen zu befassen. Wır wählen dabeı einen Weg, der 
uns dem Endziel, nämlich der Bestimmung der Feldstärke an den Leitern, näher bringt: Wir 
bilden die Entwicklungen für 9 in der Umgebung der Kreismittelpunkte I und 2 und finden 
dadurch eine allgemeine Darstellung der Konstanten in Abhängigkeit von den auf den Kreisen 
gegebenen Randwerten. 


a) Entwieklung für Mittelpunkt I. Ausgehend von Ansatz (6) entwickeln wir 
zunächst die komplexe Funktion 9-+iy in eine Reihe nach Potenzen von (£  £S,). Dabei 
können wir die (© —£,) enthaltenden Glieder einfach beibehalten. Die übrigen Summanden 
mit den Größen (© - £„) lassen sich nach der Umformung: 

= 2 I —6 
EOIERANEEEE: 
1 


sh 





ve 


. . Nun - . . r 
ın Reihen nach Potenzen von (5 A) entwickeln, dessen Absolutbetrag ın der Umgebung 
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“h 
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des Kreises 1 kleiner als I bleibt, so daß die entstehenden Reihen konvergieren. Ist einmal 
die Reihe für In (2 £„) aufgestellt, so folgen alle übrigen durch Differentiieren nach £. 
Es sind: 





FH 
‚1 7 P m 
In (£ Ch) In (£, In) \ ie -) 
SZ m “h -i 
m | 
(4). 
F 
Y Br ' F - \m 
(E-—- &,)"—=(-—1)*- > Wwrn N Bil n—1,2,3 
Bu u Bi Kr er DER a 
m— (0 


Setzen wir diese Reihen in den Ansatz (6) ein und ordnen die einzelnen Glieder nach Potenzen 
von (©  £,), so gewinnen wir die Laurentsche Entwicklung der analytischen Funktion 
p-+iy für die Umgebung des Mittelpunktes 1. Wir wollen jedoch diese Reihe nicht in 
komplexer Form hinschreiben, sondern den im folgenden allein interessierenden Realteil y 
heraustrennen. Wird der Absolutbetrag |ÖÜ —£,| mit r,, und der Winkel, den (© — £,) mit 


der positiven &Achse bildet, mit 9, bezeichnet, so werden: 


In —L)=Iinr, +i9, | 
R Feen —_ pitn.(cosnd, +i-sı },) Br 
(& c,) == Pf, -(cosnvd), + ı-sınnd, 
Ebenso zerlegen wir die in den Entwicklungen enthaltenen komplexen Konstanten: 
Inc, -W)=-ın]|n —Hl+timn+7r) 
(9), 


(un —C)"=|in —L, I" - (cosntr, —i-sinnrt,) 


wobei die von £, nach £,„, gehende Verbindungsgerade mit der positiven &Achse den Winkel z,, 
bildet. Die hierdurch gewonnene Darstellung enthält statt der komplexen Variablen (© -£,) 
die beiden reellen Veränderlichen r, und 9,, wobei natürlich entsprechende Glieder der ab- 
steigenden und der aufsteigenden Potenzreihe jeweils den gemeinsamen Faktor cosn, bzw. 
sinn”, besitzen. Wir wollen daher die Potenzen r,*” und r, ” zusammenfassen und y nach 
Art einer Fourierschen Entwicklung hinschreiben: 


7=(,+(,:0c0529, +0C,-cos4d, +... | 10 

-D,-sind, +D,:sin39, +... J 
Die Symmetrie des elektrischen Feldes zur n-Achse äußert sich darin, daß die Koeffizienten © 
mit ungeradem Index und die Koeffizienten D mit geradem Index identisch gleich Null sind. 


Zur Vereinfachung werden die Strecken 1,,1,,!, eingeführt (Abb. 2) und die Winkel r,, 
soweit sie von einem ganzen Vielfachen von %" verschieden sind, durch die Winkel y y 


ar 
iır/a9/6 
ersetzt: 


| R h—g 
,=y(h-g)'+K y=aretg 
J 
,=V(h+g)"’+%° y,=arctg —_ (11). 
h 
,=2:yh’+%k? y,— arcig 1; 





Die Koeffizienten der Entwicklung (10) lassen sich hierdurch auf folgende Form bringen mit 


p un, 





B A,-Inr, +la,2-In \ +-A,:In2g 
1 
, ‚[,eosny, cosny, Er. ‚(,Sinny, sinny, 
\ (dl, 4 ( 1)’ | + AR Be )+ b,2\ 1)’ H de Er (12) 
n | ==] 
\ Y (- +1 \ "np ( rs 
r Asp (2.9)? ya = Ion-ı 2?” N 
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mal 7 rs A, ry2nla | (4 cos2py, + ki TSIEr 
ur 2p rip °» All 1,2» 2 2» (Ag)?P 
ri \ e 2.2 p+n a Ta und ade cos (2p-+n) ).) 
— "— (2p)!(n —1)! Bere Pi duhee 
n 
(4). N, ep tn! u („Sin @» Fn)y, sin @p-+m)7.) 
u " AP) !(n —1)! Fri BRREN 
n 
| ı \', (2p+2n-1)!( y 2p+2n-2)! (-1P+r"r! 
zen 4 N (An)!l2n I)! agyrr en (2p)H2n 2! Zap +?" ) 
ion nm] n—=] 
In 
| 7 B . 
u; . 2 /,sin(2p-1Uy, , sin(2p -Dy, ı (-1pr 
nit n BE... WERT sın(Zp-I1)y, , sın(2p Yılı Br 
2D—1 rz3D ırr, 02» | + j,?? 1 .u 1,9? 1 | 4 2» | (>q)?P ) (12). 
a £ u 2.2 p+n — 2)! Ta nd Zul. )y sn@p+n re) 
ER — "Qp—V)!m—1)! BreeN Erd dere 
n= 
ı \', 2.(2p—-n u 1)" | cos(?p+n—1y, eoss2p+n EN 
‘ ni9 ' ! ep+n 2» ! 
9) — (2p- D!n-— N)! vr A de 
r, ı 5’, ._@ptn-H pr" 
r) 4 2 I)'!(2n 1)! (AgyPr?" 1 
n==3 
ıb- 
& N. „entinon pres 
a A —N!Rn—2)! (Ag)’Pr?" 
n—=] 
)). b) Entwieklung für Mittelpunkt 2. Indem wir jetzt die Mitte des Kreises 2 zum 
Ausgangspunkt für die Reihendarstellungen machen, finden wir ganz analog die Fouriersche 
Iintwieklung für 9 in der Umgebung des Mittelpunktes 2. In diesem Falle sind jedoch sämt- 
liche Koeffizienten vorhanden: 
d 
h; P=6,+e,:C08S%,+c,:c0os2,-+... \ 
Yo ; , (13). 
+d,-sind, + d,-sın2d, +... | ) 
Hierin ist wiederum die Zerlegung: 
). (© -LS,)=r,: (cos d,—+i- sin d,) 
enthalten. Die Koeffizienten lassen sich durch die Größen (11) in folgender Form darstellen 
mit 9 —=1,2,3,..: 
1 
h.-1, I 
e, a, Inr,+ |a,-In + 4A,:In—* 
k I, 
mr | 7 a 3 | 
\' COSNZ ' a \ sinn = RETIENER 
E ) _ mw. - € 
+ Pi A 1)" (2h)" Top" Ay Tr er e b„( 1)" (2h)" 1" (14). 
4 n==13 ‚a2 
N \ " | cos?ny, m) + \ "n A l)y, ‚Sin(2r -1)y, 
Asp ]..” 2" d en-ı\l j2"-1 he j,."-1 
n—=]1 | ' “1 ; 
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3. Die Bestimmung der Konstanten. Nach Gl. (10) werden die Potentialwerte auf Kreis | 


y Z0,-cosrd, + FI D,-sinrd,, 


wo die €, und D, sich nach Gl. (12) mit r, = 4, dureh die zu bestimmenden Konstanten des 
Ansatzes (6) für g ausdrücken. Nun soll der Potentialwert auf dem Kreis 1 gleich 1 sein. 
Das wird erreicht, wenn C,=1 ist und alle übrigen (,=0 und D,- 0 sind. Setzen wir aul 
den linken Seiten der Gl. (12) diese Werte (, 1, „= D, 0 für vr >0O ein, so erhalten wir 


ein erstes Gleichungssystem für die gesuchten Konstanten a,, b,, A,. By. 


Ein weiteres Gleichungssystem erhalten wir, wenn wir benutzen, daß auf dem Kreise I 

das Potential (s. Gl. (15)) 
2=Lr,:cosrd, +2 d,-sinv d, 

null werden soll. Das ist der Fall, wenn e,=d,--0 gemacht werden. Setzen wir also in 
den Gl. (14), in welehen die e, und d, durch die gesuchten Konstanten ausgedrückt sind, die 
linken Seiten gleich Null, so erhalten wir die noch fehlenden Gleichungen für die a,, by. 
4,. By. Wir haben dann die Lösung unserer Aufgabe auf die Lösung dieser beiden Systeme 
von unendlich vielen linearen Gleichungen mit unendlich vielen Unbekannten zurückgeführt. 
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Für die praktische Rechnung brechen wir die Entwieklung der Potentialfunktion 4 
(l. (6)) bei einer geeigneten Potenz ab. Nehmen wir die Glieder bis (© —{,) ", also bis 
ır n-ten Potenz, so erhalten wir genau so viele Gleichungen wie Unbekannte, wenn wir auf 
ion Kreisen I und 2 den Fourierkoeffizienten bis C,, D,, €, d. die vorgeschriebenen Werte 
Null bis auf €, 1) erteilen. 

Das Verfahren kann auf den Fall beliebig vorgegebener Randwerte ausgedehnt werden, 
n allgemeinen Falle werden die Koeffizienten C,, D,, €, d, nieht null, sondern sie erhalten 
lie Werte, die aus der Fourierentwieklung der vorgeschriebenen Randwerte hervorgehen. 

Die allgemeine Form für die Entwicklung des Potentials y, Gl. (10) und (12) bzw. (13) 
ınd (14). vereinfacht sich in unserem Falle, wenn die Ausdrücke in den eckigen Klammern 
er Gl. (12) und (14) mittels der Beziehungen: 
C,=1, C,=D,=0fürv>0, beir, =4, bzw.c,=d,=0, beir,=4, . . (15) 


eliminiert werden. Wir erhalten: 


für die Umgebung des Kreises |: 











r ww 6 
1 I A,In = HA,| z ‚eos 29, | A,| z sd +. 
| 4, Mr. 4 
(16), 
/ | % a Au j wo 
+ B,| -,)sınd, + B,|-3 —; „)sin 3%, +... 
0 ci: Au Eier 
für die Umgebung des Kreises 2: 
r, | r, 2 me me: 
7 a,In LT ( r, „.)e0s "4 a,» „.)60s2 ,+d, 5 „u) 053 en 
(1m). 
| PB. . | vn, ’ en 
+b, + „.)sin J,+-b, r „sin 29, + er Te) sIn3d, +... 


- - - - 


4. Die Feldstärke als Funktion des Ortes. Es ist bekannt, daß die Feldstärke, die man 
ihrer Definition gemäß darstellen kann durch: 


Y Ö 4 . Ö "\ , ) 
1 - . a aan SE SE ne Ei S), 
F Nr Hi I (15) 


die Maxima ihres Absolutbetrages an den Rändern annimmt. Da es uns lediglich auf die Er- 
mittlung der Maximalwerte der Feldstärke ankommt, so dürfen wir unsere Ableitungen auf 
die Berandungen beschränken. Hierdurch wird die Rechnung erheblich vereinfacht. Auch 
den Bleimantel können wir von unseren Betrachtungen ausschließen; denn indem wir das 


>> 


Innere des Kabels durch reziproke Radien „’ am Bleimantelkreis spiegeln, erreichen wir, daß 


der letztere selbst im Innern des auf diese Weise erweiterten Potentialfeldes liegt. Spiegel- 
bildlich gelegene Punkte besitzen stets entgegengesetzt gleiches Potential, und unendlich kleine 


>> 


\ . . ” > y . . r .. . ı on 
Strecken des Innern erscheinen im Außeren des Kreises im Verhältnis R vergrößert, das 


Potentialgefälle selbst also im umgekehrten Verhältnis verringert. Da sich demnach beim 
Durchschreiten des Bleimantels nach außen der Betrag der Feldstärke vermindert, kann sein 
Maximalwert nur an den Leitern liegen. | 

Iın folgenden soll unter der Feldstärke der negative, ins Innere des vom Isolierstofl er- 
füllten Gebietes gerichtete Differentialquotient des Potentials verstanden werden. 

a) Die Feldstärke in der ©-Ebene. Da die Feldstärke senkrecht auf den Äqui- 
potentiallinien steht, läßt sie sich an den Kreisen 1 und 2 als negativer Differentialquotient 
nach dem Radius darstellen, wobei der Winkel konstant zu halten ist. Wählen wir nun 
für y die Entwicklungen (16) und (17), so gewinnen wir sofort für E,5 die Fourierschen 
keihen auf den Kreisen I und 2, indem wir einfach die Koeffizienten Ü, D bzw. e, d nach r, 


bzw. r, differentiieren, r, =4, und vr, —=4, setzen. Es wird 


für Kreis 1: 


| i S 





EKo3=4A,'; +4, 75;:00os29, +A4,- — ..0084 +... 
R a a, 
(19), 
Me: ; are Bi: 
+ B,:;,:sınd, +B,:; „sin3d, +... 
. ’ı 
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für Kreis 2: 


| 2 4 6 
er i :CcoS».—- ; :Cos> N. - : .ens 2%! 
17: coSV,—+ A, j3 cos24J,-+ a, ja coS3V,-+.. 


.. 


- - y- 


fh (20) 
Ar ae 
+b: za sin ,+b,- ja sin 29, +b,- ; d -sın3d,-+.. 


h. 
2 2 2 





Diese Entwicklungen für E gelten im Gegensatz zu denjenigen für y selbst nur auf den 
Kreisen; denn im Innern des Gebietes tritt im allgemeinen zur Radialkomponente der Feld- 
stärke noch ihre Tangentialkomponente hinzu. 


b) Die Feldstärke in der z-Ebene. Zur Berechnung der Feldstärke im Kabel- 
querschnitt ist es nicht notwendig, das Potential selbst in die z-Ebene zu übertragen; denn 
der Gradient von g.;, läßt sich schreiben: 


dc 


] .. - » 
grad vio|’ dz . . . . . . . . . . . . . (21). 


grad g\>) 


Der hierbei auftretende Umrechnungsfaktor gibt an, in welchem Maßstabe die unendlich 


kleine Strecke |dz| in der Ö-Ebene abgebildet wird. Die Ableitung 2 gibt zugleich Auf: 
Ö ( < 


schluß über die Zuordnung der Differentiationsriehtungen beider Ebenen; denn als Quotient 
zweier komplexer Zahlen besitzt diese Größe ein Argument gleich der Differenz der Argumente 
von d£ö und dz. Also ist dieser Winkel gleich der Drehung, welche die Richtung von dz 
bei der Abbildung in die Ü-Ebene erfährt. 


Wir differentiieren die Abbildungsfunktion (2) nach z und bringen diesen Ausdruck, um 
Absolutbetrag und Argument abzusondern, durch einfache Zwischenrechnung auf die Form: 


dc 2.R i(?-aretg #-;+3) I») 
I Fi be ar he er rt Dar 
Die geometrische Bedeutung dieser Funktion geht sofort aus folgendem hervor: Ziehen wir 
| dc 
dz 
gekehrt proportional dem Quadrat dieser Strecke. Das Argument ist gleich dem zweifachen 


vom Punkt z die Verbindungsgerade nach dem Punkt z2=—+R, so ist der Betrag um- 


von dieser Geraden mit der «-Achse gebildeten Winkel, zuzüglich der Konstanten — . 


Da die zwischen der Feldstärke und den Leiteroberflächen bestehende Orthogonalität 
bei der Abbildung gewahrt bleibt, so entsprechen den radialen Elementen dz ebenfalls radiale 
Elemente d{S. Hieraus ergibt sich eine bequeme Zuordnung der Kreispunkte beider Ebenen. 
Bezeichnen wir mit o den Winkel, den der zum Kreispunkt z gehörende Radius o mit der 
positiven .«-Achse bildet (Abb. 1), so können wir aus Gl. (22) sofort den entsprechenden 
Winkel 9 der Ö-Ebene angeben: 

! IT 
"=06-+2.arctg ai rü d E e  e  Z ne Ze ee 23). 
Hiermit ist die Bestimmung der Feldstärke an den Leitern auf ihre Berechnung in der Ö-Ebene 
zurückgeführt, für welche die Entwicklungen (19) und (20) gelten. 


5. Die Feldstärke als Funktion der Zeit. Haben wir bisher für den speziellen Fall, daß 
das Potential auf dem Leiter 1 den Wert +1, auf den anderen Leitern den Wert Null an- 
nimmt, die Feldstärke dargestellt, so handelt es sich nunmehr darum, die zeitliche Veränder- 
lichkeit der Spannungen in Rechnung zu ziehen. Insbesondere müssen wir die durch die 
zeitliche Abhängigkeit bedingten Maximalwerte aufsuchen. Um die Rechnung zu unmittelbar 
brauchbaren Ergebnissen weiterführen zu können, machen wir die Voraussetzung, daß das 
System der Drehstromspannung symmetrisch, und daß der Kurvenverlauf genau sinus-förmig sei. 

Verallgemeinern wir zunächst das Problem, indem wir den drei Leitern von Null ver- 
schiedene Potentialwerte erteilen, so wird sich im allgemeinen die Feldstärke vektoriell aus 
den von den drei Leitern stammenden Anteilen zusammensetzen. Unter Berücksichtigung 
des zeitlichen Verlaufes der Potentiale ergibt sich hieraus für im Innern des Isolierstoffs ge- 
legene Punkte ein mit der Periode des Wechselstroms rotierendes elliptisches Drehfeld. Die 
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vektorielle Addition geht aber in eine algebraische über, wenn wir uns den Berandungen 
nähern, auf die wir gemäß unseren Feststellungen über den Sitz der Maximalwerte auch die 
‘oleenden Betrachtungen beschränken wollen. 


a)Maximalwertebeinormalem Betrieb. Wir bezeichnen mit E,, E,, E, die Feld- 
stärke an einer beliebigen Stelle der Berandung, wenn die Leiter 1,2,3 nacheinander je das 
Potential +1, die übrigen jedesmal das Potential Null haben. Natürlich sind wegen der 
Symmetrie der Anordnung durch den örtlichen Verlauf der Funktion E,, deren Ermittlung 
unsere bisherigen Betrachtungen galten, auch die Größen E, und E, bestimmt. Bedeutet 
ferner V den Scheitelwert der Phasenspannung (Leiter gegen Erde), & die Kreisfrequenz, so 
ist zur Zeit t die Feldstärke: 


/ 


a) + E,: c08 (® ! 3 a) 


/ 


‘) 


E V: B -c0os @t + E,: «os (ot 5 





Dieser Ausdruck läßt sich auch auf die Form bringen: 


E=V.- 





E(0) -coswt—+ E35) - sin ot ’ 


/ 


4 ® y hj I ® 0 u. JT Y .. ® 
wenn wir mit E(0) und E(5) die für ot=0 und &t-— vorhandene Feldstärke bei der 


/ _ 


Scheitelspannung V 1 bezeichnen. Es sind: 


| 
E(0) E, >) (E,+E,) 





Setzen wir den Differentialquotienten von E nach t gleich Null, so erhalten wir für den Zeit- 
punkt des Maximalwertes: 


(2 


tg wt= E (0) R 


und für das Maximum der Feldstärke selbst: 


MAX 


4 2) Y 7T R r v9 er) 9 Y ’ Y y ı h) Fr Tand 
Ber | E(0”’+KE | >) =V-yE?+E?+E?—E.: EK,—E, E,—E,:E, . (28). 
b) Maximalwerte bei Erdschluß. Neben dem normalen Betriebszustand interessiert 
noch besonders der Fall des einphasigen Erdschlusses, da hierbei wesentlich größere Feld- 
stärken auftreten können. Nehmen wir an, daß etwa die Phase 3 ohne Widerstand mit dem 
Erdpotential verbunden sei, so rückt der Nullpunkt des Spannungsdiagramms (Abb. 3) nach 
dem Punkt 3. Die Potentiale der Phasen 1 und 2 steigen auf den Y>-fachen Wert des nor- 
. . r > ER . ’ 
malen Betriebs; denn die Vektoren 31 und 32 stellen jetzt die Phasenspannungen der unge- 
erdeten Leiter dar. Sie schließen einen Winkel von 60° ein. Da der Leiter 3 keinen Bei- 
trag liefert, erhalten wir für die Feldstärke bei Erdschluß: 


JT 


‘) ’ 
pa) 





F.e=V:.y3- | E,(0):-cosot + E, | -sinof 


wenn wiederum 








| 
Ee (0) h + .) -E, 





Br R a Bc“ 
E.(5)=" EB, 


> 


eesetzt wird. Schließlich ist der Maximalwert bei Erdschluß: 


E, vy3\/ Bor + RZ) =VVS VER+RSHEBE,  . ;° 


e max 
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6. Die Kapazität. Stellen wir in der C-Ebene nach dem Gaußschen Satz die Ladungen 


der Leiter I und 2 als Verschiebungsfluß durch die Oberfläche dar, und benutzen hierbei die 


Entwieklungen (19) und (20) für die Feldstärke, so zeigt sich sofort, daß bei der Integration 
sämtliche Glieder mit Ausnahme des vom Winkel 9 unabhängigen Gliedes verschwinden. 
Für die Potentialverteilune (+1,0,0) sind daher die Ladungen, wenn & die relative Dielek- 
trizitätskonstante des Isolierstoffes ist, und die Potentialeinheit 1 Volt beträgt: 


€ £ 


. ! » ö (' » 8 
18.106 % Clkm, (gg. FR er 


I) ) 
In der z-Ebene sind die Ladungen ebenso groß, da von entsprechenden Leitern beider Ebenen 
oleichviele Verschiebungslinien ausgehen. 

Zwischen den Ladungen 4 und den Potentialen p eines Mehrleitersystems besteht be- 


kanntlich lineare Abhängigkeit, deren Koeffizienten « nur von der räumlichen Anordnung 


und vom Isolierstoff abhängig sind. Bei unserem symmetrischen Dreileitersystem reduziert 
sich die Zahl der Koeffizienten auf zwei: 

Pı = YIı' A r(@%+ (3) Gy: 
l’> (A, tr: +4):0, 
ls G,'Aı Tr (q, 42)" As 


Nach g aufgelöst lauten diese Gleichungen, wenn statt der Potentiale 
deren Differenzen eingeführt werden: 
"Pi P>» Ps) i ki: | 
2, —Pı Ps)’ ko to Zn ir Hl ARE Sr 
"Ps Pi pP.) Ki: 


IV 


Br k;, ut 
‘fs Pa‘ k;: I ( 
VE Ps’ k: rt 


wi 





De 


Die Größen A,,, #,, werden allgemein als „Teilkapazitäten* bezeichnet, da man sie als Kon- 

densatoren auffassen kann, die zwischen Leiter und Erde bzw. zwischen die Leiter ein- 

geschaltet sind (Abb. 4). Sie sind mit den „Potentialkoeffizienten“ «a,,, a,, verknüpft dureh 
| A. 


k,= ne, ek, —® 5 34 


Au ta, - Ai An 
Endlich interessiert noch die für symmetrische Spannung definierte „Betriebskapazität“ 


CGh,=k.,+8:E, ea % 


Da wir die Ladungen 4,. 9, für p, +1. pp, = 0d unter Gl. (28) berechnet haben, lassen 
sich die Kapazitäten sofort aus den Gl. (29) angeben: 


N € 
» ’ E, } Il: ’ ö ‘'ıl» 
FR m (A,+2-a,) u«F/km B.; Tue Fk m 
(32). 
(,, ” (A, 4,) «Fk m 





II. Numerischer Teil’). 


1. Die Konstanten der Potentialfunktion. Wie aus unseren allgemeinen Ableitungen her- 
vorgeht, ist es nieht möglich, die Feldstärke und die Kapazität als einfache Funktionen der 
Kabeldimensionen darzustellen. Bei der Untersuchung dieser Abhängigkeiten sind wir daher 
eezwungen, die Berechnungen für diskrete Zahlenwerte durchzuführen. Zunächst ist es klar, 


) 
und —, bestimmt 


daß der Verlauf des elektrischen Feldes nur von den Größenverhältnissen R f 
[1 l 


wird. Wir wollen daher R konstant halten und 
R lem 
setzen. Im folgenden sind also o und r als unabhängige Veränderliehe zu betrachten. Zur 
Vereinfachung werden wir noch das Potential bzw. den zeitlichen Maximalwert der Spannung 
5), Die numerischen Reehnungen erstreekten sieh ursprünglieh der Aufgabe der Dissertation entsprechend nur 
auf die Ermittlung des absoluten Optimalfalles. Da ein allgemeineres Interesse an den Ergebnissen angenon men werden 


konnte, wurde über den Umfang der Dissertation hinaus im Rahmen des „Ingenieurdienstes" der gesamte, praktisch 
interessierende Dimensionsbereieh untersucht und die gewonnenen Zahlenwerte mit den früheren zusammengestellt. 
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‚leieh 1 setzen und als dimensionslose Zahl ansehen, wodurch die Feldstärke die Dimension 
m ' erhält. Da die Feldstärke dem Scheitelwert V direkt, dem Bleimantelradius Z# beı sonst 
»leichen Größenverhältnissen umgekehrt proportional ist, so erfolgt die Umrechnung der von 


r 


Be ra 
ıns eewonnenen Zahlen einfach durch Multiplikation mit R° 
l 


Zur Ermittlung des Optimalfalles wurden die numerischen Einzelfälle der Tab. 1 zalılen- 
näßig durchgerechnet. Es ergeben sich für die Entwieklungskonstanten des Potentials 9 die 
'n Tab. 2 wiedergegebenen Werte. Man erhält eine Kontrolle, wenn man die Potentialwerte 
der Näherungslösung auf den Kreisen 1 und 2 ausrechnet und mit den vorgeschriebenen 
Werten 1 und ®O vergleicht. Als Maximalfehler ergaben sieh z. B. für das Beispiel 11 
o—=V0>2em, r=05cem) auf Kreis 1: —0,000.,39, auf Kreis 2: — 0,001.59. 


Tabelle 1: Radien und Koordinaten der £-Ebene. 








R=lem Kreis |] Kreis >? 
Na 0 r da 7 Mi Jı k 
a; em em em em em em em 
1 0.05 0,40 0.051.086 0.430,40 0.132.015 1.112.22 0.014.62 
> 0.05 0.50 0.044.494 0.334.82 0.133.738 1.006.7 1.158.57 
2 0.05 0,56 0.041.134 0.283.37 0.133.138 0.017.19 1.291.309 
4 0.10 0,40 0.102.564 0,435.90 0),266.667 1.133.33 0.923.76 
5 0.10 0.50 0.089.286 0.339.238 0.270.27 1.027.02 1.170.3 
6 0,10 0.56 0.082.52? 0.287.34 0.268.96 0.936.533 1,304.41 
7 0.15 0,42 0.150.455 0.424.34 0.408,77 1.152.87 0.091.22 
Ö 0,15 0.50 0,134.678 0.346.80 0,412.37 1.061.88 1.190.414? 
') 0.15 0,56 0,124.426 0.294.01 0.410.34 0.969.62 1.326.7 
10 0.20 (0.44 0.196.695 0.416.21 0.560.54 1.186.14 1.067.098 
11 0.20 0.50 0.180.996 0.357.47 0.563.338 1.112.68 1.219.75 
2 0.20 0.56 0.167.11 0.303.48 0.560.54 1.017.095 1.359.25 
13 0.25 0.46 0.241.659 0.411.24 0.725.558 1.234.79 1.156. 
14 0.25 0,50 0.228.57 0.371.42 0.727.27 1.181.78 1.259.65 
15 0.25 0.54 0.216.53 0.333.385 0.725.58 1.118.71 1.357.28 


3, Maximalwerte der Feldstärke bei normalem Betrieb. DBezeichnen wir wie früher mit 
E,,E,, E, die Feldstärken, welche auf der Berandung hervorgerufen werden, wenn die 


leiter 1, 2, 3 nacheinander je das Potential +1, die übrigen das Potential Null haben, so 


erhalten wir nach Gl. (25) als zeitlichen Maximalwert der Feldstärke für die Scheitel- 
spannung V 1 


F-yE’tE’+E'-E.5-E.5-E.E,. 
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Abb. 5. Feldstärken Eı, Es, E3 am Leiter 1. Abb. 6. Maxim. Feldst. an den Leitern bei norm. Betrieb. 
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Wenn die Potentialfunktion 9 berechnet ist, welche auf Leiter 1 den Wert +1, auf Leiter 2 
nd 3 und dem Bleimantel den Wert Null annimmt, so erhalten wir E, nach den Formeln (19, 
»). 21. 22). Damit sind auch F, und E, bekannt; denn wegen der Symmetrie ist die Feld- 
tärke E, auf dem Leiter 1, wenn Leiter 2 das Potential +1 hat, und di® übrigen das 
Potential Null, gleich der Feldstärke auf Leiter 2, wenn Leiter 1 das Potential +1 hat, und 


‚lie übrigen das Potential Null. 


Abb. 5 zeiet für das numerische Beispiel I1 (R=1em, 0=0,2cem, r=-0,5 em) die 
Werte von E,.E,, E, für den Leiter 1. Die daraus für jeden Punkt des Umfanges berechnete 
maximale Feldstärke ist in Abb. 6 dargestellt. 


Wesentlieh ist hier — und das gilt auch für alle anderen Fälle des untersuchten Be- 
reiches —, daß die örtlichen Höchstwerte F der zeitlichen Feldstärkenmaxima bei o, = 0° und 
s,— 180" auftreten, d. h. an den der Kabelmitte und dem Bleimantel zunächstliegenden 
Punkten. Tab. 3 gibt die Werte von F an den Stellen 6, 0", 30°, 150° für die gerechneten 
Fälle. 


Tabelle 3: Maximalwerte der Feldstärke bei normalem 
Betrieb (Y =]. R lem). 








0 r '® ‘3 jr‘ıs0) 

Nr. ei 
em em em em em ! 
| 0.05 0.40 300 8.80 8.01 
2 0.05 0,50 8.34 S.50 8.29 
> 0.05 0.56 s,14 s.15 8.52 
| 0.10 0,40 6.88 6.72 5.46 
5 0.10 0.50 5.97 5.00 3.84 
6 0.10 0,56 5.66 5.63 6.24 
| 0,15 0,42 6.40 6.15 1.70 
8 0.15 0,50 5.42 5.27 5.17 
Y 0.15 0.56 4.96 4.80 5.73 
10 0,20 0.44 6.48 6.25 1.69 
11 0.20 0.50 2.49 5.28 5.18 
12 0.20 0,56 4.86 1.69 2.98 
13 0.25 0,46 6.87 6,88 5.06 
14 (25 0.50 5.04 2.80 2.56 
15 0.25 0.54 5.27 5.22 60 


Zur Ermittlung der optimalen Kabelabmessungen gehen wir nun so vor: Wir halten 
zunächst den Leiterradius o fest (immer bei R=1cem) und verändern r. Für kleines r liegt 
dann der ungünstigste Punkt (größtes F) bei o,=0" zunächst der Kabelmitte. Vergrößern 
wir r, so nimmt an diesem Punkte die größte Feldstärke F“” ab, an dem dem Bleimantel 
zugekehrten Punkte o, = 180° wächst sie. Der Optimalfall ist erreicht, wenn die Maximal- 
werte an beiden Punkten gleich sind. Er wurde dureh zeichnerische Interpolation ermittelt. 


Tabelle 4: Maximale Feldstärke bei optimalen Verhältnissen und bei 








Bemessung nach V.D, E.-Vorsehrift (I 1. R lem). 

Optim. Bemessung Bemessung nach V.D. E.-Vorsehrift 
u en) (0) (IS) 

0 y F* r F% F Atk l h 
em em em ! cm em ! em ! em ! 
0,05 0.513 8.28 0.384.3 918 3.80 12.69 
0.10 0.512 5.90 0,402.6 6.85 7.18 8.098 
0.15 0.514 5.29 0,420,9 6,39 6.46 Ss.u0 
0.20 0,513 2,32 0,439.2 6.50 6.36 7.82 


0.25 0,511 3.7 0.457.5 6,06 6.62 8,38 
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Abb. % zeigt die Interpolation für die Fälle o = 0,25 und 0,05 em, Tab. 4 gibt für o= 0,05: 
0,10, 0,15, 0,20, 0,25 em die optimalen r-Werte mit den zugehörigen Werten der größten Feld- 
stärken pı» ps F* 


Tragen wir jetzt F* als Funktion von o auf (Abb. 8) und verbinden die berechneten 
Punkte durch eine Kurve, so bekommen wir als absolutes Optimum o=0,17 cm. Dazu ge- 
hört der Wert r= 0,513 em (wo die letzte Stelle unsicher ist). Der zugehörige Maximalwert 
der Feldstärke ist F* 5.26 cm ' (letzte Stelle unsicher). 


























is 0 
| ,7# | | | 
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em” f f9 4 Gascm | | 
a Fi | | 
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><P- 42500 | | ? 

(180) | | 7 
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| | 4; 

| | 

| | 1 

| | | 
| | | 
VB l l i u l a l |; L | | 
03 2W 0,45 00 0355 cm 0 0,05 00 0% 020 cm 025 

Asa?) r (Ar54.8] p 
\bb. 7. Interpolation der optimalen Abmessungen Abb. 8. Optimale Feldstärkehöchstwerte F* in 
und der zugehörigen maximalen Feldstärke Abhängigkeit vom Leiterradius % (R=|1cem). 
FW — FU8V) —_ FF(R=1cm). 


Der V.D.E. schreibt nun für Drehstromkabel gleiche Isolationsstärken zwischen den 
Leitern einerseits und zwischen Leitern und Bleimantel andererseits vor. Diese Forderung, 
die durch die Beziehung 


erfüllt wird, führt bei R = lem zu Mittenradien r<0,5cem. so daß der absolute Maximalwert 


von F in dem untersuchten Bereich stets bei o, - 0° liegt. Wir berechnen zu den Leiter- 


Tabelle 5: Maximalwerte der Feldstärke bei einpoligem Erdschluß 
(V=1, R=1cm). 








Unzeerdeter Leiter 1 Geerd.l..3 
\ 0 Y f‘® F' 30) F' 180) F630 F240 
NT. 
em em em! em em’! em ! em ! 
| 0.05 0.40 11.50 11.70 12.78 11.45 5.03 
> 0.05 0.50 11.53 11.73 13.55 11.56 3.52 
B 0.05 0.56 11.78 11.94 14,24 11,81 2,82 
| 0,10 0.40 a 5.02 8.96 1.68 5.32 
5 0,10 0.50 7.34 1.62 0.83 7.36 3.74 
5 0.10 0,56 1.32 7.62 10.62 7.40 3.04 
fi 0.15 0.4? 6.69 6,099 7.98 6.683 5.86 
N 0.15 0.50 6.01 6.37 S.SS 6.04 4.38 
ig) 0.15 0.56 5.79 6.19 9.90 5.090 3.56 
10 20 0.44 6.51 6.89 7.84 6.49 6.41 
11 0.20 0.50 5.71 6.04 8.8? 5.66 5.06 
12 0,20 0,56 5.28 5.59 10.27 5.24 4.09 
13 0.2 0.46 6.79 1.43 s.44 7.00 7.12 
14 0.25 0.50 5.92 6.44 9,42 6,06 5.96 


10.73 2.46 5.09 
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adien 0 = 0,05, 0.10, 0,15, 0,20, 0,25 em nach Gl. (33) die vorgeschriebenen Radien r und 
uchen zu diesen durch graphische Interpolation (s. Abb. 7) die zugehörigen Höchstwerte F'". 
ie Ergebnisse sind in Tab. 4 zusammengefaßt. 


An vorhandenen Unterlagen, die zum Vergleich PR — ur 
mit den von uns berechneten Zahlen herangezogen werden 8 | Va N = 
können, kommt vor allem die von R. W. Atkinson‘) y; 


auf Grund von Messungen aufgestellte empirische Formel 
in Frage. Dieser Autor findet für Kabel, die etwa den 
V.D.E.-Vorschriften entsprechen, den Höchstwert der 
Feldstärke ebenfalls an der Stelle o,=0": „Dieser Wert 
ist in allen praktisch wichtigen Fällen gleich dem Höchst 
wert in einem Einleiterkabel vom gleichen Leiterdurch- 
messer, dessen Betriebsspannung gleich der Phasen- 
spannung des Dreileiterkabels ist, und dessen Isolierungs- 











dieke gleich ist dem Abstand zwischen Leiteroberfläche | - ur 
und Kabelmittelpunkt im Querschnittsbilde.“ Verfahren oo 6 70 00 20 300 360 
wir nach dieser Rechenvorsechrift, die sich in der Form 5 £ 
cm” 
* ! Leiter 3 
V 7 
u) . 8 
FyiK: - 3: | | 7 
ÖO* In F | 1 N N N 
C 0 2 
1} 
. » . . ° 7 (0 
sehreiben läßt, so erhalten wir die in Tab. 4 als F', u | | || 
schrei ul : 2 F r ei 0 90 7 70,100 20  2W 
eingetragenen maximalen Feldstärken. Der Vergleich Asa * 
mit unseren Zahlen F“® zeiet, daß die Formel Atkin- \ 
: 2 Be : Abb. 9a und b. Maximale Feldstärke a. d. 
sons eine gute Näherung darstellt, da der Leiterradius Et 0 Su en 
o— 0,05 cm bereits außerhalb des Bereichs der praktischen Scheitelw. d. Phasensp.  —1. 


“ a Dim.: R=lem, P=(N?em, r=(W,dem. 
Anwendungen liegt. ' 


3. Maximalwerte der Feldstärke bei Erdschluß. Wir berechnen jetzt ganz analog nach 
Gl. (27) auch für den einpoligen Erdschluß die zeitlichen Maximalwerte F, der Feldstärke 
bei der Scheitelspannung V=1. Den Verlauf dieser Größe am Umfang der beiden ungeerdeten 
sowie des geerdeten Leiters ermitteln wir wieder für das Beispiel I1 und erhalten das in 
Abb. 9a und b dargestellte Ergebnis. Da sich gegenüber dem normalen Betrieb die Spannung 
leiter gegen Mantel auf den y 3-fachen Wert erhöht hat, während die Spannung zwischen 
den Leitern die gleiche geblieben ist, so zeigt F, bei o, = 180° ein starkes absolutes Maximum. 
Dieses ist auch in allen übrigen Fällen vorhanden. Der Höchstwert auf der Innenseite des 
Leiters, der nur bei extrem kleinen Radien r als solcher in Erscheinung tritt, liegt an der 
dem geerdeten Leiter 3 benachbarten Stelle o, = 30°, 





Tab. 5 enthält die Werte F, des ungeerdeten Leiters 1 bei o, =0", 30°, 150°, 330°, so- 
wie in der letzten Spalte den stets bei o, 240° gelegenen Maximalwert des geerdeten 
Leiters 3. Wie aus Tab. 5 hervorgeht, liegt bei Bemessung des Kabels nach der V.D.E.- 
Vorschrift der absolute Maximalwert von F, stets bei o, = 180°. Für die Leiterradien 0 0,05, 
0,10, 0,15, 0,20, 0,25 em interpolieren wir auch hier wieder graphisch die zu den vorge- 
schriebenen Radien r gehörenden Werte F,“*® (s. Tab. 4). 


4. Die Kapazität. Da uns die Konstanten A, und a, bekannt sind, können wir aus den 
Gl. (32) die Teilkapazitäten und die Betriebskapazität ermitteln. Wir setzen e= 1 und er- 
halten die in Tab. 6 unter k,,, k,s, ©, eingetragenen Zahlen. Wir wollen diese exakten 
Werte noch mit den sich aus der bekannten Näherungsformel ergebenden vergleichen: 
L. Liehtenstein’) spiegelt die Leiter in der z-Ebene am Bleimantel und erhält für die 
Potentialkoeffizienten: 


1S R’ — r? km ) R’+-r’+r?:- 7 km 
a, —o. n . a,,l) -]n - - = 4 
& Ro uf ö 3 31° uf 


6%, R.W. Atkinson: a.a.0. 


”), L,. Liehtenstein: a.a.0. 
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wobei er der Bedingung Rechnung trägt, daß das Potential am Bleimantel verschwindet. 
Wir berechnen nach diesem Verfahren unter Benutzung der Gl. (30) die Kapazitäten und 
tragen ihre prozentualen Abweichungen von den exakten Werten in den letzten drei Spalten 
der Tab. 6 ein. Wir bemerken allgemein ein starkes Ansteigen der Fehler in den Teil- 
kapazitäten bei wachsendem Leiterradius. Bei der Berechnung der Betriebskapazität heben 
sich die Fehler nahezu auf, so daß eine sehr gute Annäherung resultiert. 








Tabelle 6: Die Kapazitäten bei ee 1llR= lem). 
(L (L «L, 
\r VO / hrı Kia Gy Akıy. A hs ACH 
em em uFikm uF/km uF/km n 07 97 
7 7 7 
| 0.05 0,40 0.014.075 0.002.831 0.023.469 0.4 0.5 — 01 
> 0.05 0.50 0.017.058 0.001.942 0.022.883 0.3 0,9 0,0 
> 0.05 0.56 0.018.420 0.001.525 0.022.094 0,3 - 1.2 0.0 
| 0.10 0.40 (018.700 0.004.876 0.033.529 > 1.4 0.6 
5 0.10 0,50 0.021.045 0.003.370 0.032.055 — 1.7 1-33 0.4 
ni 0.10 0.56 0.024.213 0.002.685 0.032.269 1.6 4.5 0.1 
7 0.15 0.42 0.023.085 0.006.843 0.043.614 1.8 >26 1.3 
OÖ 0.15 0.50 0.026.774 0.005.057 0.041.944 1.3 + 6,6 0.4 
) 0.15 0,56 0.030.140 0.004.050 0.042.288 — AB 1-97 0.3 
10 0.20 0.44 0.028.452 V.OOO,.LTS 0.055.065 8.0 40 I 
11 0.20 0.50 0.032.306 0.007.191 0.053.878 8.6 —- 10.5 1.0 
1? 20 0,56 0.037.211 0.005.724 0.054.384 8. — 16.7 0,» 
15 0.25 0,46 0.035.452 0.012.138 0.071.867 14.4 5.7 1.2 
14 0.25 0.50 0.039.205 0.010.003 0.069.485 14.9 - 15.4 2.4 
15 0.25 0.54 0.043.799 0.008.537 0.069.408 15.2 21.1 1.8 


C. Schluß. 


Die der Lösung unserer Aufgabe zugrunde gelegten Idealisierungen: Homogenität des 
Isolierstoffes und glatte zylindrische Leiteroberflächen, bedingen eine Abweichung der Er- 
vebnisse von der Wirklichkeit. Die erste der beiden Bedingungen kann nun für die in Frage 
kommenden Hochspannungskabel als praktisch erfüllt angesehen werden, da bei diesen für 
die Leiter- und Gürtelisolation und auch auf die Zwickel wohl durchweg das gleiche Material 
verwendet wird. 


Schwerwiegender erscheint die Erhöhung der Feldstärke, die durch die Unterteilung der 
leiter in Einzeldrähte hervorgerufen wird. Dieser Einfluß läßt sich in erster Annäherung vom 
Hauptproblem absondern und getrennt behandeln, wobei zweckmäßig die von W. Deutsch’) 
für das Einleiterkabel mit unterteiltem Leiter entwickelte Lösung benutzt wird. Wir könnten 
etwa so verfahren, daß wir für ein dem Dreileiterkabel entsprechendes konzentrisches Kabel 
die Feldstärkeerhöhung berechnen und dadurch einen prozentualen Zuschlag zu unseren 
Werten gewinnen. 151 


8) W. Deutsch: ETZ 1911, S, 1175. 
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KLEINE MITTEILUNGEN 


Fourier-Analyse von Kurvenzügen, die 
stückweise aus Ausschnitten von Sinus- 
linien zusammengesetzt sind. Abb. 1 bringt 
ein Beispiel für eine Art von Kurvenzügen, die in 
der Elektrotechnik neuerdings häufig anzutreffen 
sind. Es bezeichne (vgl. Abb. 1) 


m Yn sin A (ot DB re 


den Verlauf der den n-ten Abschnitt bildenden 
Sinuslinie und / das Verhältnis der zugehöriren 
Kreisfrequenz r=4@ zur Kreisfrequenz © der 





1 y 
- FM, . - 
but 9 Sm. m SM NRom-ı.m ie 
,. n (5b). 
we=?2 | 
._ > * » J < \ 
NSm-ı.m'COSNSm-ı1.my 


d.h. die in Gl. (2) geforderte abschnittweise Inte 
eration reduziert sich auf eine einfache Summen 
bildung, in die alle Sprungstellen zwischen 0 und 
2 ı mit den zugeordneten Sprüngen s des Funktions- 
wertes und s’ von dessen Ableitungen eingehen?!) 








Abb. 1. Kurvenzug, der stückweise aus Ausschnitten von Sinuslinien zusammengesetzt ist. 


(rundwelle des Kurvenzugs. Für die Fourier- 
Koeffizienten folgt dann mit ®f=r und Gl. (1) 
In 
Aa,„=\yeosnrdı 
0 
: (2a) 
p Sm mt 
22 Yo ı sinA(z "y)-ecosnrdr 
Im = = ; 
Sm 1m 
und 
IT 
b,c=\ysinned.r 
u 


(2b), 
Sms md 


u 
= 2 Yın \ sinA(r —y,.)sinnzdı 
mal }_ 





Sm—-ı m 
wobei die & gemäß Abb. 1 die Übergangsstellen 
bedeuten und &,,=0 sowie S,..r+41 =27T zu Setzen 
sind. 


a) Für »==% wird 


\sinZ(.r —y)-cosnzdı= 


n? — A2 . 


(OA). 
| 


’ ch». ’ 
COS Al2 — Y)-Ccos n.ı sınz(l.£2-Yy)sinnz? 


\ ı J 


\sink(e - v)sinnzedr= 





P- A“ a 
r " ' (3b). 
cos (2 -y)-sinnz . sin (ar —W)-Cos nz, 
Setzt man noch 
Yın sin A (Sm-ısm m) \ (4a) 
Yn-ı sinAlm-ı.m Yn-ı) —Sm-ım 
} m COS AlSm-ı.m Ym) 
1 (4b) 
y > Mr a u : k ’ * “ 
Yın ı COS AUS mn - um Y m-ı) 7 S m-ı.m 


wobei Sy, _1.m den Sprung des Funktionswertes 
y und $S’y_..m den Sprung der Ableitung y’ an 
der Stelle &,,_4.m bedeuten, so erhält man nach 
Einsetzen von Gl]. (3) in (2) die einfachen Aus- 
drücke: 


r 
1 \' 
2 vor Ne e 
TEE Pr \S nm COS NSm-ı.m . 
u“ I va). 
m), 


( 
> ‘7 < x 
PNSm-m’SIMNSm-ımy 


b) Für „=% ergibt sieh 


\sin 4 (x p)*COS Azdı | 
1 ; . — ee Bo (ba). 
= .; 1 — Ar-sinAy—+sinAk (2 — y) sin r) | 
2A ; 
\sinA(r y)sinArd.ı | 
1 7 j 1 ee (65). 
= — ıAT-Cc0SAwy cos A (ır yv)-sin Au | 
I) / / 5 
Setzt man zur Vereinfachung 
S mm +1 Sm-.m Im | 
Ynı sin yn = Ymo (7) 


a er 
} m COS A ym = j J mo 


so folgt aus Gl. (2) im Zusammenhang mit (4) 





r 
y 
ad, + TI J;. \ PP +9 
A’: 3 . En"Y/ mo 
Di, 
m==| 
(Sa). 
. 
WA 
Syn 1. m'>IN A Sn 1 0 
m y4 
N 
y 
II \ 
! : 
2 I > = Jd [77 HY nu 
2,iN\ 
m | S| 
( ’ 
r 
Bi % \ 
r . 
u zul Sm- em" ISIN Acın-ı.m 
u J 
mMm== 2 





Für den Koeffizienten n=0, der bei Wechsel- 
stromkurven die Gleichstromkomponente darstellt. 
ergibt sich aus der Beziehung 


On 
aön—iy dx TE EB EB (Y) 
m 
ähnlich der Herleitung von Gl. (5a) der Wert 
. 
| y 
d,° 2ı= ru \ Sm 1: Mm at 
a y 
m=—=2 
(Aus dem Forschungs-Institut der AEG.) 
Berlin. R. Feinberg. 423 


ı) Vel. G. Koehler und A. Walther: Arch. Elektro- 
teehn. Bd. 25 (1931), S. 147. 
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Spannungsverteilung in der Guriplatie 
einer Rippendecke. |)ie Ungleichmäßigkeit 
der Spannungsverteilung in Gurtplatten von Trä- 
eern ist von Th. v. Kärmän und G. Schnadel 
mit Hilfe der Airyschen Spannungsfunktion und 
des Prinzips der kleinsten Formänderungsarbeit be- 
handelt worden. 

Die Integration wurde von den oben Genannten 
und ihren Schülern mit Erfolg durchgeführt in dem 
Falle der unendlich langen Träger mit periodischen 
Momentenflächen bzw. in dem Falle endlicher langer 
Träger mit starren Endquerschnitten. 

Es läßt sich nun ein weiterer, verhältnismäßig 
leicht lösbarer Fall angeben, wenn man die Airy- 
sche Differentialzleichung für Volumen- und Ober- 
tlächenkräfte erweitert und Periodizität nicht in 
der Riehtung längs der Spannweite, sondern quer 
zur Spannweite voraussetzt. 

Außbere Kräfte in den Gurtplatten müssen an- 
resetzt werden, um die Wirkung der Stege (oder 
Rippen) wiederzugeben, und Periodizität hat man 
bei vielen Rippen und gleichen Rippenabständen 
und gewisser ziemlich allgemeiner Symmetrie der 
Belastung. 

Man erreicht dann eine vollständige Erfüllung 
der Randbedingungen an den Trägerenden. 


ds _ 


y 2 








Pas - 








6er] 6 dF- an BE 


























P 
Z 
„A y EAg 
ua az 114 t 
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D -_— b -— 7) 


Auf die Gurtplatte werden von 
ausgeübt, die als 


Gurtplatte: 
den Stegen Schubspannungen 


Kräfte bezogen auf die Volumeneinheit, einge- 


( 


führt werden können. 


Die Gleichgewichtsbedingungen lauten dann: 


( ( 
A N r,.Y) 
O. O4 
(1). 
N a 
() ( OT 
. 0 | 
nr; Or 
Setzt man nun: 
0: «D Oo 
do Öö ) u, 67 o 9 
\ / Op? 2 dr? | 
(2) 
N2«cd 
( : 
T . Op \ N dr| 
0r04Y 


so wird die Airy sche Differentialgleichung zufolge 
der Komptabilitätsbedingungen («u Poissonsche 
Zahl): 


| | «b fi ! Wi) 
+ #7 +) 
O’0p O0» (N) 
} Bu [Ki = - | 
O7“ 04° 


Der Natur der Aufgabe nach ist hier f als eine 
in 4 periodische Funktion mit der Periode 5 und 
als in. .r beliebig vorschreibbare Funktion anzusehen, 
von der nur verlangt werden soll. daß sie überall 
endlich und nur an der Stelle der Rippen, d.h. 
über die Breite d von Null, verschieden ist. 


Math. und Mech. 
Man kann z. B. entsprechend Abbildung setzen: 
2 ’ In um d 
A X,(a 11 cos | y)) von y „bis *® 
i 2 f d d 
A=0 füry< ‚und > 
2 
dann wird: 
an ; 
Op: | COS ! | X,dz 
) | | 19 } 
und damit nach G]. (3): 
un 2n\? f 
— — uNX, cos yv\luXı 03 
/ Ä 7’ a). 
Be, 1, 
von Y bis -- > 
- d .d 
Jj=9 für y =und>+7- 
Zur Integration setze man nun in bekannter 
Weise: 
2 2enny — ER Inn 
D = 2 3.608 De, Dr (4). 


dann erhält man durch Einsetzen in (3) für die Z,, 
die Differentialeleichungen: 


h 
b f 2 
— | J-cos(vy)dy . ©). 


„ Er 
y? pn 


1 _ 1] — N 


) y? ee 


bekannte 
etwa zu be- 


rechte Seite als nun 
werden. 


Hier kann die 
Funktion von .r angesehen 
zeiehnen mit: 


h 
gs 4 
. Joa y)ay—md. : .: 2: iR 


Die Integration von (5) ergibt dann: 


— Zn,p-+ A„ Sin (ve) 


I 


B„&oj (vr) 


a (6). 
D,rzM%oj (vr) | 


C„rSin(vr) 


so daß man, um die Randbedingungen der Gurtplatte 


_ on RD 
an den Träzeren len zu befriedigen [z. B. 0x - 
q | . Ri) Yy? 
N” oD / 
( 
0 =U % — I) für r 00 ‚die rie 
p aröy ir +5 lie rich 


tige Zahl von Konstanten zur Verfügung hat. 4,. 
B,. C„ und D„ sind also vollständig durch Zu.p 
bestimmt. 

Gurtplatte und Steg. Damit ist für ge- 
gebenen Schubspannungsangriff der Stege an der 
Gurtplatte die Aufgabe gelöst, dagegen ist eine 
weitere Betrachtung nötige, wenn dieser Schub- 
spannungsangriff erst aus dem gegenseitigen Zu- 
sammenhang von Gurtplatte und Stegen gefolgert 
werden muß. 


In der Literatur ist bisher dieser Zusammenhang 
aus dem Minimum der Formänderungsarbeit be- 
rechnet worden. Er läßt sich jedoch nach einem 
Vorsehlage von E. Reißner. der an anderer Stelle 
ausführlicher von ihm diskutiert werden wird, ein- 
facher ermitteln. Berechnet man auch den Steg 
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‚ls Platte, so ist die Definition der tragenden breite 
ınd die Übergangsbedingunz eindeutig. Behandelt 
nan dazegen den Ster nach der elementaren Bie- 
ınestheorie. so gibt es verschiedene Möglichkeiten, 
son denen hier die einfachste wiedergegeben sei. 


Wir führen zunächst mit Kärmän den Begriff 
ler tragenden Breite: 


h 
\oxdy 

b; . v (7) 
Ox)y=V 


oder die damit zusammenhängende reine Zahl: 


h 
\ Ox dl Y 
F, bö 6 
‚ F, hd hdlo)y—o 


ein. 


Es ist nun. wenn M (.r) die Momentenfunktion, 
welche auf ein Rippenfeld entfällt, bedeutet, nach 
der elementaren Biegungstheorie: 


a 6M 1 1) (S) 
x)/y=\) F, h 1 j I 4 ’ . . ° . ° . 
oder 
h 
\ Ox d Y. 
(ox) | 1 a. A 
y-0 hd (ox)y=0 F, h' 


also: 
bh 
1 4öf 6M 
> (x), N) 1 h ” | ÖOx d Y - F, h . . . (8a B- 


u 


wo der Faktor '/, wegen der Mittelbildung über 
die Rippenbreite hinzutreten muß, wenn man die 
(+l. (2) und (6) benutzt. 


Diese Bedingung stetigen Spannungsüberganges 
vom Stee auf den Gurt liefert die einfachste Be- 
ziehung,. Man muß dann, wie wohl zulässig, die 
geringe Unstetigkeit im Dehnungsübergang (infolge 
der Querkontraktion in der Gurtplatte) vernach- 


lässigen. 


Da beide Seiten Funktionen von allein sind. 
können sie beide mit der Periode gleich der Spann- 
weite Z in Fouriersche Reihen entwickelt und 
die Koeffizienten der linken Seite durch Vergleich 
paarweise bestimmt werden. 

Damit ist die gestellte Aufgabe grundsätzlich 
zelöst. Die noch erforderliche Durchführung von 
Zahlenbeispielen soll an anderer Stelle gegeben 
werden. 


Charlottenburg. H. Reißner. 460 


I, Formel (8) ergibt sich wie folgt: 


(Fe + Fı)z=F,; = ‚, Fr =bJ, Fa=hd, u p, 
Schwerpunktsabstand =; n " a 
Widerstandsmoment - Fız+F5 = + [= :)”) : 

>49), Ox -T 


Formänderung und Beanspruchung 
dünner kreiszylindrischer Ringe und 
Rohre. Bei Berechnungen von Maschinen und 
Baukonstruktionen treten des öfteren Fälle auf. 
bei denen sich der betrachtete Konstruktionsteil 
mit enüzender Genauirkeit als ein dünner kreis 
zylindrischer Ring ansehen läßt (umlaufende Ringe: 
Behälter). Ferner kommt der Fall vor. daß z. B. 
ein Maschinenteil als ein aus solehen Ringen zu- 
sammengzesetztes Gebilde aufzefaßt und seine Be- 
rechnung auf die Ringeberecehnung zurück geführt 
werden kann. Wichtige ist dabei immer besonders 
die Kenntnis der Formänderunzen,. Momente und 
Auflagerkräfte an den Rändern. 




















e kalcm auf dem Umfang 
M cmkajem auf dem Umfang 
7 
= Xel & 7 7 
122 i 
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—— er 0 ! 
RI Ar n 
29 N — {} l 
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— nt / TB | 
DA *® 
Y, I FE | I 
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Abb. 1. 


Die Gleichgewichtsbedineunze für ein Element 
eines solchen Ringes (s. Abb. 1). 


dQ os-I2-.dg—+yqyr.dse.dpyg.. (W 
führt auf die Differentialzleichung 


A -) +) 
ay_isy__ 129 -().. (>) 


dz“ 'r252 E s3 0 u - . . a). 


Vorausgesetzt ist dabei. daß s klein gegenüber r 
ist. Die Lösung erhält eine für die praktische 
rechnerische Verwendung zut brauchbare Form. 
wenn man hyperbolische Funktionen heranzieht. 
Das allgemeine Integral ergibt sich aus dem parti- 
kulären Integral 

















yr? 
Y 4 z (3) 
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Abb. 2. Bezogene Ausbiegung bei <-- I unter Belastung 
durch p, Rand s=#% frei. 
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Math. und Meeh. 
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/} 4% 74 oo Kol? +(0s?& 
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Abb. 3. Bezogene Ausbiegung bei < l unter Belastung », 
Rand <= vollkommen eingespannt. 
3115 
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I k Sin? + sın?a 
! 1 Ye I —— 
7; Sina Rolx-sına casa 
d Hl 
4 
Bea 7 
4 5 
MM Ay47% 05 j ... m DPF E 7 Be 
f: 0 Bei 


bu ” 
5 


Abb. 4. Bezogene Ausbiegeunge bei < l unter Belastung 
dureh p, Rand <= 0 gehalten. 















IH Rand & :0 rer: 
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130 Hd" Sara na HN 
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. ;  „Iindarsum:a 
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| Ah 
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\bh.>. Bezogener Winkel bei < I unter Belastung dureh » 
und bezogene Ausbieeunge beis= | unter Belastung dureh M. 
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\bb.6. Bezogener Winkel-bei,s = | unter Belastung dureh M. 





und dem allgemeinen Integral der reduzierten 
Gleiehune dann in der Form 


Y Ucosaöstojas—(Ü,sinaZfoja: | 
(4) 


’ u 2 - l v * FE 2} ” yr? 
(„cosasemnai- U‚sinas£emafi- 4 | 


“ Es 


mit a j 3//Yrs und den aus den Randbedin 
gungen zu bestimmenden Festwerten (,....(,. 
Aus dieser Gleichung der elastisch verformten 
Meridianlinie lassen sich die Umfangs- und Längs- 
spannungen für jede Stelle berechnen. Bestimmt 
man die Inteerationskonstanten für die drei Be 
lastungsfälle 


I. Belastung am Rande Z=1 dureh Randkraft 
p je Längeneinheit des Umfanges, 

II. Belastung durch Innendruck g je Flächen 
einheit, 

Ill. Belastung am Rande £== 1 durch Randmoment 
M je Längeneinheit des Umfanges, und für die 
verschiedenen Kinspannungsmöglichkeiten 
des Rohrendes 2=0, so bekommt man die 
folgenden Ergebnisse für die radiale Aus- 


j “ dy 

bieeune y und den Winkel der Meri- 
G=z 

dianlinie mit der Achse am Ende € l in 


Abhängigkeit von der die Abmessungen 
zusammenfassenden kennzeichnen 
den Größe. 

l. Belastung dureh p». 


l. Rine bei & 0 frei: 


‚v„'yr pZimastoja—sinacosa 
4 Z «) . = mp Pr - - 4 
sV s E Sin’a —sin’a 
» Zın? Mn sın? yı 
7 > | a r p oım”a—sSın’a 
s? E Sina —sin?’a 
2. Ring bei <=0 vollkommen eingespannt: 
„v, ryr p Sinastofa—sinacosa 
y=2y3 = nn, 
sVs LE ol? a—- cos? a 
r p Zin?a--sin?a 
Y 23 


s? E Soj?a-- cos?a' 
3. Ring bei 2<=0 aufliegend: 


os trVr p Sin? a--sin?’a 


sYs E Sina Soja — sinacosa' 
r » Sinaftoja--sinacosa 

ZinaNtoja -sinacosa' 
Il. Belastung dureh g. 


r De Br rg ' 
I. King bei &=V frei: Y= 5: Y=d. 
sE 


2. Ring bei £=0 vollkommen eingespannt: 


' 


r? q (Koja — cosa)? 


Y r 3] . ’ 
s E Stoj’a + cos’ «a 

„Yo rYr g sinasoja— cosakin« 
71 2 ] 5) . 


s] S E Ktoj? (L 7 eos? ( 
3. Rinz bei 2=0 aufliegend: 


r? 4 (Sina sine) (Noja -cosa), 


Y ; =; - - 

sE ZinaNoja - sindeosa 
to trvVrg sınacına 
Y -) >» 


sy s E Einakoja— sinacosa 
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Ill. Belastung dureh M 


1. Rinz bei 2=0 frei: 
2V3 r MSin®a--sin’a, 
U ” «) . . nz Ps . > “ 
7 s’ E Siın?a — sın?a 
ı 433 |r ME&inaektoja—sinacosa 
7 _ eo) . — r . 5) . 
s®yVs E Sin?’a — sin?’a 
2, Rinz bei <=0 vollkommen eingespannt: 
‚r MSina-sin’a 
’ - 14 | ie) 1} v } I „2 j 
s? E toj? a — cos’ a 
1/33 |r MEinasoja- sinacos«a 
'Y J ' Tr - . 
s®Vs E to)” a —+-cos’a 
3. Ring bei < 0 aufliegend: 
‚—=2y3 r MSinastoja--sinacosa 
! pe «) > Zus - ° . 
s? E Sina Mo) a — sin acosa 
Wr r M Zin?a— cos?a 
aut | 


s? Vs E Sina Soja — sina cos: 


Außer bei den Gleichungen für y im Falle II 2 
and 113 gilt folgendes. Die in diesen Ausdrücken 
auftretenden Funktionen von «@ nähern sich, wie 
man leicht feststellen kann. mit wachsendem « 
dem Grenzwert eins. Für € =x, d.h.. dar und s 
nicht gleich null sein können, für /=&®, sind 
daher 4 bzw. % zleich dem von a freien Teil der 
Ausdrücke. Diese von « freien Glieder bedeuten 
also Ausbiezung Yı—=«» und Neigung =» am 
linde = I für unendliche Rohrlänge. Die den 
Kennwert « enthaltenden Funktionsausdrücke sind 
für die Fälle I und Ill in den Schaubildern dar- 
vestellt. 

Bei dem längs des Randes 2=0 vollkommen 
eingespannten Rohre ist an diesem Rande das 
Monent M, je Längeneinheit des Umfanges 


I. bei Belastung durch p: 


Yrs cosa Sina --sinastoja 


j B P Noj’a + cos?a 


Mo 


Il. bei Belastung durch g: 


rs Sin?a + sin?a 
MH. == 4 = —; 
> B Ktol?a —- cos? a 


Ill. bei Belastung dureh M: 


y 1 cos aNv) «a 
77T Kof?a-+ cos?u 


und die Auflagerkraft >, je Längeneinheit des 
Umfanges 


Il. bei Belastung durch p: 


a cos. aNoj a 
Po = n J fo | Di) ; 
Soja -—+- cos’a 


Il. bei Belastung durch g: 


I\rs Sinaftoja+-sinacosa 
) 


Po .. N 6) 
v3 oO]? a —+- cos” a 


Ill. bei Belastung durch M: 


2V 3 „,cosa Sina — sinastoja 


= \ Pen — 
Pu Vrs Koj?a— cos?’a 


Beachtlieh ist, daß bei dem freien am Rande & | 
durch » belasteten Ring das Verhältnis der Um- 
fangsspannung o, am belasteten Rande zu dem 
Werte ogr. der sich bei gleiechförmiger Verteilung 
der Spannung ergeben würde, mit abnehmendem « 
dem Grenzwert 4 zustrebt und ihm sehon für mäßig 
kleine @« sehr nahe kommt. Aus dem Ausdruck für 
diesen Verhältniswert 


0 Sina oja — sina cos a r 
- 2a (5) 
gl zu a sın’ a 


erhält man nämlich. wenn man die Entwicklung 
der Funktionen von a nach Potenzen von @ benutzt 
und von jeder Reihe die ersten vier Glieder bei- 
behält, nach einigem Umformen 


0 1 I \25 
| l | a FR „a F r = .) a” (6) 
Gl 1 d’rı)*+d dtrı)*+ dh od), 


oder, wenn man /, r und s einsetzt: 


0, \ .) i* n I ’ 
4 -+ 0,1143 -—_— — 0,0026 . (6a). 
Ol y. r’ 5° 


Berlin-Charlottenburg. 


Max Jobst Gereke. 147 


Die Berechnung von Eigenschwingungen 
einer gleichseitig dreieckigen Membran 
mit Hilfe von Sechsecksnetzen. |)ıas )iffe- 
renzenverfahren hat in den letzten Jahrzehnten eine 
dauernd steigende Beachtung bei der Lösung ellip- 
tischer Randwertprobleme gefunden, weil es gegen- 
über anderen Lösungsverfahren die Vorteile besitzt. 


I. schon nach verhältnismäßig kurzer Rechnung 
einen guten Überblick über den Verlauf der lL,ösungs- 
funktion zu zeben, und 


2, viel anschmiegungsfähiger an die Randbedin 


eungen zu sein). 


Der zweiteenannte Vorteil tritt noch in weit 
stärkerem Grade in Erscheinung, wenn man nicht 
nur rechtwinklige Gitternetze verwendet, sondern 
jeweils dem betreffenden Problem angepaßte Netze. 
Ich möchte im folgenden zeigen, wie man durch 
Verwendung von „Sechsecksnetzen” Probleme be- 
handeln kann, die ohne Benutzung des Differenzen- 
verfahrens der mathematischen Behandlung, etwa 
mit Hilfe von Reihenentwieklungen. kaum zugäng- 
lich erscheinen. Als Beispiel werden die Eigen- 
schwingungen einer am Rande eingespannten Mem- 
bran der konstanten elastischen Kraft « und der 
konstanten Dichte o betrachtet. die im Ruhezustand 
ein eleichseitizres Dreieck D der Seitenlänge «a in 
der r-y-Ebene bedecke. Bedeutet v(r,y) die Maximal- 
ausbieerung an der Stelle z,y. so ist vw(r.y) Lösung 
des Eigenwertproblems 


"u, du 
#7 4 
02° 04° | 
(1). 
o i h 
= v?u u im Inrieern von D | 
„0 auf dem Rande von D . . 2.2...) 


Dabei sind die Frequenzen » noch zu bestimmen. 


I), Das zeigte schon ©. Runge: Zeitschr. f. Math. u. 
Phys. 56, 1908, S. 275 bis 232 bei Berechnung der Torsion eines 
Stabes von kreuzförmigem Querschnitt. 
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Zur Lösung dieser Aufgabe benutzen 
ecks- oder Sechsecksnetze. Um in diesen 
die entsprechenden Differenzenausdrücke 
zu erhalten. verwenden wir die bis 
höherer Ordnung in Ah richtige Formel 


Netzen 


auf Glieder 


(3). 


wobei / der Mittelwert der /-Werte auf einem Kreis 
mit dem Radius A um den Punkt r,, ist: 


hecosy.y+-hsingy)dyg. 


wir Drei- 


für Au 


\Y_ 


2, die Rechenarbeit beim Sechsecksnetz erheblich 
seringer ist als beim Dreiecksnetz, da in (4) weniger 
Koeffizienten auftreten als in (6) ?). 


Die Rechnung verläuft sehr einfach, wenn man 
nieht allzu feine Netze benutzt. Bei 16 inneren 
Punkten (Abb. 4) hat man für jeden der 16 inneren 
Punkte die Gl. (5) anzuschreiben: 


on, Us 0) | 
’ 


OMU, HU, U. u = | ee FR 


Diese 16 Gleichungen für die unbekannten «-Werte 
haben eine nichttriviale Lösung, wenn die Koefti- 
zientendeterminante verschwindet. Das ist eine 














4 

J 0 ’ 
LIYiH 2 
Abb. R 


Bei einem Sechsecksnetz (Abb. 3) entspricht dem 
Mittelwert das Mittel der Werte von f in den drei 
Nachbarpunkten, also bei einer Bezeichnungsweise 
der Punkte wie in Abb. 3 


und nach (3) ergibt sich 


| 
If >= f, + l; I l; 


Ba er 6 
3? Io “ 


(4]. (1) entspricht also die Differenzengleichung 


0 In f, I ? Je . . . . . . R . (D). 


wo A=,,,(3 0) gesetzt ist. 


Analog findet man für ein quadratisches Netz 


(Abb. 1) den bekannten Ausdruck 
| h 
1/7 > 7. 4 rhrhrm 4.) 
und für ein Dreiecksnetz (Abb. 2) 
I/ > Per l/ ] ls | l; 1 7 1 ft Te 6%) - (6). 
3? 


Das Sechseceksnetz ist im allgemeinen günstiger 
als das Dreiecksnetz (wie sich auch bei unserem 
Beispiel bestätigte), weil bei dem gleichen Problem 
und bei gleicher Anzahl innerer Punkte. also bei 
gleicher Unbekanntenanzahl 


l. die zugehörige Maschenweite A beim Sechsecks- 
netz kleiner ist als beim Dreiecksnetz. also bessere 
Näherungswerte erzielt werden, 





Abb. >. 








Abb. B. 


Gleichung 16. Grades für o, aus der man Nähe- 
rungswerte für 16 Eigenwerte (nicht notwendig die 
ersten 16 Eigenwerte) erhält. Die Rechenarbeit 
läßt sich und das ist für größere Anzahl innerer 
Gitterpunkte wichtig sehr vereinfachen durch 
folgende erlaubte Ansätze: 


a) Bei dem „symmetrischen Ansatz” v;= uw;’ hat 
man nur 10 Unbekannte, also eine Gleichung 
10. Grades zu lösen, von der man einige Wur- 
zeln aus 


b) dem „vollsymmetrischen Ansatz“ kennt, bei 
dem man die Funktionswerte in allen symmetrisch 
gelegenen Punkten gleichsetzt, also z. B. in Abb. 4: 

U, — U, n 


u le ZUM. 





Ara 





77 








Abb. 4., 


2) Dreieeksnetze hat bereits H. Mareus in seinem 
Buche: „Die Theorie elastischer Gewebe und ihre An 
wendung auf die Bereehnung biegsamer Piatten“, Berlin 
1932, verwendet. 
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Dann sind nur 5 Funktionswerte %, Ur, U, Uz, Ug 
ınbekannt. man hat mithin eine Gleichung 5. Grades 
‚u lösen. Die Lösungen dieser Gleichung sind auch 
,ösungen der a) entsprechenden Gleichung. 


ec) Der .„antisymmetrische Ansatz” u; = ur. 
eb =uy—=wm0 liefert eine Gleichung 


, Grades. Die sich aus a) und e) ergebenden 
Werte o sind genau die 16 aus (7) zu berechnen- 
den o-Werte. Alle in a) bis e) genannten Glei- 
chungen lassen sich durch 0?=r in eine Gleichung 
von halb so hohem Grade überführen, so daß man 
z. B. bei b) nur eine quadratische Gleichung zu 
lösen hat. 


Den folgenden Zeichnungen wurden Berechnungen 


mit 49 inneren Gitterpunkten (also a—=8hy 3) zu- 
erunde gelegt. 


a) 4 > 


R i 120 
b)) „eh=- 
. a? 
i 0 
e) /. j 
’ a? 
h L 215 
dl ) ben N 6} 
a? 
b 300 
e) A, Ar —, 
a? 
2 330 
f ) Ay 10 a? 
Berlin. Lothar Collatz. 428 


(Aus dem Institut für angewandte Mathematik an 
der Universität Berlin.) 


Knotenlinien der ersten zehn Eigenschwingungen. 
Symmetrischer 
Typus Typus 


Antisymmetrischer 


a) 
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BUCHBESPRECHUNGEN 


Dr. EUGEN JAHNKE, weil. Prof. a. d. Berg- 
akademie in Berlin, und Dr. Dr. FRITZ EMDE, 
Prof. d. Elektrotechnik a. d. Techn. Hochschule 
Stutteart, Funktionentafeln mit Formeln und 
Kurven. 2.. neubearb. Aufl. (7. Tausend.) XVII 

330 8. m. 171 Textfig. Leipzig und Berlin 1933, 
Verlag von B. G. Teubner. Preis geb. 16 M. 


Die Funktionentafeln von Jahnke und Emde sind 
so bekannt und werden so vielfach benutzt, daß 
sie einer besonderen Empfehlung wirklich nicht 
bedürfen. Endlich ist jetzt die langersehnte neue 
Auflage herauszekommen, für die Herr Emde allein 
zeichnet. Das Buch ist in fast allen Abschnitten 
erweitert: vieles ist neu hinzugekommen. Das soll 
hier nicht im einzelnen aufgezählt werden; hervor- 
heben möchte ich aber die schönen zeichnerischen 
Darstellungen, insbesondere der komplexen Funk- 
tionen, die außerordentlich dazu beitragen, die 
komplexen Funktionen der Anschauung näher- 
zubringen. Wer im numerischen Rechnen erfahren 
ist. wird die Leistung von Herrn Emde und seinen 
Mitarbeitern bewundern; der Dank aller derer, die 
das Werk benutzen werden, mag Herrn Emde für 
die unglaubliche Mühe belohnen, die er im Dienste 


der angewandten Wissenschaften auf sich ge- 
nommen hat. Dank gebührt auch dem Verlag 
B. G. Teubner, der das Buch der Tradition getreu 
vorzüglich ausgestattet und trotzdem den Preis in 
mäbigen Grenzen gehalten hat. Trefftz. 


Dr.-Ing. ERNST LEHR VDIL Spannungs- 
verteilung in Konstruktionselemen- 
ten. Auswertung der bisherigen Forschungsergeb- 
nisse für die praktische Anwendung, im Auftrage 
des Vereins Deutscher Ingenieure ausgearbeitet. 
59 8. u. 179 Abb. im Text u. auf 9 Taf. Berlin 193 
VDI-Verlax G.m.b.H. Preis 7.50 M. 


Das Heft enthält eine wertvolle Zusammenstel- 
lung dessen, was bisher an experimentellen Metho- 
den zur Spannungsermittlung entwickelt worden ist 
und eine Darstellung der wichtigsten Ergebnisse. - 
Im ersten Kapitel werden nach einem einleitenden 
Abschnitt über den Spannungszustand und seine 
Veranschaulichung die experimentellen Methoden 
der Spannungsbestimmung beschrieben (Dehnungs- 
messer, Spannungsoptische Untersuchungen, Sonder- 
verfahren für Torsionsprobleme). - Im zweiten 
Kapitel sind die wichtigsten zahlenmäßigen Ergeb- 





3uchbesprechungen 


Ztsehr. f. angew. 
Math. und Mech. 
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nisse teils zeiehnerisch. teils in Zahlentafeln darge 
stellt. Die Wiederreabe Ist sehr übersichtlich. vor 
allem auf den beieerebenen neun Tafeln, die die 


Spannungsverteilungen für praktisch besonders be 
leutsame Fälle enthalten. Ein sehr umfangreiches 
L,iteraturverzeichnis ergänzt (diese Darstellungen. 
Im Vordererunde der Betrachtunzen stehen natur 
remäß die zweidimensionalen. d. h. die ebenen und 
die  rotationssyvmmetrischen Probleme. Wenn das 
Heft (die anrerende Wirkung hat, die es erzieien 
möchte, so ist zu hoffen, daß auch die Untersuchung 
ler ungleich sehwierizeren. eigentlich dreidlimensio 
nalen Fälle die erwünschten Fortschritte macht. 


Der Verfasser und seine Helfer haben sieh nicht 
umsonst erobe Mühe zereben. Nicht nur der prak 
tische Ingenieur wird aus ihrer Arbeit Nutzen 
‚\ehen. auch die rechnende Festiekeitslehre (deren 


l,eistuneen,. wie mir scheint, in den Ausführungen 
les Verfassers etwas zu schlecht weekommen) wird 
aus (dem reichen Anschauunesmaterlal Gewinn 
haben. Trefftz. 46 


Dr..Ingz. Dr. phil. R. POZDENA, Meter und 
Niloreramm. Entstehung und Sicherung des 
internationalen metrischen Maßsystems. Mathem.- 


Phys. Bibl.. Reihe 1, Bd. 76. 45 S. m. 15 Fig. 
leipzig und Berlin 1934, Verlag B. G. Teubner. 
l’reis kart. 1.20 M. 


Das Buch verfolet, wie der Verfasser im Vorwort 
das Ziel. .„.das Zustandekommen und die 
Sicherung des metrischen Maßsystems zu schildern, 
las von den Bewohnern der ganzen Kulturwelt täg- 


saot,. 


lieh gebraucht und vielfach als ‚ganz einfach‘ und 
fast „selbstverständliceh' hingenommen wird“, 
Die erste Hälfte des Buches enthält eine histo 


rische Darstellung der Entwicklung des Meßwesens. 
Nach (der Darstellung der frühesten Ansätze zur 


Normalisierun®e der Maße und der Versuche, ein 
physikalisch zegebenes „Naturmaß“ (z. B. in der 
länge des Sekundenpendels) zu finden, wird über 
lie wiehtiesten Gradmessuneen und deren Aus- 


wertune®e berichtet 
les metrisehen 
Hälfte (les 


und die schließliche Einführung 
Systems geschildert. Die zweite 
Buches handelt von der Sicherung (des 


metrischen Mabsystems, d. h. vor allem von («der 
/urückführun®e «der Meterlänere auf Liehtwellen- 


läneen mittels interferometrischer Methoden. 


Das anreerend weschriebene Büchlein, das auf 
wenieen >eiten eine Fülle des Wissenswerten ent- 
hält. kann sehr empfohlen werden. 


Trefftz. 462 


Dr. Dr.-Ing. GEORG SCHEFFERS, o. Prof. a. d. 


Teehn. Hoechseh. Berlin. Wie findet und 
‚elehnet man Gradnetze von Land 
und Sternkarten? Mathem.-Phvs.  Bibl.. 
Reihe 1. Bd. 85/86. 98 S. m. 12 Taf. u. 27 Abb. 
leipzie und Berlin 19534 Verlage B. G. Teubner. 
l’reis kart. 2.40 M. 


Über (die Methoden. nach denen die Gradnetze der 


l,andkarten hergestellt werden. herrscht im allge- 
meinen eine erstaunliche Unkenntnis. Auch in der 
\usbillun®e «der Lehramtskandidaten findet «dieses 


Kapitel nieht immer die Beachtung, die ihm wegen 
der gebotenen Möglichkeit der Verbindung von Uni 
versitäts- und Sehulmathematik zukommt. Das 
Büchlein von Scheffers, das er der Petrischule in 
leipzig als «deren ehemaliger Schüler zu ihrer 
Hundertjahrfeier gewidmet hat, behandelt die Kon 


struktion «der verschiedenen Kartennetze in sehr 
anregrender Weise, Die Gesichtspunkte (Flächen 
treue, Winkeltreue u. del... naeh denen ein Gral- 


netz konstruiert wird, werden besprochen und die 
Konstruktionen im einzelnen angegeben und aus 
eführt. Zahlreiche historische \usführungen 


lem man eine 
. besonders reiz- 


lieses 


machen die Lektüre Buches, 
erobe Verbreitunz wünschen möchte 


voll. Trefftz. 462 
KARL MENNINGER, Zahlwort und Zif 
fer. Aus der Kulturgeschichte unserer Zahlsprache. 


unserer Zahlsehrift und des Rechenbretts. X + 364 5. 


Breslau 1034. Ferdinand Hirt. Geh. 7 M. geb. 9 M. 


Wer als Praktiker oder Theoretiker der Wirt 
schaft oder Technik die Zahl ständig gebraucht, ist 
sich oft gar nicht bewubt, welche gewaltige kultu 
relle Leistung die nach und nach gewachsene Kunst 
les Rechnens bedeutet. Der Mathematiker findet 
hier ein zeeienetes Thema für allgemeinverständ 
liche Vorträge bei besonderen Gelegenheiten. Es sei 
an die schöne Rede erinnert, die der einstige Ver 
treter der angewandten Mathematik an der Univer 
sität Göttingen, Karl Runge, bei der Jahresfeier 
der Universität 1915 gehalten hat: Mathematik uni 
Bildung (Jahresbericht der Deutschen Mathematiker 
vereinigung, Bd. 24, S. 400 bis 415). Ref. hat öfter 
in Festreden über solche Dinge gesprochen. Eine 
reichtliebende Quelle für derartige Vorträge bietet 
(las anzuzeigende Buch. Mit einem staunenswerten 
Fleiß und Geschick hat der Verfasser aus gedruck- 


ten Werken und aus vielen Museen des In- uni 
\uslandes zusammengetragen, was sich über Zahl 
zeichen. Zahlwörter und das Rechnen mit den vier 


einfachen Recehnungsarten sagen läßt bis zur Ein 
führunz der „neuen”“ Zahlscehrift, „der indischen. 


dem endeültigen Siez eines fremden Kulturgutes 
nach tausendjähriger Wanderung, aber auch «dem 


Sier des Geistes, den er mit der reinen Stellenoril 
nune in der Janzen Geschichte der Zahlschrift end 


lich erlangt“. 170 schöne Abbildungen ergänzen 
den Text. Sehr praktisch ist die Art, wie auf sie 
auch an späteren Stellen verwiesen wird: (B 87 
d. h. Bild auf S. 87. Verfasser hat nieht mit der 


„mathematischen Brille“ die Entwicklung rein ge 
danklich betrachtet: sieht doch diese Brille nach 
seiner Meinung nicht immer kulturgeschiechtlich 
richtige. Er vermeidet daher auch die üblichen Kunst 
wörter, wie System, dezimal u. a. und gebraucht 
nur deutsche Ausdrücke, wie Zehner-Stufung, Zwan 
z\eer-Bündel usw. Mit besonderer Vorliebe ist er 
den Zahlen- und Geldnamen nachgegangen. Das auf 
S, 260 erwähnte Dollarzeichen, das er als vielleicht 
aus dem ]b,. dem Zeichen des Pfundes &, entstanden 
angibt. ist nach der eingehenden Untersuchung des 
kürzlich gestorbenen amerikanischen Historikers (er 
Mathematik Cajori aus p* (peso) entstanden (The 
Kvotution of the Dollarmark. Popular Seience Monthly 


Dez. 1912). Für eine neue Auflage möchte Ref. 
empfehlen, dem Leser nieht mehr M. Cantors 
Vorlesungen über Geschichte der Mathematik als 


„auszezeichnetes Quellenwerk“ zu nennen: neben 
der mit Recht genannten zuverlässigen Geschichte 
der Elementarmathematik von Tropfke wären 
die schönen Bändehen 1 u. 34 der Lietzmann-W it 
tingschen Mathematisch-physikalischen Bibliothek 
anzugeben: Löffler, Ziffern und Ziffernsvsteme. 
Kin ausführlicher Sachweiser erleichtert (den Ge- 
brauch des auch dem Leserkreis dieser Zeitschrift 
sehr zu empfehlenden volkstümlichen Buches. 
Frankfurt a. M. W.Lorey. 431 


. STEINITZ und H. RADEMACHER, \or- 


lesungen über die Theorie der Foiyeder. 


Die Grundlehren d. math. Wiss. in Einzeidarstell. 
mit besonderer Berücksichtigung d. Anwendungs- 
veb. Hrse. von R. Uourant. Gemeinsam mit 
\W. Blaschke. M. Born u. B. L. van der Waerien. 
bi. 41.) Berlin 1934, Julius Springer Verlag. XIl 
351 S. Preis zeh. 27 M, geb. 28.80 M. 
In der bekannten, meisterhaft klaren, strengen 


und eeeistreichen Darstellung, die alle Schriften von 
Steinitz auszeichnet und (die sieh in gleicher Weise 
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n denvonRademacher stammenden Paragraphen 
vorfindet. wird eine systematische Theorie «er kon- 


vexen Polyeder des dreidimensionalen Raumes ge- 


reben. Sie gipfelt in dem „Fundamentaisatz der 
konvexen Typen“: Jedes K-PolyeJler kann 
ıls konvexes Polyeder realisiert 
werden. Ein K-Polyeder ist ein rein kombina 
torisch aus „Ecken“, „Kanten“ und „Flächen“ dureh 
\ngabe der wechselseitigen Inzidenzen aufgebauter 
Komplex, der gewissen einfachen Bedingungen (z. B. 
‚usammenhängend zu sein und die Eulersche Cha- 
rakteristik der Kugeltläche zu haben) genügt. Zum 
sleichen Typus gehören zwei Polyeder, wenn sie 
lenselben kombinatorischen Aufbau besitzen. 
Nach einem historischen Überblick über die 
olyedertheorie wird der Begriff des Polveders axio- 
matisch kombinatorisch deiiniert, wobei der Funda 
mentalsatz der Flächentopologie kombinatorisch be- 
wiesen wird. also auch Polyeder höheren Zusammen- 
hanzes abgehandelt werden, ohne auf ihre eben- 
tlächige Realisierbarkeit Rücksicht zu nehmen. Im 
Iritten Abschnitt werden drei Beweise für den 
"undamentalsatz der konvexen Typen gegeben. Bei 
len rein geometrischen Beweisen wird Wert darauf 
releet, die Theorie der konvexen Polyeder allein 
auf Grund der Axiome der Verknüpfung und An- 
ordnung aufzubauen. Zum Beispiel wird als Hilfs- 
satz allein mit diesen Axiomen «er ‚Jordansche 
Kurvensatz für Polygone der Ebene und der ent- 
sprechende Satz für ebentlächige Polveder des eukli- 
dischen Raumes bewiesen. Als weiterer Hilfssatz 
wird gezeigt. daß ein einfaches (d. h. sieh nicht 
selbst durehdringendes) ebentlächiges Polyeder des 
euklidischen Raumes stets orientierbar ist, die pro- 
jektive Ebene sich also nieht ohne Selbstdarehdrin 
rungen als ebenflächiges Polyeder realisieren läbt. 
- An den dritten Beweis schließt sich die strenge 
Begründung dafür an, daß je zwei Polyeder von 
rleiehem kombinatorischen Aufbau (Typus) stetiz 
ineinander überführbar sind und daß es ck ver- 
schiedene Polyeder von vorgegebenem Typus gibt: 


% ist die Anzahl der Kanten, und kongruente 
Polveder werden als nichtverschieden betrachtet. 
H. Seifert. 472 


Ferner sind bei der Schriftleitunz folgende 
Bücher eingegangen (ausführliche Besprechung 
bleibt vorbehalten): 


X. L. BOWLEY, Prof. a. d. Universität London. 
Grundzüge der mathematischen Öko 
nomik. Beiträge zur quantitativen Wirtschafts 
forschung herausgegeben von Eugen AUTSCHUL. 
Band I. 128 S. m. 19 graph. Darstellungen. Leipzig 
1934, Verlag Hans Buske. Preis geh. 8 M. 


HEINRICH LIEBMANN, o. Prof. a. d. Un. Heidel 
berg. S ynit hetische Geometrie, Teubners 
mathem. Leitfäden Bd. 40.) 109 S. m. 45 Fig. 
Leipzig und Berlin 1934, Verlag B. G. Teubner. 
Preis kart. 5.60 M. 


Aerodynamie Theory, A General Review of Pro 
eress, Under a grant of the Guggenheim Fund for 
the promotion of aeronauties. hrsg. v. William Fre 
deriek DURAND. Vol. I: Mathematieal Aids (W. 
F. Durand), Fluid Mechanies (W, F. Durand und Max 
M. Munk), Historical Sketch (R. Giacomelle und E. 
Pistolesi,. 394 S. m. 151 Fig. Berlin 1934, Verlag 
Julius Springer. Preis geb. 20 M. 


Dr. HELMUT HASSE, o. ö. Prof. d. Math. a. d. 
In. Marburg. Aufgabensammlunge zur 
höheren Algebra. (Sammlg. Göschen Bid. 1082.) 
175 S. Berlin und Leipzig 1934, Verlag Walter 
de Gruyter & Co. Preis geb. 1.62 M. 


Dr.-Ing. E. h. WILHELM NUSSELT VDl, 0. 
Prof. a. d. Techn. Hochsch. München, Tech 
nische Thermodynamik. (Sammlg. Göschen 
Bd. 1084.) 144 S. m. 65 Abb. Berlin und Leipzig 
1934. Walter de Gruyter & Co. Preis geb. 1,62 M. 


NACHRICHTEN 


Paul Ludwik f. 


\m 28. Juli ist in jähem Tode in verhältnismäbig 
juneen Jahren ein Mann dahingegangen, der auf 
dem Gebiete (der Metall- und Werkstoffkunde 
Bleibendes geleistet hat. P. Ludwik war zu Schlan 
in Böhmen im Jahre 1878 geboren, studierte an (der 
Teehnisehen Hochschule in Prag und kam dann als 
Assistent an die Technik nach Wien, wo er 1904 
das Doktorat erwarb und sieh schon ein Jahr dar- 
auf habilitierte. 1910 wurde er außerordentlicher. 
1918 ordentlicher Professor für mechanische Techno- 
logie und Materialprüfungswesen.  Berufungen an 
die Teehnische Hochschule in Charlottenburg und 
an «das Kaiser-Wilhelm-Institut für Metallforschung 
lehnte er ab. da er sieh nieht von Wien trennen 
wollte. 1923 wurde er zum ordentlichen Mitglieil 
ler Wiener Akademie der Wissenschaften gewählt. 

(roß sind seine Verdienste um (die wissenschaft- 
liche Metallkunde. In seinem 1909 erschienenen 
Büchlein „Elemente der technologischen Mechanik“ 
behandelt er hauptsächlich vom experimentellen 
Standpunkt die allgemeinen Gesetze des Metall- 
tlusses beim Walzenzuer und stellt die Fließkurve 
eines Werkstoffes dar. Wichtige sind seine Ver- 
suche über das Altern von Stahl, über Ermüdungs- 
voreänge und Dauerfestiekeit. Auch um die Klar- 
stellung des Härtebegriffes hat er sich verdient ge- 
macht, die Kugeldruckhärteprobe stammt von ihm. 

Trotz «dieser reichen Forschertätiekeit war er 
nicht einseitig, sondern ein Mensch von allgemeinen 
Interessen, wie die Herausgabe einer kleinen philo- 


sophischen Schrift bezeugt. Er war beliebt bei 
seinen Schülern und Kollegen und dabei von einer 
außberordentlichen Bescheidenheit. obwohl ihm äußere 
Auszeichnungen in hohem Maße zuteil wurden, 
Freunde und Fachgenossen werden ihm ein ehrendes 
Andenken bewahren. K. Wolf, Wien. 461 


Gesellschaft für angewandte Mathematik u. Mechanik. 
Hauptversammlung in Bad P’yrmont. 

Programmäßie fand am 14. und 15. September 
die Hauptversammlung der Gesellschaft in Bad P’yr- 
mont,. im Rahmen des Mathematiker- und Physiker 
tages statt. Gegenüber dem im Heft 4 veröffen! 
lichten Programm sind die folgenden Veränderungen 
eingetreten: Die Herren G. Vogelpohl. 
W. Müller und ©. Weber waren an der Teil 
nahme verhindert. so daß ihre Vorträge ausfallen 
mußten. Statt dessen kamen hinzu die nachträg- 
lich gemeldeten Vorträge: R. Hermann. Aachen. 
Zur Zähigkeit wässeriger Lösungen (Freitag nach- 
mittag); H. Schliehting. Göttingen. Kanal- 
einlaufströmung: E. W. Spannhake. Karlsruhe. 
Neuere Kavitationsversuche am Massachusetts In 
stitute of Techn. in Cambridge USA. und ihre Deu- 
tune (Sonnabend vormittag). Die Vorträge waren 
erfreulich gut besucht. das Interesse der Teilnehmer 
äußerte sich in lebhafter Diskussion. 

Die Mitzrliederversammlung fand am Sonnabend. 
15. September, vormittags 9 Uhr, unter Leitung (des 
zweiten Vorsitzenden, Herrn Trefftz, statt. Zu- 
nächst wurde der im Laufe des Jahres verstorbenen 
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Zuschriften an den Herausgeber 
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Dann erstattete 
Vertretune des Geschäftsführers 

Ks wurden die folgenden Be 
Vorstand wurde 
ebenso (die Kassenprüfer A. Betz 
(öttineen) und H. Fromm (Berlin). An 
saıtzuneseemäß ausscheidenden Berren 0 Föpp!| 
Braunschweie), R. Grammel (Stuttgart). E. Meiß- 
ner (Zürich) wurden die Herren Th. Pösch] 
Karlsruhe) und ©. Biezeno (Delft) in dden wissen- 
schaftlichen Ausschuß eewählt. Der Vorstand wurde 
bevollmächtiet, den Eintritt der Gesellschaft in den 
RTA und in den Reiehsverband der Mathematiker- 
vereine zu vollziehen. Der Jahresbeitrag wurde un- 
verändert zu 3 RM festgesetzt. Ort und Zeitpunkt 
der nächstjährigen Tagung sollen vom Vorstand im 
Einvernehmen mit den physikalischen Gesellschaften 
und der Deutschen Mathematikervereinierunz festge- 
setzt werden. 


Mitzlieder der Gesellschaft gedacht. 
Herr Trefftz in 
den Jahresbericht. 
eelabt: 

wiedereewählt. 


Der bisherige 


sehliüsse 


Auflösung des Deutschen Verbandes Technisch- 
Wissenschaftlicher Vereine. 

Der Deutsche Verband Technisch-Wissenschaft- 

licher Vereine, dem auch die Gesellschaft für ange- 


Stelle der 


wandte Mathematik und Mechanik angehörte, hat 
semäbß der in Nr. 27 der RTA-Nachrichten veröifent 
lichten Erklärung am 1. August d. J. offizieil seine 
Auflösung beschlossen. Sämtliche Arbeiten 
Deutschen Verbandes werden von der Reichsremein 
schaft der technisch-wissenschaftlichen Arbeit (RTA 
weitergeführt. 


(les 


Nicht von der Auflösung betroffen sind die selb 
ständigen Ausschüsse und sozen. angeschlossenen 
Orreanisationen des Deutschen Verbandes (Deutsches 
Nationales Komitee der Weltkraftkonferenz. der 
Deutsche Ausschuß für das Schiedszerichtswesen. 
die Siemens-Ring-Stiftung. die Technisch-Wissen 
schaftliche Lehrmittelzentrale und die Gesellschaft 
zur Herausgabe der Illustrierten Technischen Wörter- 
bücher). 

470 


Persönliches. 


Hf. Priv.- Dozt. Dr. R. Iglisch ist beauftragt 
worden. an der Techn. Hochschule in Aachen 
die Mathematik in Vorlesungen und Übungen zu 
vertreten. 


ZUSCHRIFTEN AN DEN HERAUSGEBER 


Biegung krummer Stäbe. Ileır G.Lindner 
untersucht in einem in Bd. 14. Heft 1 dieser Zeit- 
schrift erschienenen Aufsatze die Spannungsvertei- 
lunz in krummen Stäben unter drei verschiedenen 
Annahmen über die @Querschnittsverformung. Zu 


der in der Einleitunz zerrebenen Berründunz und 
zur >Sehlußfolzeerunze möchte ich folgendes be- 


merken: 

Herr Lindner hält es für eine Schwäche der 
iibliehen Grashofsehen Annahme, daß bei 
seharfwinkliger Krümmung (diese Annahme zu & 
eroßen Spannungen führt, weil dieses Ergebnis 
mit der Erfahrung nieht im Einklange steht. Wir 
dürfen aber nieht außer acht lassen, daß die auf 
(Grund dieser Annahme abgeleitete lineare Span- 
nunesverteilunz sieh auf rein elastische, dem 
Hookeschen Gesetze unterworfene Verformungen 
bezieht. (Auch die beiden anderen von Herrn 


,indner untersuchten Fälle beziehen sich auf 
solehe Verformunzeen.) Nach der strengen 
Klastizitätstheorie sind aber zu- 


rückspringende Ecken Unendlich- 
keitspunkte (Th. WyvBß. Die Kraftfelder in 
festen elastischen Körpern. S. 94). Daß in der 
Wirkliehkeit in den Eekpunkten bei mäßiger Be- 
lastune keine KEinrisse entstehen, das ist dem 
Spannungsausgzleiche infolge plastischer Verfor- 
mung umd (der Stoffstruktur zu verdanken. Übri- 
vens sind für das ebene Problem strenge Lösun- 
ven bekannt (Handb. der Phys. Bd. 6. S. 184). 
So bleiben im Falle «der reinen Biegung des Kreis- 
rinesektors in Übereinstimmung mit der Gras- 
hofscehen Annahme die Querschnitte eben, und 
die Spannungen werden beim Grenzübergange zum 
Sektor mit wo kleinem inneren Halbmesser im 
Krümmungsmittelpunkt ®. 


Die im Schlußabsatze abzeleitete Folgerung aus 
den enerzetischen‘ Betrachtungen, wonach die 


parabolische Krümmung. d. h. die lineare Span- 
nungsverteilung von (den drei untersuchten Fällen 
der streng theoretischen am nächsten käme, ist 
nieht riehtiz. Wir können nämlich die drei Fälle 
vom Gesichtspunkte der Variation des Spannungs- 
ustandes betrachten. wie es W., Hovezaard ge- 
tan hat bei der Untersuchung des Spannungszu- 


standes vom Balken (.„ZAMM“ Bd. 10. 8. 15). 
Dann hat aber nieht die lineare 
Spannungsverteilung mit dem größ 


ten Verdrehunegswinkel %. sondern 
die Grashofsche. hyperbolische mit 


dem kleinsten Verdrehunegswinkel gı 
das Vorreceht, weil die innere elastische Ar- 
beit A mit 9 proportional ist und diese nach dem 
Castiglianoschen Prinzip ein Minimum sein 
muß (Handb. der Phys. Bd. 6, S. 71). Folge- 
dessen ist Herrn Lindners Aufsatz im Gegen- 
teil ja gerade eine Rechtfertieung der Grashof- 


schen Annahme. Von den vielen Versuchen. die 
für die gute PBrauchbarkeit dieser Annahme 
sprechen, möchte ich hier nur an die Versuche 


von K. Böttcher verweisen  (Forschungs- 
arbeiten. Heft 337. 
Budapest, E. Herrmann. 443 
Erwiderung. Nein Beitrax dient nicht dazu. 


die grundlegenden Berechnungen von Grashof zu 
bekämpfen, sondern sie zu ergänzen und zu för- 
dern. Eine unendlich große Spannung, die prak- 
tisch unschädlich bleibt, bedeutet einen Wider- 
spruch, der gelöst werden muß. Bei der Erörte- 
runz von Versuchen, die mit der Theorie nicht 
stimmen, findet ©. v. Bach: es scheine fast. als 
ob zu den Radien jeweils ein kleiner Zuschlag zu 
reehnen sei. Die hyperbolische Form entspricht 
diesem Sinne. 

Wenn. z. B. in der Festiekeitsmaschine, eine 
bestimmte Durchbiegunz verursacht wird, so 
werden sich die Spannungen nach dem Prinzip ver- 
teilen, dab die Summe der inneren Arbeiten ein 
Minimum wird; in diesem Falle aber nicht all 
vemein eilt die hyperbolische Form. Wenn je 
doch. z. B. am Eisenbahnzuzrhaken. eine freie Kraft 
wirkt, wird sie sich nach der „tiefsten Lage“ ein- 
stellen, und die Spannungen werden sich so ver 


teilen, daß die Durchbierung und die Arbeiten 
erößer ausfallen: in solehen Fällen besteht also 


die parabolische Form. 
Übrigens ist die hyperbolische Form nicht eine 

rashof'sche Annahme. 
Stuttzart. 


G. Lindner. 449a 
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